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Introduction 


SCOLE stands for the "Spacecraft Control Laboratory 
Experiment". The objective of the SCOLE Program is to 
provide an example configuration and control objectives 
which enables direct comparison of different techniques in 
modeling, systems identification and control. The "SCOLE 
Design Challenge" was formulated in 1983 by L. W. Taylor and 
A. V. Balakrishnan. The details of this challenge are reprinted 
at the end of this document. 

Annual SCOLE Workshops have been held for specialists 
to share and compare their research results. This proceedings 
is a compilation of the material presented at the 5th 
Workshop held at Hilton Lodge at Lake Arrowhead, California 
on October 31, 1988. 
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ABSTRACT 


A Perfect Model Following (PMF) controller for the linearized SCOLE 
project is described and evaluated. The plant includes both the flexible and 
rigid body modes comprising 8 modes and 16 states. The PMF controller is a 
special case of the output following or Command Generator Tracker concept 
(CGT) which is presented and discussed. Results obtained show that all state 
tracking causes some violation of the control moment constraints. 

However a redesign for the tracking of only selected outputs should be 
feasible. 
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I, INTRODUCTION 

The Spacecraft Control Laboratory Experiment (SCOLE) consists of a 
large antenna attached to the space shuttle orbiter by a flexible beam. The 
purpose of the experiment is to evaluate control laws for flexible 
spacecraft as proposed by Taylor and Balakrishman [1]. 

This report presents the results and the methodology used to evaluate a 
control law by the Perfect Model Following (PMF) method* PMF is a special 
case of the Command Generator Tracker concept (CGT) developed by Broussard 
[2]. In Section II of this report both concepts are outlined. 

The SCOLE system used is a linearized model which consists of flexible 
and rigid body modes. The flexible model state variable vector consists of 


ten states representing 5 second order modes. 

T , 

1 ■ lq^ 4 X q 2 4 2 4jl 

The control is the 8 element vector 

u T = IT , T , T , F , F , T , T , T ] 
sx sy sz rx ry rx ry rz 


Where 


T . 

T . 

T 

sx 

sy 

sz 

T . 

T , 

T 

rx 

ry 

rz 

F rx' 

F 

ry 

* 


F , F m Forces applied at reflector 

rx ry 

The output vector for the system is as follows: 

r T - l*,. V V v V * t - V * t . * t - V i r J 

where 


6 , 0 , Ui * Shuttle attitude, 
s s r s 

. . * 

4 , 0 , - Shuttle attitude rate, 

s s $ 

♦ , ©r# * Reflector attitude. 

... 

, ip *» Reflector attitude rate 
r r r 

r , r - Deflection of tip with respect to the base, 

x y 
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The flexible system, then, is represented as follows: 
x ■ Ax + Bn 


(4) 


y “ Cx (5) 

Where A is a (10 x 10) matrix, B is a (10 x 8) matrix, and C has dimensions 
(14 x 10). Matrices A, B and C are given in the appendix. 

In order to define complete motion we have to add the rigid body 
equations. This rigid body model is represented as follows: 






0 -x y 


F 

RB 

- T + T + 

Z 0 -x 


rx 

F 

RB 

“$ ~r 



ry 

^RB_ 


-y x 0 


0 


($) 


Where I is a (3 x 3) matrix representing the combined inertia motion given 
s 


by 


1.132508 x 10 
7555. 

-115202. 


7555. 

7.007447 x 10* 
-52293 


-115202. 
-52293. 
>.113962 x 10* 


Vector T and X and F , F are the controls from Eq. (2), and x, y, and z 
“s r rx ry 

are 18,7136, -32,4363, and -129,621 respectively. The (3 x 3) matrix in the 
RflS of Eq. (6) corresponds to the moments due to the applied forces. 

The total displacement Yjq> then, is expressed as the displacement of 
the flexible model plus the rigid body displacement. In equation form we 
have 

7jj) = y + RBV (6a) 

where y is from Eq. (5) and RBV (Rigid Body Vector) is given by 

(RBV) * ^bb'°RB' ^RB' W e RB* *RB' °* °* ^RB’ 0 RB'^RB'^RB' ®RB' 

In order to couple the two models it is necessary to express Eq. (6) in 

state variable form, and then form an augmented system with Eqs. (4) and 

(5). Toward this end we define the rigid body state variables as follows: 
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**1*18 * ^ X 2*RB " ^RB 
* X 2 J RB * ^RB 
**3*RB " (x 4 ) RB “ °RB 
* X 4*RB “ ®RB 
* X 5*RB " * X 6*RB * ^RB 

(x 6*RB “ ^RB 

Combining Eq. (6) and (7) reanlta in the following system: 

X RB “ ^ X RB + B RB U 
y RB * C RB X RB 

Whore is a (6 x <) matrix, is (6 x 8) and is (14 x 6). 

matrices are given below: 


(7) 


( 8 ) 

(9) 

These 


0 


1 0 0 0 0 


0 0 



0 


0 

0 

0 

0 


0 

0 

0 

0 

0 


0 0 
1 0 
0 0 
0 0 
0 0 


0 

0 

0 

1 

0 


( 10 ) 


4/J 



0 

0 

1 

0 

1 

0 

8.844xl0~ 7 

^ -8.46xl0 _1 ° 

1 

1.43 lxl 0*" 8 

1 

5.741x10~ 7 

0 

0 

1 

0 

1 

0 

-8.467x1O" 10 

-7 

1.427x10 

1 

1.035xl0~ 9 

1 

-1.846x1O" 5 

0 

0 

1 

0 

1 

0 

1.431xl0~ 8 

1.035xl0~ 9 

1 

1.408xl0' 7 

1 

4.433xl0 -6 

0 ' 

0 

1 

0 

1 

0 

-4 1 

1.149x10 1 

8.84xl0 -7 

1 

-8.467xl0~ 10 1 

1.431xl0~ 8 

o 1 

0 

1 

0 

1 

0 

-8 1 

-9.037x10 ° 1 

8.467xl0 -10 

1 

1,427xl0~ 7 

1 

1.035xl0 -9 

0 ^ 

0 

I 

0 

1 

0 

4.491x10~^ ^ 

1.431xl0 -8 

1 

1.035xl0~ 8 

1 

1.408x10~ 7 


( 11 ) 


4 12 





1 0 0 0 0 0 
0 0 1 0 0 0 
0 0 0 0 1 0 
0 1 0 0 0 0 


0 0 0 1 0 
ooooo 
ooooo 


c 


RB 


0 


1 

0 

0 


0 

0 

0 

0 


0 o 

0 0 

1 0 

0 0 


0 

0 

0 

1 


0 

1 

0 

0 

0 

0 

0 


0 1 0 0 0 0 

oooioo 

0 0 0 0 0 1 


( 12 ) 


The total system, which we call the plant is then given by: 


x = A x + Bu 
p P P P P 

y = C x 
*P P P 


There the augmented 

Fa 


matrices 


» 

i 


0 



B and C are given by 
P P 


A =■ 
P 



^ J 


(16x16) 


(13) 

(14) 


(15) 
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B 


B 


B, 


RB 

[C 


(16x8) 

C RB ] (14x16) 


(16) 


(17) 


x is the (16x1) state vector given by 

4 - ! *bb' 

Where x is the flexible model state vector of Eq, (1) and x D _is the rigid 

KB 

body state vector of Eq. (7). 

Up is the plant control vector defined by Eq. (2), i.e. 


Finally the adaiissible controls must satisfy the inequality constraints 
|F| < 800 lb (18) 

|T| < 10.000 ft-lb (19) 


Where F represents the forces applied at Reflector (F and F ), and T the 

rx ry 

moments applied at Shuttle and Reflector (T , T , T and T , T , T 

sx sy sx rx ry rx 

respectively). 

Our objective will be to find u (t). within limits, in order to force 

P 

to decay from 20°(0.349 rad) to 0° in about 10 seoonds. 
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II. METHOD 


A. Command Generator Tracker (CGT). 

In this section, we will outline tie basic elements of the CGT and 
state the conditions for perfect model following. 

Given the continuous linear system 


i (t) * A x (t) + B u (t) < 20 > 

P p P P P 

y (t) - C x (t) < 21 > 

P P P 

where x (t) is the (nxl) plant state vector, n (t) is the (mxl) control 
P P 

vector, y p (t) is the (qxl) plant output vector, and A p , are matrices with 
the appropriate dimensions. Our objective is to find the control u p (t) such 
that the plant output vector y p (t) approximates reasonably well the output 
of the following model: 


x (t) = A x (t) + B u (t) < 22 ' 

m m ju pi hi 

y (t) - C x (t) < 23 > 

m mm 

Vhere x (t) is the (n xl) model state vector, u (t) is the (mxl) model 
mm “ 

input or command, y^t) is the (qxl) model output vector, urd A^. B^ are 

matrices with the appropriate dimensions. It is assumed that: the pair (A^, 

E ) is cortrollable and output stabi1 izable, the pair (A , C ) is 
P * v 

observable, and B has full rank. 

P 

When perfect output traoking occurs (i.e. when y^ ■ y m for t 1 0) we 

dc-fire the oorrespcrdix g r-tete and control trajectories to be the ideal 

state x*. and ideal control u*. trajectories, respectively. By definition, 
P P 

the ideal plant satisfies the seme dynamics as the real plant. Also, the 
output of the ideal plant is defined to be identically equal to the model 


output. Mathematically, 


i*(t) * A x* + B u*, for all t l 0 
P P P P P 

y*(t) = y (t) “> C x * <L x_ 
p m p p 


m 


(24) 

(25) 
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The basic assumption of the CGT theory is that the ideal trajectories 
are linear functions of the model state and model input, i.e. 






x_(t) 


In Eq. (26) we have restricted v to be a constant inpnt; otherwise 
derivatives of the model input may be required. Matrices S^, i» j ■ 1. 2 
are found as follows. Concatenate the ideal plant state equation with the 
ideal plant output equation described by equations (24) and (23) and get 


(27a) 


Substitute Eq. (26) into Eq. (27a) and obtain 


(27b) 


Now we differentiate the first equation in (26) to obtain 

i = S.. x + S u 
p 11 m 21 m 

and since u is a constant input, we have 
m 

* 

* n • 

x = S,. x 
p 11 m 


(27c) 


Now we substitute the equation for the model dynamics into (27c) to obtain 


x =S..A x + S., B u 
p 11 m m 11 m m 

We combine (27d) with (23) to obtain 


s n A 
11 m 


B„ B 
11 m 


(27d) 


(27e) 


Mow equate the right-hand sides of (27b) and (27e) and get 


4/6 


iimiii in mi un ,i mb ill ill i 



S 1t A 
11 m 


S„B 
11 m 


B 


11 


21 


12 


22 


and noting that z and a are arbitrary we obtain 
iii n 


Si, A 
11 m 


S„B 
11 a 


B 


L C P 


'11 


'21 


'12 


22 


(27f) 


If we define 


vH 

rH 

a 

1_ 

°12 

. 

A 

P 

B 

P 

f-i 

<s 

a 

_1 

°22_ 


C 

L p 

- 1 

o 


-1 


(27g) 


then the equations to be solved are 


s_ 

» ft.* 


A + 

C 

(28a) 

11 

11 

11 

m 

12 m 


S 12 

-°u 

s n 

B 

m 

l 

(28b) 

s_ 



A + 

C 

(28c) 

21 

21 

11 

ID 

22 ra 


S 22 

-°21 

s u 

B 

ID 


(28d) 

The 

existence of 

the inverse 

requires that the number of controls, m. 


equals the number of outputs, q. If m > q a pseudo - inverse may be required, 
while the case m < q might not have a solution. 

For the SCOLE model, as seen in Section I, m < q. However, as will be 
seen later, this requirement is not necessary if the so-called conditions of 
perfect model following (PMF) can be satisfied. 

Going back to the CGT controller, Bq. (26), the ideal control is given 
by 

% - s n x . * s 22 W) 


4)7 





The existence of S.. is assured provided (1) u is a constant, (2) 

lj » 

m > a, and (3) the product of the i-th eigenvalue of and the j-th 

eigenvalue of A does not equal unity for all i, j [2]. 

D 

• * 

If at some time x ■ x then u * u will assure perfect tracking, 

p p p p 

However, if y t y at t * 0, we may achieve asymptotic tracking provided a 
p m 

stabilizing output feedback gain is included in tbe control law. To see 

this, consider the error equation 

e a x — x * A x + B u “A x — B u * A e + B (u — u) 
p p pp pp pp PP P PP P 

Now, with an output stabilizing gain, the control law is 

u* + * (y - y ) 

P P m p 

= u + K C e 
P P 


then the error equation becomes 

6 = (A - B K C )e 
P P P 


(30) 

(30a) 

(30b) 


and the error will approach zero provided that K is an output stabilizing 
feedback gain. 

Note that the time for y to approach y depends upon the eigenvalues 

P m 

of the matrix (A - B K C ). 

P P P 

B. Perfect Model Following (PMF) 

The perfect model following conditions are a special case of the CGT 
when (1) the state vector is available, and (2) it is assumed that x*(t) * 


x (t), i.e. C - C 

m P n 


i. 


Therefore from (26) the PMF conditions imply that 


s n " 1 


(31) 


S 12-° 

The ideal plant input, given by Eq. (26), is, as in the CGT controller, 

u ■ S„, x + S-„ u 
p 21 m 22 m 
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even if the number of output*, q, is larger than the number of control*, m. 

We should point out, however, that the PMF has certain disadvantages 
which can make the design more difficult, namely 

(i) We assume that the state vector is available. This is not a 


realistic assumption. 

(ii) In a PMF controller all the plant states are forced to follow all 

the model states by definition. This fact shows that the dimension 

of A is equal to the dimension of A . If dim A is large, like in 
m P P 

the SCOLE project, computation of and b® difficult* 

If at some time x « x , then u =* u will assure perfect tracking* 

p p P P 

However if x £ x at t » 0 a state feedback stabilizing gain is required to 
p m 

achieve asymptotic tracking* Using a method similar to the one that led to 
Eq. (30b) it can be proved that the error equation for the PMF case is 


o = (A - B K)e 
P P 


4/9 


and the control law will be of the form 


u + K (* - x ) 

p n p 


where lit t 


stabilizing gain. Note that in this case the 


• * 

error is given by « » x -x - x - x because by definition x ■ x in the 

P P ■ P p n 

PMF theory. 

A block diagram of the system (Plant. Model, and PMF) is shown in 
Fignre 1. 


420 

















III. DESIGN PROCEDURE 

A. Approach: PMF vs. CGT. 

As a preliminary controller design for the linearized SCOLE project we 
designed a PMF controller so that all states of the process follow all the 
states of the reference model. In this manner all 14 outputs defined in (3) 
should have acceptable responses. 

For a CGT controller it is required that n ^ q if equation (28a) to 
(28d) have a solution. It should be noted, however, that it is possible to 
design a CGT controller if appropriate changes are made in the dimension of 
the output vector or in the number of plant outputs that are to follow an 
equal number of model outputs. These considerations are discussed in Section 

IV. 

Here we present the design criterion of a PMF controller. Specifically, 
in the subsections to follow, we present the desigr considerations and 
method of evaluation of matrices A m' S 21* S 22* and K. 

B. Plant: 

The plant for the SCOLE problem is given by 

x - A x + B u (34a) 

P P P P P 

y =* C x (34b) 

P P P 

where A is (16 x 16), B is <16 x 8) and C is "(14 x 16) . All these 
p P P 

matrices are given in equations (15), (16) and (17). The controls u are 

P 

constrained by (18) and (19), and the objective of the problem is to have 

4__ >0 in about 10 seconds. 

Kd 

C. Model, 

The model is given by 

i - A x + B u (35a) 

m mm mm 

y » C x (35b) 

m mm 
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According to the theory developed it the lust section for a PMF 

controller design we should choose A and B so that the system 

9 mm 

specifications and Eq* (32) are satisifed. Thus dim A^ « dim A^, and dim B^ 

» dim B . 

P 

For the present design matrix B was selected to be equal to B * It 

m p 

should be noted that the set of equations (32) have a solution when the 

column vectors of the difference matrix (A — A ) and of the matrix B are 

m p » 

linearly dependent on the column vectors of the matrix By choosing 

B = B this condition is satisfied for Eq. (32b). For Eq. (32a) to have a 

m p 

solution A and A should have their odd numbered rows equal, 
m p 

The design criteria for A^ was to have 0 in about 10 seconds. 

Similar to the plant, the model will have 8 modes (16 states). Modes 6, 7, 
and 8 correspond to 0^ and \p ^ respectively. The first five modes are 

the flexible modes of the system. The model matrix is given by 



The A ., represent second order systems of the form 
ml' $ 


A A 


mi 



-2 


<d . 
ni 


4 23 



where £ is the damping ratio, and is the natural frequency of the i-th 
mode. In order to keep the overshoot at a reasonable level we will fix the 
damping ratio to £*.707 . Therefore, the design criteria will be the 


selection of the w 


ni' 


u» will affect both the 2% settling time of each mode as well as the 
n 

control law a (t) through matrix S-, (see Eqs. (32a) and (29)). The effect 

p LI 

of the to (and thus of A ) on the control u (t) was tested by selecting 

ni s m p 

appropriate natural frequencies (by the design procedure stated below) for 

A , running computer simulations, and observing the values obtained for the 
m 

u^'s. Using this empirical method we found out that the larger the values 

of the cd . . the larger become the u , • Therefore we selected the smallest 
m's p's 

natural frequencies that would satisfy the specifications and the 

constraints on u (t). 

P 

The 2 % settling time criterion was applied as follows: For a second 

order system the 2% settling time, T , is given by 

s 


s £w 

n 

Each mode of matrix A will exhibit a settling time given by the last 

m 

equation. Now, our objective in this report is to have the sixth mode (<l D _) 

ho 

decay in about 10 seconds. Therefore, T^ i 10 sec* for this mode. Then, 
with * 0.707 the frequency of this mode should satisfy 


or. 


<d 2. 0.567 rad/sec. (36) 

n 

The other 7 natural frequencies were selected as espleired before .? i* 

order to satisfy, if possible, the constraints on the control. 

The final A is shown in the appendix, 
m 
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C. PMF Design 


For a PMF we should find matrices S., and S_* which should satisfy 


< A m " V " B p S 21 

B p S 22 = B m 


C.l. Evaluation of S 


21 


Define 


AA = A - A 


Then, Eq. (37a) becomes 


(37a) 

(37b) 


(38) 


B S * AA (39) 

P 

where AA is the (16 x 16) matrix given by Eq. (38). In Eq. (39) we have 

dropped the subscript of matrix S 21 for clarity. Note also that matrices 

and S have dimensions (16 x 8) and (8 x 16) respectively. 

To solve for S in Eq. (39) we will make use of the structures of 

matrices B and AA. Both of these matrices in the SCOLE project have theix 
P 

odd rows equal to zero. Therefore Eq. (39) becomes 


B S 
P 


0 


0 

BS 2 


aa 2 

0 


0 

BS. 


AA 

4 


4 

0 


0 

• 


• 

• 


e 

• 


9 

« 


9 

0 


0 

BS 1* 


-1 

VO 

H 

3 

_ 1 


* AA 


(40) 
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where BS,., BSBS - £ are the even rows of the product B S, and AA_, 
2 4 lo p 2 

AA^ f . • . the even rows of AA. Eq. (40) can be written as 


BS„ 


AA„ 

2 


2 

BS 


AA. 

4 


4 

♦ 


• 

a 


• 

• 


a 

• 


♦ 

BS, ^ 


AA., 

16 


16 


(41) 


The left-hand side of Eq. (41) is the product of a reduced F^ matrix 

times the „ matrix. This reduced B matrix is the original B matrix 
21 p P 

without its odd rows. We will call it Then, 


B, 


B PR 


16 


(42) 


(8x8) 


where the B., , i = 2, 4,....,l6 are the even rows of B . 

l 9 s p 

The right-hand side of EQ. (41) is the corresponding reduced AA matrix. 


We call it AA^. Therefore Eq. (41) can be written as 

b pr 8 = aa r 


(43) 


Note that the dimensions of Bp^, S, and AA^ are (8x8), (8 x 16) and (8 x 
16) respectively. For the SCOLE project natriy Epp 
therefore equation (43) has a unique solution given by 
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is non-singular; 






( 44 ) 


S 21 * “PR “r 
I n snmmary we have 

(i) Matrix is given by Eq. (44) and is unique, 

ii) Matrix S-. depends on matrix A . 

21 *n 

The computation of Eq. (44) was done by a computer program, and S 21 is 
presented in the appendix. 

C.2. Evaluation of S 22 

^22 * 8 solution of Eq. (37b), i.e. 

B p S 22 " B m 

Because in our model we selected B = B , we have 

p m 


B S = B 

p 22 p 


(45) 


The dimensions of B end S„„ are (16 x 8) and (8 x 8) respectively. Using 

P 22 

the method developed to evaluate S^ in section Cl, we can express Eq. (45) 
as a set of 64 equations in 64 unknowns (the elements of S 22 } * The solution 
of this set of equations is unique and, therefore, matrix S ^ is the (8 x 8) 

identity matrix, 1,e * S 22 = *8* 

D. State Feedback Stabilizing Gain, 

A feedback stabilizing gain matrix, K, was designed using the pole 

assignment method. 

The open-loop plant is given by 

x *= A x + B u 
P P P P P 

and the feedback systems becomes 

x (A - B K)x 
P P P P 

when K is the (8 x 16) feedback gain matrix. 

The poles of the matrix (A - E F) were located to the left of the 

P P 

poles of the model in order to achieve fast tracking in the event that the 


4 27 


initial conditions of the plant and model states are different, i.e. z £ z 

pa 

at t = 0. 

The actual design of the time invariant matrix K was achieved by means 
of a computer program using the method given by Brogan [3] . 

Matrix K is given in the Appendix. 

E. Computation 

To simplify computations the plant state, x , the model state x , and 

p m 

ti e conrerd i* ore concatenated into one system of 40 "new" states; the 

ti 

first 16 states are the plant states, the next 16 states are the sodel 

states, and the last 8 states are the command, u . 

m 

The commands are unit step functions and are modeled as follows: 

u = 0.0, u (0) - 1.0 
m m 

Therefore the augmented system becomes 


x = A x + B u 
P P P P 

x “A x + B u 
m m m m 

u = 0.0 

m 


(46a) 

(46b) 


(46c) 


u = S 01 i + u + K (x - x ) 
p 21 m m m p 

Substituting (46d) into 46a) gives 


(46d) 


x « (A - BK)x + B(S + K)x + Bu 
p p p mm 

where B = B = B 
p m 

The concatenated system then becomes: 


(46e) 


x = Ax 


where 


m (40x1) 
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and 


(A - BK) 
P 

0 

0 


B(S+K) 

A 

a 

0 


(40x40) 


A is a (40 x 40) matrix. 

The state and model output are given by 
y - C x 


y = C x 
'm o 


where C 


P 

T 


C « C and is given by equation (17). x and x are: 
pm pm 

[x 1 x 2 ....x 16 J 
tx 17 x 18*'“ x 32^ 


where the arc given by Eq. (47) e The rigid body states are: 


T 

(x HB > p “ lx ll X 12- X 16 ] 

T 

^ x RB^m “ tx 27 X 28- X 32* 

The total displacement becomes 


y TD 


y + y 


RB 


(49) 


where y is the orginal output, and y^g is the rigid body output. 

The actual computation of Eqs. (47), (50a) and (50b) was carried 


(50a) 

(50b) 

(51a) 

(51b) 

(52a) 

(52b) 

(53) 

out by 


a computer program developed by the authors using an integration subroutine. 
The integration step chosen is 0.01 seconds because the absolute value of 
the largest eigenvalue of matrix A, Eq. (49), is less than 10. 
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IV. RESULTS 


The results of a simulation of the system are shown in this part of the 

report. The initial conditons and the command values used are the following: 

Plant I.C's : x .(0) - 10.00, i«l,2,...,10 
Pi 

Model I.C's : x .(0) - 10.50. i-l,2,...,10 

mi 

Rigid Body Plant and Model: 

*RB (0) " °- 349 rad (2 ° 0) 

♦rb <0) - 0 - 0 
9 m (0) - 0.0 

W 0) ■ »•» 

V 0) - 0.0 
V»> - 0-0 

Command : u ,(t) * 0.0, t 2 0, i«l,2,...,8 
ml 

Several other runs with different sets of initial conditions were made 
which are not presented in this section. However, some comparisons with the 
results presented here are made. 

Here we present four sets of plots (1) output plant and model of 
flexible modes, (2) output plant and model of rigid body, (3) total 
displacement, and (4) controls. 

Figure 2: The eight plots presented here correspond to the original 

(flexible system) output vector of Eq. (3), not including the reflector and 

shuttle attiti'de rates. These plots represent, then: shuttle attitude 

6 , ip ), deflection of tip w.r.t. base (£ , £ ), and the reflector attitude 
s s x y 

]p^) . In each plot we show the behavior of the plant and model 

together. 

It can be seen that because the initial conditions of the plant and 

model states are different (x £ x at t * 0) the plant outputs achieve 

P m 
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asymptotic tracking. At about t “ 10 seconds x^ approaches x^, and the plant 

outputs achieve perfect tracking. Other runs (not presented here) with the 

same initial conditions show that perfect tracking occurs for all t 2 0, as 

expected. Therefore, in both cases PMF is accoaqplished. 

It should be noted that it is possible to decrease the time for y p to 

approach y by making the model natural frequencies larger. However, this 
m 

consideration must be taken carefully because a faster system requires more 

control effort. This, and other considerations of the model design were 

treated in Section III of this report* 

Figure 3: This figure shows plots of the rigid-body attitude, 4gg. 

From this plot we see that our objective was achieved, namely, to have 

A + 0 in about 10 seconds. A computer print-out (not included here) shows 
RB 

that in about 10 seconds ♦gg decays to a value less than 2% of its initial 
condition value (*gg(10) - -0.675xl0 -2 ) . Also, it can.be seen that the 
plant approaches the model assymptocially, as expected. 

The other rigid body modes (0^ and if^) are not presented here because 
they stay at zero level for all t^O. This is expected because their 

initial conditions were set to zero. 

Figure 4: In this figure the plots of some elements of the total 
displacement vector are presented. The total displacement vector is given 

by Eq. (6a), i.e. 

y = y + RBV 
*TD 7 

where y is the original output vector (flexible modes) and RBV the rigid 
body vector. The sir plots presented here correspond to the following 
elements of the above vector: 

y TDl “ *s + *RB 
y TD2 " 6 + °RB 
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7 TD3 = ^s + ^RB 
y TD9 = *r + ♦RB 
7 TD10“ ®r + ®RB 
7 TDll“ + ^RB 

The first 3 are the shuttle attitude plus rigid body, and the last 3 
the reflector attitude plus rigid bo&f. 

These plots show that the total displacement approaches zero in a short 
period of time; in all cases in about ten seconds. The following table shows 
the time, in seconds, to approach zero for each output presented in Fig. 3. 



t TDl 

y TD2 

y TD3 

y TD9 

y TD10 

y TDll 

Time 

5.5 

7.0 

7.0 

10 

7.5 

5.0 


This data was taken from a computer print-out of the simulation. 

The plots of the rate elements of the vector are not shown here . 

« * a • 

These are: rates of shuttle attitude plus rigid body (♦ + i , 0 + 0 , 

s KB s RB 

• • * • • . 
ip + ip DD ) * and rates of reflector plus rigid body (♦ + i + 

s kb r RB r RB r r 

9 

• The behavior of these elements are similar to the ones presented here. 
The total displacement elements y^^ and y TD8 are t ^ ie ^ e ^^ ect ^ on °* tip 
w.r.t. base (£ ^ and £ )• Because the corresponding elements of the rigid 


body vector are zero (see Eq. (6b)), then y, 


TD7 


y 7 8nd y TD8 


y g where y ? 


and y g are the original flexible modes presented in Fig. 2. 


Figure 5: This figure shows the plots of the controls, v ft). It is 

P 

clear that not all the controls are within the permissible limits given by 

Eqs. (18) and (19). The worst case corresponds to the first control, T , 

sr 

which at t = 0 attains a value of -186,000 ft-lb; as tir»e passes T 

si 
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decreases very fast, and at about 7 seconds its valc-e is less than the 
10,000 ft-lb limit. For t 2 7 seconds its valne decreases steadily. 


The second and third controls (T jy and T sj ) at t = 0 have values of 
15,000 ft-lb and 27,000 ft-lb respectively and for t 2 1 sec regain within 
limits. 

The sixth control (T ) at t = 0 has a valne of 10,900 ft-lb. For 

t 2 0.5 it remains within limits. Control numbers 4, 5, 7 and 8 (F^, F ry » 

T and T ) remain within limits for all t 2 0. 
rz rz 

As we see the first 3 controls present very large deviations fror' tie 

minimum values desired. A solution to this problem would be a new design of 

matrices S.. and K as seen from the control equation given by 

u = S„, x + K (x - x ) + u 

p 21 m m p m 


(54) 


(55) 


As we recall, S 21 is the solution to 
(A - A ) = B S 

m p p 21 

and it was shown, due to the structures of matrices (A^ — A^) and B^, that 

this solution is unique. Once matrix A^ is selected, S 21 is unique. Thus, 

to redesign S.,, requires a different A . Fc-wever A controls the settling 
9 21 m m 

times of the states (see Section III of this report). Therefore, we see the 
difficulty of reconciling the limits on the controls and settling time (our 
objective with respect to • 

In relation to redesigning matrix I, the state feedback stabilizing 

gain, we see from Eq. (54) that when the initial conditions of the model and 

plant are equal, i.e. x (0) = x (0), then the second term it Fq. (54) is 

p m 

zero and K has no effect on o^. To see the effect of K on u^ when x^(0) 4 1 

x (0) we did a simulation with x (0) = x (0) and the result showed a 
m pm 

decrease in the values of u , but not enough to meet the constraints. For 

P 

example, the first control # at t = C, reached a value of u 
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— 130,000 ft-lb for the case x (0) = x (0). This value is much better than 

p m 

the value obtained when x (0) = 10.0 and x (0) = 10.5; however, in loth 

p m 

cases is ^ ar ftom being within the limits <, 10,000 ft-lb. 

Therefore, matrix K has little effect on u^ in general. 

Another approach tried was to irclide v saturation element in the 

system. This element would limit u to its constrained value whenever u is 

P P 

off limits. If we call this saturation value u , we see that 

P 

u S £ u = S_, x +K(x - x ) + u 

p p 21 m m p m 

and the plant and model systems will be 

• s t S u S 

x = A x + B u 
P P P P P 

i = A x + B u 
m mm mm 

$ s s 

where x E plant value when u = u . It is clear from above that x £ x 
P P P P P 

and, because x is independent of v , it follows that the plant will not 
m p 

follow the model, thus destroying the purpose of the PKF controller. / 

simulation was done with the saturation element to see if at some time t > C 

the plant would approach the model, and the results in fact showed 

divergence. Therefore, this approach was disregarded. 

As seen from the plots all u . 0 as t -*■ ®. This should be the case 

pi 

because here we have u . (t) ® 0. A simulation with the commands equal to the 

mi 

unit step (u . (t) * 1 (t)) showed again that x + x . We conclude, then, 
mi p m 

that because of the constraints on u it is very difficult, if not 

P 

impossible, to design a PMF controller that will satisfy both the limits of 


Up for all t and the objectives of the system. As explained in the next 
section, we believe that a CGT controller would be able to satisfy both. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

A PHF controller design for the SCOLE project was presented l-ese. T J-« 
objective was to have 0 in about ten seconds, and to maintain tic 

controls. jip(t) within specified limits (given by Eq*. (18) and (19)). 

The first of the objectives was net as we can see from the results. 
The rigid body outputs, as well as the flexible a ; ode outputs decay to zero 
in less than ten seconds. 

However, it was not poss ible t o r eet el 3 tie constraint* on the 

controls. The first three controls, the moments applied to the shuttle (T^, 

T , T ), reach unacceptable values (-186,000, 15,000, and 27,000 ft-lb 
sy sz 

respectively) at t * 0. They decay, however, to values within the limits in 

7 sec, 1.0 sec, and 1.0 sec respectively. 

It was explained in the last section of this report the reasons behind 
the behavior of these choc’s: an attempt to improve the controls would 
require a slower model, raking the response of the plant also slower (the 
model following action). A slower system would not meet the first objective. 
The main problem encountered in the PMF design, as explained, was the fact 
that matrix S, 1 is unique, and thus, we do not have too much freedom to 
select model parameters to meet both objectives* 

This preliminary design, Icvever, vas attempted to prove that it is 
possible to design a direct model following controller for the linearized 
SCOLE project. 

A more realistic approach is to design a CGT controller using tie 

method presented in Section II. Unlike the PMF, the CGT does not require die 

A = dim A . This is an advantage because we have the freedom to choose 
m p 

selected plant outputs to follow an equal number of model outputs. In a PMF 
we are forced to have all the plant states to follow the model states. 
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However; a solution to Eq. (28) of Section II requires the nunber of 
controls, m, to be equal to the number of outputs, q; or at least m > q. In 
this project we have m > q (m = 8 and q ■ 14), Tt could overcome this 
problem by re-defining the output vector of Eq. (3). Two recommendations 
could be explored. 

(i) To have the elements of y to be the sum of position and rate as 
follows: 


+ V e . + V K + V v y V V V V S + % ] 


(8x1) 


A corresponding matrix should be defined accordingly. In this 
case we would have m - q = 8. 

(ii)Ve could have fewer plant outputs follow an equal number of mcdt.1 
outputs. In this case we have m > q. The selection of the 
outputs that are to follow the model outputs will depend on the 
objective of the design. In the SC01F project *r agfroprJate 
selection would be the inclusion of in the output vector in 

order to meet the objective. Therefore the dimension of y might be 
as low as 1. 

Because the CGT is an output foil over # {i d because we do not assume 
that the state vector (or part of it) is available for feed-back, an output 
stabilizing gain should be designed to make the error of Eq. (30a) decay to 
zero in case the initial conditions of plant and model are different. 

In conclusion, this preliminary approach has shown the feasibility of 
designing a model following type of controller for the SCOLE project. A more 
suitable controller should be possible using the CGT method. 
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1. INTRODUCTION 


This report summarizes participation in the NASA LaRC Spacecraft Control Lab¬ 
oratory Experiment (SCOLE) under NASA Grant #NAG-l-720 in the study of 
various approaches in analysis and control of systems having characteristics of Large 
Space Structures (LSS). The work is in continuation of an ongoing project originated 
under Contract NCC1-108 entitled “Development of Decentralized Control Algo¬ 
rithms for Flexible Structures”. Principal investigators for this continuing effort 
have been Professors Steve Yurkovich and Umit Ozgiiner, and Prof. K. A. Ossman 
contributed to the effort as a Research Associate. Several graduate student research 
assistants, namely P. Cooper, A. Tzes, H. Fu, and M. Cheung, have assisted in the 
research program. M. Cheung, F. Khorrami and A. Iftar contributed to controller 
designs in COFS work. 

The primary objective of this research has been twofold: first, to assess the perfor¬ 
mance of various control methods in simulation tests and actual experimentation, 
and second, to determine the effectiveness of various techniques for on-line and off¬ 
line system identification relative to control applications on the SCOLE apparatus. 
Justification for the objectives of this work lies in the inherent need for reliable and 
robust control designs in active vibration control, and in the necessity for accurate 
system models in most large space structure applications. Throughout our work 
for this program and others, we have recognized that traditional approaches, such 
as finite element techniques a/one, may not in some instances be accurate enough 
for use in designing a vibration control system. For these reasons we have been 
pursuing approaches which realize a model directly from experimental data , consid¬ 
ering techniques for system identification and parameter estimation. Despite recent 
efforts along these lines, little ground has been broken toward solution of the many 
problems posed, and the recognition of the immediate need for such research is 
apparent. 

Indeed, these points are especially important for application of adaptive control 
techniques, a major emphasis in this work. Large space structures are excellent 
candidates for adaptive control because the modal data, in many cases, is not suf¬ 
ficiently accurate to design a fixed controller for the structure. An exact model is 
difficult to synthesize because the low resonant frequencies are closely packed, the 
damping is light and often unknown, the model parameters will change when the 
system is put into orbit, and the forces such as those due to gravity and temperature 
gradients will change continually throughout the orbit. Also, the system parameters 
will vary during a large angle slew. With the highly complex mechanical structures 
of the future and under current study, there is a need for robust adaptive controllers 
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whi '!i can meet the control objectives in spite of model uncertainty, nonlinearities, 
and spillover due to residual system interaction. 

The overall group effort at Ohio State within the Department of Electrical Engineer¬ 
ing for work at LaRC is summarized in Figure 1.1. While the focus of this report 
is on work related to SCOLE, we have begun preliminary studies on the MiniMast 
configuration at LaRC; details of these studies will appear in a later report for con¬ 
tinuation of this project. In Figure 1.1 the dashed lines represent work in progress 
or projected studies, and the various acronyms which appear here and throughout 
this report represent, respectively, Variable Structure (VS) Control, Self-Tuning 
Adaptive Control (STAC), Maximum Entropy/Optimal Projection (MEOP) Con¬ 
trol Design, and Model-Reference Adaptive Control (MRAC). Not represented in 
the figure, but also forming a portion of the research in this effort has been anal¬ 
ysis and control design for the envisioned flight article in the Control of Flexible 
Structures (COFS) program. Because of the direction that program has recently 
taken, we have used this report as a vehicle for describing our findings and controller 
designs for the model of the COFS mast, an we will utilize the expertise gained on 
other configurations that may be proposed as a replacement? 

The remainder of this report is organized as follows. In the next section we briefly 
formulate and review the basic problem as related to SCOLE. Following this is a de¬ 
scription of the various controller design approaches taken in this research. We have 
grouped these design procedures into three main categories: Nonadaptive Designs, 
Direct Adaptive Designs, and Indirect Adaptive Designs. In each case extensive 
simulation studies have been undertaken in preparation for experimentation. After 
this is a section describing the actual experimentation completed at LaRC within 
the last twelve months. Following this is a brief summary of the work we completed 
relative to the COFS mast model. Finally, we conclude with comments on future 
directions in the research and in the laboratory experiments relative to SCOLE and 
the MiniMast. 


2. SCOLE APPARATUS 


The SCOLE is comprised of three basic structures, 83.8 by 54.0 inches shuttle 
platform, a 120 inches long mast and a hexagonal reflector panel. There are 15 
sensors and 12 actuators distributed on the structure for the experiment, and our 
controller design studies have focused on different combinations of these actuators 
and sensors, for vibration suppresion in either the mast portion or the reflector 
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Figure 1.1: Summary of LaRC Effort 
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panel. The sensors include nine accelerometers and two, 3-axis rotational rate 
sensing units. Actuators include seven reaction wheels, two gas jets thrusters and 
three Control Moment Gyros (CMG), although the algorithms we currently consider 
do not employ the CMG’s due to their availability status in the laboratory. 

The model we employ in all of our studies for the SCOLE configuration is a ten 
mode finite element model with modal displacement equation 


d q 

dt q 


0 I 
-n 2 D 



( 2 . 1 ) 


where 



1 


-2(u>i 1 

n 2 = 

U)\ 

D - 

N 

3 

CM 

1 


Wo . 


— 2(u>io J 


As noted above, inputs and outputs are selected as needed for various control algo¬ 
rithms and control objectives. Finally, the ten modal frequencies for this baseline 
model are given below in Table 2.1. 


3. NONADAPTIVE CONTROL APPROACHES 


Three fixed (nonadaptive) control schemes have been investigated for implementa¬ 
tion on the SCOLE apparatus. The first two, simple rate feedback and an observer 
with pole placement design are relatively simple in nature, and we omit details of 
those algorithms here. In fact, our primary intention in considering these straight¬ 
forward controller design schemes was to form a baseline comparison for more so¬ 
phisticated designs to follow. In the section later in this report on experimentation 
we discuss the results of the rate feedback scheme and supply a few more details 
there. 

The third nonadaptive scheme we have considered is the Maximum Entropy/Optimal 
Projection (MEOP) approach to designing optimal low-order controllers for high- 
order systems with parameter uncertainties, developed recently by Hyland and 
Bernstein [1], [2]. This technique presents a unique direct method to the design of 
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Table 2.1: Modal Frequencies 


Frequency 

Mode Type 

Axis 

.055 

Pendulum 

y axis 

.088 

Pendulum 

x axis 

.165 

Pendulum 

y axis 

.251 

Pendulum 

x axis 

.254 

Pendulum 

y axis 

.565 

Bending 


.638 

Bending 


1.51 

Torsional 

z axis 

2.94 

Bending 


4.38 

Bending 



optimal, robust, reduced-order compensators in that the compensator design and or¬ 
der reduction are performed simultaneously. This is in contrast to the conventional 
two-step approach to reduced order compensator design in which a compensator 
order reduction is preceded by the LQG design or a LQG design is preceded by 
plant order reduction. 

The fundamental idea of the MEOP approach is to use stochastic modeling for the 
uncertainty in order to improve the system robustness with respect to parameter 
variations and higher order unmodeled dynamics. In [3] a study was conducted com¬ 
paring this design philosophy with other competing asymptotic LQG-type design 
methods, where, for an example representative of lightly damped flexible structures, 
MEOP design fared favorably with regard to stability robustness to uncertainty in 
modal frequencies. The ME design portion for robustification is coupled with an 
oblique projection method to reduce the size of a compensator in some optimal fash¬ 
ion. Such a scheme is particularly appealing in the area of large flexible structure 
control where the size and robustness of the compensator are critical. Thus, the 
MEOP design methodology generalizes the highly sensitive classical LQG design to 
a robust, low-order dynamic compensator design. 


3.1 MEOP Overview 

As a brief introduction to the MEOP approach, consider an n th order linear time 
invariant system 
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x = Ax + Bu + wi , 


y = Cx + «> 2 , 


( 3 . 1 ) 


under the usual assumptions of controllability and observability, where uq, u >2 are 
zero mean white noise processes with noise intensity matrices Fj and V 2 , respectively. 

It is required to design an order robust, zero set-point compensator 


x c = A c x + Fy , u = -Kx e (3.2) 

to minimize the cost functional 
1 t T 

J = Km E J o (x T (t)Rjx(t) + u T {t)R 2 u(t))dt , (3.3) 


where x G R n , x c G u G 7£ m , y G 'll 1 , A G 7£ nxn , B G ft nxm , C' G ft' xn , 
j4 c € 'R r>cXn ‘, F G 1Z ncXl , K € R. mxn , Ri G 7^ nxn (positive semi-definite), and 
i ?2 G 'R rnXm (positive definite). 

The first order necessary conditions for the quadratically optimal, steady state, 
robust, reduced-order dynamic compensation are the existence of non-negative def¬ 
inite matrices P € R nxn ,Q G 7l nXn , P G ft nXn ,Q G TT*™ satisfying the following 
coupled Lyapunov and Riccati equations: 


0 


0 


0 

0 


PA, + A T t P + £ Af PA, - PjR~ 2 }P, + R x 

1=1 

+ £(A, - Q.V^CifPiA, - Q,V~'C 1 ) + tZpBR; 1 B t Pt ± 

t = l 


A.Q + QA* + Y, A,QAj - Q.V 2 ~ X Q T S + V, 

1=1 

+ j^(A, - BiR^P.MAi - B,R^P.) t + t±QC t V 2 ~ 1 CQtI 

» = 1 

PAqs + A t qs P + PjR^P. - 

ApsQ + QA T PS + Q.VfSQj - t l QC t V 2 ~'CQtI , 


where the projection operator r G 72” xn is given by 


(3.4) 


(3.5) 

(3.6) 

(3.7) 




t= e row 

k=l k 


and Hk[QP] represents the k th eigenprojection of QP, with t ± = I n — r. In (3.4)- 
(3.7), Ai G 7Z nxn , Bi G 7Z nxm , G 7Z lxn are the uncertainty matrices of the plant 


| 
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[2], from which 

A, = A + i ££., A] , B, = B + 1 -4.B, 


C. = C + 1 £?., C, A, 


R„ = R2 + Z: st Bf{P + P)B, 


V 2 , = v, + T.U CAQ + <?)c, T 




P5 = Bj’P + n=, 5f(P + PM, 
Aqs = A, - Q,V{, l C, , 


Q, = QCJ + e:=, MQ + Q)CJ 
Aps = A, - SjiZ^P, , 


where n is the number of sets of uncorrelated uncertainties. 

In terms of the solution to the above MEOP equations, the compensator dynamics 


are specified according to 


A c = r(A, - Q,V 2 - l C, - B,R^P,)G t , 

(3.8) 

F = TQ.V-' , 

(3-9) 

K = R^P,G t , 

(3.10) 

where the operators G € 7Z." cXn ,r E TZ ncXn must satisfy VG T = 

I nc , and G T T = r. 


The complexity of the above expressions is apparent; indeed, this is a trade-off one 
experiences in designing with this method. 


3.2 SCOLE Design 

In this design and simulation study, we are interested in controlling the reflector 
panel, where the control objective was to damp vibrations. A full order system 
dynamics with three inputs and five outputs model was used for the designs to 
follow. The three inputs are the reaction wheel actuators located at the hub while 
the five outputs include three rotational rate sensors located at the hub and two 
accelerometers located at the center of the reflector. 

The performance index to be considered for the SCOLE configuration is the infinite 
time state regulator problem given by: 

J= [ [x T (t)Rix(t) + u T (t)R.2u(t)]dt 

Jo 

where x 6 iZ 20 ; u E R 3 \ R\ € P 20 * 20 , non-negative definite; and JZ 2 € R 3x3 , positive 
definite. The choice of state weighting matrix was made first by transforming the 
original system into its open loop internal balanced realization form. Then the state 
weighting matrix of the balanced system was taken simply equal to the diagonalized 
controllability or observability gramian, and the state weighting matrix for the 
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original system, i2 ls was obtained by transforming the state weighting matrix from 
the balanced system. The input control weighting matrix was chosen to be diagonal 
and with equal weights. 

A time invariant state estimator is considered here. The variance matrices at the 
input and output are assumed to be diagonal and stationary, and are chosen such 
that the poles of the estimator are faster than those in the controller. 

The motivation of MEOP design is to improve the robustness of the OP design. We 
are interested in robustifying the modal frequencies in the model. The uncertainty 
matrix for the modal dynamics was taken as: 


A\ — a 


0 0 
-n 2 d 


9 


and the maximum entropy design parameter, a, was chosen to be 10%. The uncer¬ 
tainty matrices for both the input and the output were set to zero. 


3.3 Results and Simulation 

Both the MEOP and OP design methodologies were applied to the SCOLE configu¬ 
ration. Table 3.1 and Table 3.2 show the robustness and performance characteristics 
of the OP and MEOP designs. The second columns in each table, for the parameter 
|e|, use an improved robustness measure due to Yedavalli [4] for structured pertur¬ 
bations. The parameter, |e|, measured the guaranteed stable absolute upper bound 
for each component which is subjected to perturbation in the plant model dynamics. 
This robustness measure was applied to the closed loop dynamics, but only the sub¬ 
matrix, — ft 2 , of the original, open loop system was assumed to be perturbed due to 
modal frequency variations. The corresponding perturbation on the submatrix, Z), 
was assumed to be negligible due to the natural low damping on the structure. For 
instance, |e| is equal to 0.0154 for the order MEOP design, implying that for each 
diagonal element in ft 2 the maximum variation is bounded by an amount 0.0154 to 
ensure system stability. This corresponds to a —17 to +15% change in frequency 
for the first mode. The Yedavalli robustness measure is a fairly conservative one 
since it gives only a single absolute stability bound for all the components subjected 
to perturbation. However, this is used here as an indicator that the MEOP designs 
have improved the robustness of the OP designs. . 

The third columns show the stable range of the relative percentage change in all 
modal frequency related components in the open loop dynamics matrix from the 
closed loop realization. The significance of this information is in reflecting the actual 
robustness with respect to modal frequency variations. The final columns indicate 
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Table 3.1: Optimal Projection Design 


Order 

|e| 

r ¥(%) 

Cost 

20 

0.0138 

-20 to +4 

0.229 

12 

0.0141 

-30 to +20 

0.231 

10 

0.0153 

—45 to +30 

0.231 

8 

0.0140 

-9 to +30 

0.235 


Table 3.2: MEOP Designs 


Order 

M 

*(%> 

Cost 

20 

0.0148 

-25 to +40 

0.407 

12 

0.0156 

—50 to +50 

0.311 

10 

0.0154 

— 50 to +50 

0.319 

8 

0.0154 

—40 to +40 

0.322 


the steady state cost performance of each design. Clearly there is a trade off between 
the steady state performance and the order of the compensators. The inconsistent 
costs in Table 3.1 between the full order design and the reduced order designs may 
be accounted for by the possibility that the solutions have not yet converged in the 
design. 

An initial condition was imposed on the SCOLE configuration for simulation pur¬ 
poses. Modes 1, 2, 6, 7 and 8 were initially displaced by 0.1 inch. Figure 3.1 shows 
the outputs for the open loop system, where y\ and t/ 2 are the rate measurements 
from the center of the reflector while y 3 and y 4 rate measurements from the hub. 
The closed loop version is depicted in Figure 3.2, in which an 8 th order MEOP 
compensator design is simulated. In Figures 3.3, 3.4, and 3.5, ui, u 2 and u 3 are the 
torque inputs at the hub on x, y, z axes, respectively. With the state space repre¬ 
sentation of the structure in modal coordinates, the vibrational energy of mode i is 
given by 


Ei(t) = + u)*(t) , (3.11) 

where is the modal frequency of mode i and Xi is the general modal displacement 
associated with a given mode in a given direction. 

Figure 3.6 shows the time profile of the vibration energy of the structure for LQG 
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and robustified LQG designs as compared with the open loop case. In Figure 3.7, 
the robustified reduced-order designs (10 t,, and 8 th order) are compared with the 
open loop case. For the designs indicated, all possess excellent vibration damping, 
since 90% of the initial vibrational energy was dissipated in the first second. 


3.4 Discussion 

The MEOP design approach is a very useful tool for control design, particularly 
for large space structures. The order of the compensator can be reduced signifi¬ 
cantly, as compared to the full order LQG design, to one iinplementable in typical 
hardware configurations without much loss in performance. The inherent robust 
characteristics in the compensator design are vital in applications where system pa¬ 
rameters may vary when the structure is, for example, deployed in orbit. The issue 
of spillover or the robustness with respect to unmodeled dynamics is another impor¬ 
tant. subject in large space structures control, and should therefore be investigated 
in this light. 

In the aforementioned designs and simulation, all the reaction wheel actuators are 
assumed to have no physical constraints. From the simulations, the required control 
effort in the first few seconds is very near the limit of the actual hardware. Further 
simulations have shown, however, that when the actuator limit is observed and 
the control gains are constrained, vibration damping is achieved with only slight 
degradation in performance. 
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Figure 3.1: Open loop outputs 
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Figure 3.3: Torque input on r-axis at the hub 
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Figure 3.4: Torque input on j/-axis at the hub 
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Figure 3.5: Torque input on c-axis at the hub 
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Figure 3.6: Vibrational energy profile (full-order designs) 
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Figure 3.7: Vibrational energy profile (reduced-order designs) 
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4. DIRECT ADAPTIVE APPROACHES 


Simply put, the model reference adaptive control approach considers a problem 
where specifications are given in terms of a reference model which describes how 
the process output should respond to a given command. Parameters of the reg¬ 
ulator are adjusted in such a way that the error between the model output and 
the process output, say y m - y, is made small; the difficulty arises in determining 
the adaptation rule. One approach is to use sensitivity-related feedback gains as 
a function of this error. For example, if K(t) represents the vector of adjustable 
controller parameters, then one method for determining the adaptation is to choose 
the controller parameters according to the rule 

-jf = -Hym - y)V K {(y m - y)} , (4.1) 

where k determines the rate of adaptation, and the elements of V/c{ • } are the 
sensitivity derivatives of the error with respect to A'(<). The justification for such a 
rule lies in the assumption that the adaptation parameters are slowly varying rela¬ 
tive to other system quantities, so that the parameters are changed in the direction 
of the negative gradient of the error. Since the controller parameters are updated 
directly, such a scheme has come to be known as direct adaptive control . 

Important issues involved in these investigations are: 

1. The choice of k is not a trivial problem. If chosen small, such an adaptation 
scheme usually performs well, although instabilities may result with a naive 
selection for this adaptation parameter. 

2. For complicated systems, the elements of V# may not be determined exactly, 
thus requiring a certain degree of approximation. 

3. Results are available for extension of the standard theory to nonlinear systems, 
multivariable systems, and nonminimum phase systems. 

The presence of unmodeled dynamics becomes an important consideration in adap¬ 
tive control schemes for flexible structure control. In simulation studies, therefore, 
we typically investigate spillover effects by designing control for reduced-order mod¬ 
els and testing the design on full-order models with simulated disturbances and 
varying degrees of higher-order dynamical effects. 
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4.1 MRAC Schemes 


To introduce the particular type of model reference adaptive schemes being pursued 
for this portion of the work, consider again the general model (2.1)—(2.2), rewritten 
in slightly different notation as 


A p x p + Bu 


(4.2) 


where 



so 


that 


y 
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P 

. *r . 

+ 

. B * . 


0 

fnin 

-n 2 

D 

[o c 2 

] 


(4.3) 


(4.4) 


(4.5) 


where x p € 9R 2n , u £ 5? 3 , y € 3? 3 , and for our purposes here n = 10. 

The basis of z is chosen such that a decoupled set of second order systems or modes 
is formed, fl 2 and D are diagonal submatrices containing the modal frequencies 
and damping respectively of each mode. The modal damping is chosen to be 1% 
for all modes. 


The elements of the input-output vectors are as follows: 


u = 


y axis reaction wheel 
x axis reaction wheel 
z axis reaction wheel 


y = 


y axis rate gyro 
x axis rate gyro 
z axis rate gyro 


The input - output devices are physically colocated on the structure itself, resulting 
in 



(4.6) 
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Table 4.1: Modal Frequencies Used for Control Law Design 


Frequency 

Mode Type 

Axis 

1.51 

Torsional 

z axis 

2.94 

Bending 


4.38 

Bending 



Remarks 

The control laws developed in the following subsections utilize only 3 modes of the 
full order model chosen via a balanced realization analysis, with one input per mode 
required. The modes upon which the designs are formulated appear in Table 4.1 The 
reduced order model is identical in structure to (4.4)-(4.5), except that notationally 
we shall replace x p and x p by x and i, respectively, where now x E W 2n , u E 9? n , 
y E S? n , and n = 3, and where now Q 2 and D are 3x3 matrices. Further, an 
observer is required to generate the state estimates for the state based MRAC law. 
The implementation details will be oinited here For brevity; the equations developed 
in the following subsections assume that the states are directly measured. 


4.1.1 Design 

The MRAC law utilized in the designs requires a model for generation of a state 
error vector e. The model is as follows: 

x m — -f" Bn , (4. {) 


with 

- 

X m — • 

% m 

so that 


im 


0 

TiXn 


2*m 


0 

& m 


ft 2 

D m 


im 

+ 

b 2 


(4.8) 


(4.9) 


where x m E J? 2n , u 6 S 3 , and n = 3. The modal frequencies are those chosen for 
the reduced order plant given in table 4.1, and the damping for all modes is chosen 
to be 5%. 
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The basic assumptions for the adaptive control law implement ation are: 


• The influence matrix B 2 of the plant is not time varying. This allows the 
construction of an adjustable system or combination controller/plant where 
the parameters can be varied. 

• Only vibration suppression control is considered. Therefore, the reference 
input u r is always zero. 

The error between the model states and the plant states is 

(4.10) 

(4.11) 

The generic control law from Hyperstability theory for this formulation is given 
according to 

u — J <f>\idt + <f> 2 X • (4-12) 

The control elements 4> can be chosen from the large class of hyperstable systems. 
For our implementation, 4> is chosen for the the class of Proportional-Integral (PI) 
control laws and Relay—Integral (RI) control laws. Thus, for these classes <j> is 
defined as 


<t>\ = F a v[G a i\ T 

(4.13) 

4>2 = F' a v[G' a x} T , 

(4.14) 

where, for v G 9R 6 , 


v = We 

(4.15) 


and where F a ,F^G a , and G' a are constant positive definite matrices. 



Decoupling PI Implementation 


For the decoupling PI design F a and F a are chosen as 

F = [F. £] 

= [ Onxn 1 ] , 


(4.16) 

(4.17) 
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and G a , and G a are chosen such that 

G a = OtI 2nX 2n ? (4-18) 

G t 1 = a7 2 n x2n > (4.19) 

where a and a are the integral and proportional constants, respectively. Finally, 
the error transformation matrix W is chosen as W = / 2 nx 2 n« 

Given the above definitions, the resulting plant-controller system is decoupled , so 
that in fact 

Xi = 2 Wi(di)ii + w]x x + a [ e t \\x\\ 2 dt + a^-pll 2 , (4.20) 

Jo 

where and e t are individual elements from the vectors x and e, and w t and rf, 
are the natural frequency and damping, respectively, for each mode taken from fi 2 
and D . We see that each mode of the system is dependent on the corresponding 
individual error vector and the Euclidian norm of the state vector. 

Since we have achieved a decoupled modal control, it would be advantageous to 
apply greater control energy to the slow modes or those modes which are “most 
observable” in the output space. This can be accomplished by altering the choice, 
of F a and F a to 

F = [ 0„ xn B^T ] (4.21) 

where T is a diagonal matrix of modal weighting elements. The resulting decoupled 
system with modal weighting is, therefore, 

Xi = 2 Wi(d % )x, + w*Xi + 7 ,a [ e,||i|| 2 d< + 7 i a , e i ||i || 2 , (4.22) 

where 7 i represents the approiate (diagonal) element from the weighting matrix T. 

Simulation studies for this controller design consist of time plots for the input and 
output vectors. The simulations were performed by initializing all the modes with 
a value representing typical output vibrations. For purposes of comparison with 
the simulation results for controller designs in this section and those to follow (all 
MR AC approaches), Figure 4.1 gives the free response (no control) as simulated for 
the three rate gyro outputs. Figure 4.2 depicts the simulated rate gyro responses 
for the three axes with the MRAC PI decoupled algorithm, while Figure 4.3 depicts 
the corresponding input torques of the reaction wheel actuators. We note that the 
control was activated at time < = 1.5 seconds. 
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Figure 4.2: Rate gyro outputs: MRAC PI decoupled 
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Figure 4.3: Reaction wheel inputs: MRAC PI decoupled 
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Coupled PI Implementation 

The second variation on the MR AC design approach we have taken involves a 
proportional-integral type control law without the decoupling structure of the pre¬ 
ceding case. In this case one choice of F a and F a for an input-output coupled control 
is 

F = [ 0 nxn B\ ] (4.23) 

with G a , G a and W chosen as in the previous design. The resulting control law is 
given by 

« =Fa fg e\\x\\ 2 dt + Fa’ellxll 2 (4.24) 

Thus, because of the structure of F the controller-plant equation can be written as 

x = W diag x + Dp diag i + B 2 Bja [ e||i|| 2 dt -f B 2 Bl<x c\\x || z . (4.25) 

Jo 

We have found that this non-decoupled implementation provides superior perfor¬ 
mance over the decoupled implementation for situations in which sensors/actuator 
pairs are not physically colocated. Figures 4.4 and 4.5 represent the output and 
input responses, respectively, for the MRAC PI non-decoupled algorithm. 
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Figure 4.4: Rate gyro outputs: MRAC PI non-decoupled 
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RI implementation 

Development for the Relay—Integral control is essentially the same as described 
for the PI control scheme, with the difference being that the sign of the error 
vector is utilized rather than the error vector itself in the proportional part of 
the control law. Since the decoupled control is useful in the SCOLE context, we 
refer to the development of equations specified for the decoupled PI control for our 
implementation. The final decoupled RI control law, then, is 

Xi= 2w,(di) X, +w?x, + a f e,||x || 2 dt + a's0rj[e,]||x|| 2 . (4.26) 

JO 

Simulations were performed by initializing all the modes with a value representing 
typical output vibrations. Figures 4.6 and 4.7 are the input and output responses, 
respectively, with the MRAC RI algorithm. 
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Figure 4.6: Rate gyro outputs: MRAC RI 
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4.2 Adaptive Model Following VS Control 


The next adaptive technique considered is adaptive model following in which the 
position and velocity outputs are forced via variable structure control to track given 
reference position and velocity paths. This technique for adaptive model following 
was first introduced by Young [5], and employed for flexible structure control by 
Ozguner and Yurkovich [6]. 

A variable structure controller drives the state trajectory to a chosen surface in 
the state space then forces the trajectory to slide along this surface to the origin. 
This is accomplished by switching the feedback gains between sets of feasible values 
whenever the state trajectory crosses the hypersurface in the state space. Variable 
structure control is robust to variations in parameters when in the sliding mode; 
however, the trajectories may “chatter” while sliding to the origin along the surface. 

As mentioned above, in model reference control the objective is to force the plant 
output to track the output of some pre-specified reference model. One way to 
accomplish this is to switch the feedback gains to force the tracking error to converge 
to zero along a desired sliding surface. A block diagram of the model following 
variable structure controller is shown in Figure 4.2. 


4.2.1 Design 

Consider again the model for the SCOLE apparatus given in (2.1)-(2.2). For pur¬ 
poses in this design we assume velocity measurements (colocation) given by 

V — b T q ■ (4.27) 


Note that y and y can be expressed as 
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(4.28) 


The desired position trajectory 9 and velocity trajectory 6 are generated from the 
reference model 


' 9 ' 


' 0 
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' 9 ' 
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a i 

«2 


9 


(4.29) 


where a x and a 2 are design parameters. The tracking error vector is defined as 
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Figure 4.8: Model Following Variable Structure Control 



The feedback control law is given by 


(4.30) 


u = K t e + K l fydt + K 2 y + 6 , ( 4 - 31 ) 

where A%, A”x, K 2 , and 6 are feedback switching gains chosen to force the position 
and velocity error to converge to zero along the switching line 

a = c T e = Cxex + c 2 e 2 = 0 . (4.32) 


It was shown in [6] that the tracking error would converge to zero along the given 
switching line if the feedback gains are chosen according to the following switching 
rules (assuming c 2 > 0): 

€\<r > 0 => K] > ax /{b T b) ; €\<j < 0 ==' hi < Q-\j(b b) ; 


4 








e 2 o- > 0 => K 2 > (cj + c 2 a 2 )/c 2 (b T b) ; e 2 <7 < 0 =t- K 2 < (c, -I- r 2 n 2 )/c 2 (b T h) 

ycr < 0 ==> A'i < ai/(b T b) ; 
ycr <0 =t- A' 2 < a 2 /(b T b) ; 

<r < 0 => 6 < {b T Q 2 q)/{b T b) . 

4.2.2 Simulation Results 

The adaptive model following variable structure controller was applied to the re¬ 
flector portion of the SCOLE configuration in simulation studies, with the control 
objective to damp vibrations. Using balanced realization model reduction tech¬ 
niques, a two mode (Four state) model was developed for each sensor/actuator pair. 
There are five such pairs for the reflector: the accelerometers and cold gas jets (x,y 
directions) mounted in the center of the reflector panel and the rate sensor and 
reaction wheels mounted at the hub (x,y,z directions). The reference model was 
chosen to be second order with a natural frequency u) n = 1 rad/sec and a damping 
ratio ( = 0.1. 

Figure 4.9 shows the control input (unconstrained) for the cold gas jet in the in¬ 
direction, the desired position and velocity trajectories generated from the reference 
model, and the true position and velocity trajectories. Figure 4.10 shows the control 
input for the cold gas jet and the true position and velocity trajectories when the 
physical constraints on the jet are taken into account. In both cases, “chattering” is 
apparent in the true velocity trajectory which is characteristic of variable structure 
control. As expected, the response is degraded somewhat when physical const raints 
are introduced. 

Figure 4.11 shows the control input (unconstrained) for the reaction wheel in the y- 
direction, the reference velocity and position trajectories, and the measured position 
and velocity. Figur 4.12 shows the control input and the true position and velocity 
when the physical constraints on the reaction wheel are considered. 

4.2.3 Discussion 

These simulation studies have indicated several avenues to pursue in future studies 
for the application of switching-type control for flexible structures. First, we have 
seen here and in other studies ([6]) that this technique supplies the desired amount 


ycr > 0 => A'i > a x /(b T b ) ; 

ycr > 0 => K 2 > a 2 /(b T b) ; 
(7 > o =4> 6 > ( b T n 2 q)/(b T b ) 
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of damping as specified by the reference model. The issue of control spillover re¬ 
mains an important question, however, since the chattering effect may introduce 
undesirable effects, depending on the type of actuation hardware used. To remedy 
this situation where needed, time varying sliding surfaces [7] or sliding regions may 
be utilized. Secondly, the actual equipment limitations, taken into account in the 
simulation studies here in terms of amplitude limits, will effect performance. Allow¬ 
able switching speeds must also be considered; indeed, the reflector jets on SCOLE 
respond at frequencies up to about 40 or 50 Hz. 

The simulation example of this report has indicated that the model following scheme 
for flexible structure vibration control shows promise and deserves further investi¬ 
gation. We believe, however, that the true power of the variable structure control 
approach lies in its application to nonlinear plants; many works to date have illus¬ 
trated this point. For this reason, then, investigations focusing on the application 
to nonlinear slew maneuvers in, for example, configurations such as SCOLE are 
warranted. 
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Figure 4.9: Jet input (unconstrained), reference and measured outputs 
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Figure 4.11: RW input (unconstrained), reference and measured outputs 
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Figure 4.12: RW input (constrained), reference and measured outputs 
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5. INDIRECT ADAPTIVE APPROACHES 


5.1 Design 


In indirect adaptive control, the system parameters are estimated using some on¬ 
line identification technique such as recursive least-squares, then the control law is 
computed based on the most recent parameter estimates. For the SCOLE configu¬ 
ration, an indirect adaptive LQ controller was designed to quickly damp vibrations 
in the mast caused by movement of the shuttle body or reflector. This indirect LQ 
controller is discussed in [8] and will only be briefly described here. 

There are four actuators (reaction wheels) and four sensors (accelerometers) mounted 
on the mast. The controller design was based on a reduced order version of the 
model (2.1)—(2.2) which included only the four bending modes. This model was 
discretized to give the state-variable formulation 

x(k + 1) = F d x(k) + G d u(k) (5.1) 

Vacc{k) = H d x(k) + Ju(k ) (5.2) 

Upot{k) — Cx(k) , (5.3) 

for appropriate dimension matrices F d , G d , H d , J. It is assumed that the feedfor¬ 
ward matrix J resulting from the acceleration measurements is known but F d , G d , 
and H d are unknown. It is most likely that for actual testing on the SCOLE, the 
accelerometer data will be integrated for use in the identification process thus elim¬ 
inating the need to know the feedforward matrix. For the purposes of parameter 
identification, the system was put in the input/output form 


Vacc{k) = S T 4>(k - 1) -(- Ju(k) , (5.4) 

where 

9 T = l~Ai — A 2 (Bi + A\J) (I?2 + A^J)] (5-5) 

d> T (k-l) = [yj cc (fc -1) yL(fc-l) u T (*-l) u T (k-2)} . (5.6) 

The matrix 9 contains the unknown system parameters and <f>(k — 1) is a regression 
vector of past input and output values available in real time. Denoting 6(k) as the 
estimate of 9 at time k, the identification algorithm is given by 

9{k) = 9{k — 1) — aP{k — l)f(9(k — 1)) 

P( k- l )<f>(k - l)[y a cc(*) - 9 T (k - l)<f>(k - 1)] 

' 1 +4> T (k -1 )p(k- mV- 1 ) 


(5.7) 

( 5 . 8 ) 
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(5.9) 


P(k) 


p( , . P(fc-l)^(fc-1)(A r (fe-1)P(A:-1) 

1 ’ vl-x + 4?(k - l)P(k - 1)^-1) 


The parameter estimator is a recursive least-squares identifier with data normal¬ 
ization and an additional correction term which allows prior parameter information 
in the form of upper and lower bounds to be included. 

The control law which is based on the most recent estimate of 0 is a state feedback 
control law of the form 


u{k) = -L(k)x{k). (5-10) 

The state estimate is generated from an adaptive observer described by 

x(k + 1) = F(k)x{k) + G{k)u{k) + M{k){y acc {k) - Hi{k) - Ju(k)) , (5.11) 


where M T (k) = [—.di(fc) — ./^(fe)] and 


F(k) = 


- A x (k ) /' 
-A 2 {k) 0 


G(k) = 


Bi(k) 

B 2 (k) 


H = [I 0...0]. 


The feedback gain matrix L(k) is computed from one alteration of a Riccati differ¬ 
ence equation in the manner 

L(k) = [G T (k)R k G(k) + I}- 1 G T (k)R k F(k) (5.12) 

R k+ 1 = Q+ L T (k)L(k) + (F(k)-G(k)L(k)) T R k (F(k)-G(k)L(k)) (5.13) 


5.2 Simulation Results 

With one exception to be discussed later, the output (position) was computed using 
the ten mode model in order to include the effects of the six unmodeled modes. 
Only one set of input/output pairs, that of the reaction wheel and accelerometer 
(x-direction) mounted toward the top of the mast, will be included in this report. 
For all cases, the initial state was chosen to be x T (0) = [0 0 0 0 0 .1 .1 0 .1 .1] (in) 
so that only the four modeled bending modes were initially excited. 

For comparison purposes, a set of baseline plots were first generated. In Figure 5.1, 
the open-loop output position response is shown. In Figure 5.2, both the LQ control 
input without identification and the corresponding output position are shown. In 
this case the identification loop was bypassed and 8(k ) was set equal to the actual 
9 for computing the control law. The effect of one of the unniodeled pendulum 
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modes on the output position is very apparent. Also, (he control inpul greatly 
exceeds the physical constraints on the reaction wheel (1.25 in-lb.) for the SCOLE. 
In Figure 5.3, the LQ control input and the output position are shown when these 
physical constraints are taken into account. 

Next, identification is introduced into the loop. All entries in 0(0), the initial guess 
for 0, are in error by varying values between zero and twenty percent from the true 
values. In the first run, unmodeled dynamics were not included, that is, the output 
position was computed using the four mode “control” model. The unconstrained 
and constrained control input along with the corresponding output position response 
are shown in Figures 5.4 and 5.5, respectively. In the absence of unmodeled modes, 
the adaptive controller performs well in vibration control of the mast. 

Finally, in order to assess the efFects of unmodeled dynamics, the adaptive controller 
was applied to the ten mode model of SCOLE. In Figure 5.6 the unconstrained 
control input and the output position response are shown. The unmodeled dynamics 
cause the closed-loop system to go unstable and the input and output grow without 
bound. In Figure 5.7 the constrained control input and the corresponding output 
position response are shown. The forced bounds on the control input lead to much 
better performance of the adaptive closed-loop system. 

The results, thus far, indicate that a reasonable control strategy for vibration control 
in flexible structures would be to shut off the estimation process after the parameters 
are suitably identified and then use a fixed LQ controller. The estimation error 
would be continually monitored and identification would resume if this error exceeds 
some pre-specified bound indicating a change in system parameters has occured. 


5*04 












506 








Figure 5.3: LQ Control with constraints 
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Figure 5.5: Adaptive LQ with constraints (no unmodeled dynamics) 
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Figure 5.7: Adaptive LQ with constraints (full model) 
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6 3 Robust Indirect Adaptive Control 


The robustness of the adaptive controller with respect to unmodeled dynamics can 
certainly be improved. One option currently being investigated is the incorporation 
of a variable dead zone into the identification scheme for the system parameters 
[9]. The motivation for adding a dead zone into the identification scheme is clear 
if we express the tracking error as e(k) = e,(A;) + e m (k) where c,(k) represents an 
error in identification of the system parameters and e m (k) denotes the error due to 
disturbances or unmodeled dynamics. If the tracking error is due mainly to identi¬ 
fication error, then certainly parameter estimation should continue. However, if the 
tracking error arises mainly from unmodeled dynamics or disturbances, continuing 
with the identification could move the parameter estimates further from the true 
values and thus degrade system performance. The variable dead zone is just an 
attempt to extract a reasonable approximation of the identification error from the 
measurable tracking error using both on-line information and prior system infor¬ 
mation. Thus, whenever the tracking error is due mainly to unmodeled dynamics, 
the identification scheme is “shut off.” A variable dead zone has the advantage of 
simplicity; that is, it can be easily incorporated into most existing adaptive con¬ 
trollers. However, designing the dead zone is not a trivial problem and this choice 
will affect the performance of the system. 


The design of the dead zone has thus far been investigated for some simple exam¬ 
ples which have characteristics of flexible structures control problems. By way of 
example, then, consider a two mode flexible structure described by 


d_ 

dt 
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y{t) = [0 0 0.5 1] 
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(5.14) 


(5.15) 


The reduced order model used for designing the adaptive controller is given by 

— f 91 1 = 0 1 1 f 9i 1 . [ 0 

dt Qi -4 -.04 <7i + 0.5 


y{t) = [0 0.5] 


9i 

4\ 


+ 

J 

+ «m(t) 


«(0 


(5.16) 

(5.17) 


The estimate for the identification error e,( k ) used in the least-squares algorithm 
was: 






c,(k) = 


1 


(5.18) 


f 0 if |c(*)|<rf(fc) 

i e(fc) — d(k) if e{k)>d(k) 

[ e(k) + d.(k) if e(fc) < -d(k) 

where d(k) represents an upper bound on the magnitude of the modeling error e m (A*) 
and is given by 

m(k) = (0.951)m(fc — 1) + u(k) (5.19) 

d(k) = (1.026) \m(k) | . (5.20) 

The control law chosen was the adaptive LQ control law previously discussed for 

SCOLE. 

Figure 5.8 shows the control input and the resulting output when the parameters 
of the reduced order model are assumed to be known exactly. The output was 
computed from the full order model. 

Figure 5.9 shows the control input and resulting output when estimation is included 
but a dead zone is not incorporated into the least-squares algorithm. The initial 
parameter estimation errors ranged from 10—20%. As in the case of SCOLE, the 
unmodeled dynamics cause the adaptive control loop to go unstable resulting in an 
input and output which grow without bound. 

Finally, the dead zone described above was incorporated into the estimation scheme. 
As shown in Figure 5.10, both the control input and the resulting output converge 
to zero. 
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Figure 5.8: LQ control (unmodeled dynamics present) 
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Figure 5.9: Adaptive LQ control without the dead zone 
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Figure 5.10: Adaptive LQ control with the dead zone 
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5.4 Discussion 


Initial results on the incorporation of a variable dead zone into the least-squares 
estimation scheme look promising. Current work centers on the design of the dead 
zone for SCOLE. Furthermore, a better choice of modes will be made in the re¬ 
duced order model for the SCOLE to reduce control and observation spillover. As 
mentioned, one of the unmodeled low order pendulum modes (0.08 Hz) appears in 
the output position response even in the non-adaptive case. 

Another approach for improving robustness is also being investigated. The empiri¬ 
cal transfer function estimation (ETFE) method will be used to identify the system 
parameters [10]. Using finite length discrete Fourier transforms (DFTs) of the in¬ 
put/output data, the ETFE method provides both an estimate of the plant transfer 
function and a frequency-dependant upper bound on the modeling error between 
the estimated transfer function and the actual plant transfer function. This identifi¬ 
cation scheme has been shown to possess certain desirable robustness characteristics 
for use in adaptive control [11]. Because it requires DFTs of the input/output data, 
the ETFE method will require more computational time than the variable dead 
zone identification scheme. However, Tzes and Yurkovich have shown that with 
modest computational load (less than 10 ms) the ETFE scheme may be carried out 
on-line, with a recursive least squares estimator, to accurately identify frequency 
response characteristics for a slewing flexible beam [12]. Moreover, it is shown in 
[12] that the frequency response information thus obtained is useful in designing an 
LQ-adaptive controller with frequency weighting; details of this work will appear 
in a future report. 
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0. LABORATORY EXPERIMENTATION 


0.1 Visits to Langley 

May 1987 

In May of 1987 we visited the LaRC facility and worked with Mr. Jeff Williams for 
two days in the laboratory preparing for control implementation. Our progress was 
limited due to unforeseen problems with the computer and compiler systems. We 
were, however, successful in familiarization with the current set-up, and have since 
completed the designs described herein, and have completed coding for use on the 
structure. 

December 1987 

The next visit to Langley was in December of 1987 for five full days of imple¬ 
mentation and experimentation. All control designs and necessary software were 
transported via floppy disks, so that the on-site work began by downloading the 
various control law programs and their data files for compilation on the the SCOLE 
computer, a Charles River 68000 based system. A variety of differences between the 
Vax 11/785 FORTRAN code and the Unix/Green Hill FORTRAN compiler were 
identified so that first two days were consumed in making the code changes, com¬ 
piling, and testing the simulator code. During the software preparation, hardware 
problems were also delineated. For example, the reaction wheels at the reflector/- 
base had not been used before, and several wiring checks had to be carried out. 

The third day began with experimentation involving the algorithm with the least 
computational requirements, the observer/pole placement design. The initial plan 
was to run this algorithm at 100 hertz but the computational requirements caused 
the sample time to be exceeded before the control law was complete; the algo¬ 
rithm did, however, operate at 40 hertz without a time exceeded warning. We 
had over estimated the Floating Point Operations per Second (FLOPS) capability 
of the computer, primarily because the previous hardware configuration was such 
that the floating point operations were done by a Sky floating point board. This 
board has capabilities of 3 micro second operations or 330,000 FLOPS. Since then, 
however, the computing environment for the SCOLE set-up has been changed to a 
68020/68881 processor, and the compiler has also been changed. Because of this un¬ 
foreseen problem, some time was spent benchmarking various operations; it appears 
that the best expected performance is now 100,000 FLOPS. Utilizing the Sky board 
again, with the current configuration, did not result in increased performance. In 
the interest of time, we terminated our efforts at increasing FLOP speed at that 
point. 
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Because of the unexpected slower speed capabilit.es, most, of the control laws brought 
were too slow in their present form, since most of the algorithms used continuous 
time observer equations. The algorithms can, however, be made more efficient, 
and current efforts are directed in this vein. The inirect adaptive control designs, 
which are already in discrete time form, may require more FLOP capability than 
the SCOLE computer can provide at. the present time. 

Two algorithms were marginally acceptable in terms of achievable sampling period, 
specifically, the pole placement and MEOP algorithms, although the maximum 
achievable sampling rate for the MEOP design was 10 hertz. Before experimentation 
with these control laws, various tests were run to determine if the amplifers were 
functioning properly and if wire phasing was correct. Due to concerns with the 
amplifiers, the reaction wheel actuators were run at 5 in-lb or 30% of the maximum 
rating. 

Both of the aforementioned algorithms showed some undesirable behaviour. To 
determine some rationale for the poor results, the outputs from the observers were 
recorded, indicating that several of the modes which should be pure sinusoids were 
heavily distorted. This in turn indicated that noise and observation spillover may 
be responsible for poor observer convergence. The SCOLE Sensor conditioners are 
two pole low pass filters with a cutoff frequency of 10 hertz. Our state space model 
of the SCOLE structure upon which the observers were designed includes modal 
frequencies of a maximum 5 hertz. Observation spillover from the unmodeled modes 
in the 5 to 10 hertz range may account for some of the distortion. Channel noise 
seemed to be significant relative to the signal input data. A data sample of channel 
noise was collected and brought back to determine noise power and spectrum. 

On day five, system identification experiments/data collection and a simple control 
law (described below) were priorities. To identify various input output relationships, 
white noise inputs were applied to the various actuators individually. All the outputs 
were collected for each actuator run. A thousand data points were collected at 
sampling rates of 20 or 50 hertz. During this data collection, we noted that the 
low authority reaction wheels mounted near the top of the mast provided little 
excitation of any of the modes. This data is currently being used for investigation 
of various system identification techniques, primarily on-line schemes such as the 
above-mentioned ETFE approach. 


0.2 Software 

The software to implement the control algorithms described in previous sections was 
developed with the following goals. One, the code must be implementable on OSU 
computers (Vax 11/785 and MicroVax II) simulating all experimental hardware 
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dependent subroutine calls (A/D, D/A, timers, and so on) with code to run a state 
space model of SCOLE. The control law code is validated by this process since the 
SCOLE computer debugging facilities are limited. Also, the simulation produces 
data for comparison to actual data output. The same code used for the simulations 
is used to run the experiment by linking in the runtime libraries instead of simulation 
libraries. 

Secondly, the software is designed with standard hooks to the control law subroutine. 
This allows one standard interface to the experimental hardware to be designed, 
implemented, and validated. To implement a control law, only two subroutines 
must be written, interfacing to the main logic through sensor-in actuator-out data 
arrays. All other tasks are performed by the standard interface. These tasks are 

• Set sample rate, number of samples, sensor bias levels, and so on; 

• Check if control law calculation time exceeds sampling interval; 

• Sample analog inputs; 

• Remove input bias, reorder the data array to standard model order, and 
convert to usable values; 

• Integrate all sensor values; 

• Scale and limit output values; 

• Convert to DAC array ordering, output analog outputs; 

• Save all input and output values. 

6.3 Results 

To complete our most recent visit, we chose to implement a simple rate feedback 
law, consisting of a proportional feeback term for each rate sensor and reaction 
wheel actuator in each plane. The feedback gains were chosen by exciting the 
modes in each plane individually and choosing appropriate values according to the 
observed behaviour. The final data collected was done by exciting the system in all 
three planes and running all three actuators simultaneously. Table 6.1 shows the 
actual feedback gains. From the table, we note that the X direction actuators were 
relatively sensitive and could not tolerate large gains. Figures 6.1 and 6.2 show the 
actual sensor outputs and actuator inputs, respectively, for the rate feedback law 
to an initial disturbance. 
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Table 6.1: Feedback Gains fo r A iate Feedback Control Scheme 



X Axis 

-5 

Y Axis 

Z Axis 

-100 

-50 


The results from this experiment imply that the output from the X axis rate feed¬ 
back gyro is always positive, not a sinusoid oscillating about zero. This would 
essentially indicate that the position error is positive and unbounded, suggesting a 
failure either in the rate gyro itself or in the measurement process 1 . This probably 
accounts for some of problems encountered with the more sophisticated algorithms, 
and may account for the low feedback gain required on this axis during the rate 
feedback experiment. The performance in the other axes was satisfactory. 


1 Jeff Williams has been made aware of the problem and is investigating. 


52 ) 



Seconds 


X Axis Rate Gyro 



i. 2. 3. 4. 5. 8. 7. 8. 0. 10. 


Seconds 

Figure 6.1: Actual sensor outputs, Rate Feedback design 
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Figure 6.2: Actual actuator inputs, Rate Feedback design 
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7. COFS CONTROLLER DESIGN 


7.1 Introduction 

In spite of its cancellation, the original Mast Flight System provides (and will proba¬ 
bly continue to provide) an interesting model that can be used to test in simulation 
different control approaches. Finite element models of varying complexity exist, 
experiments have been performed on the individual actuators. 

The basic element in the Mast Flight System is a 60.7 meter long, triangular cross 
section truss structure. The truss has 54 bays with LDCM actuators distributed 
along the length. There are four primary actuators located at the tip of the beam 
so as to produce both x and y axis forces and z axis torques. The actuators at the 
intermediate stations are smaller devices and can be utilized to check decentralized 
control issues 1 . 

Our effort in the Control Research Lab of the Ohio State University was concerned 
with three issues: 

1. Design of an optimal decentralized controller. 

2. Design of a decentralized controller through overlapping decompositions and 
loop-shaping. 

3. Design of low-order controllers through fixed—order, direct optimal projection 
techniques. 

Throughout this section we will be referring to certain data which was supplied to 
us at different times from NASA LaRC. These may not be the same versions that 
were finally accepted and certified or utilized in NASA’s final design effort. Be that 
as it may, the data used here is still very representative of such a flight system. 


7.2 Design of an Optimal Decentralized Controller 

7.2.1 Background Theory 

Consider 

x = Ax + 53 BiUi ;.t(0) = x 0 (7.1a) 

1 = 1 

x In the analysis to follow identical actuators have been assumed. 
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y, = C; i 

with cost 




,00 v v 

J= ( x t Qx + 2 ^ x T N l u l 4 - Y, u\R x u t ) dt 

Jo i =1 t=i 

and the following feedback structure constraint: 
ut = Ki yi ;i = . 

It can be shown [13,14] that the solution of the above optimal control 
reduces to solving the following system of non-linear algebraic equations: 


(7.1b) 


(7.2) 


I 


AjP + PA C + Q = 0 
A C L + LA T C + X 0 = 0 


(7.3) 
problem 

(7.4) 


and 


V*. J = Bf PLCj + RiKiCiLCi + Nf LCf = 0 ; i = 1,..., v 


T 1 atT t r<T 


where 


A c = A+ Y B.KiCi 

i = l 

Q = Q + Y, CjKjR x K t C\ + Y (N,K X C, + CjKjNj) 


i~ 1 

Xq ~ £0 *o 


t=i 


(7.5) 

(7.6a) 

(7.6b) 

(7.6c) 


Many approaches have been introduced in the literature to solve the above system 
of equations. The approach pursued in our software package 2 is a gradient method. 

Extension to the basic problem can be done by considering the following decentral¬ 
ized dynamic stabilizing compensator: 


Z{ — GiVi 


(7.7) 

Ui = "I - N{y t 

; * = 1 , • • -, v . 

(7.8) 


2 The software package DOLORES was developed by Professor U. Ozgtiner, Mr. Khorranii and 
Mr. Tien for another project 
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The problem of dyrv.r output feedback is reduced to the previous problem by 
augmenting the state space of the system and the stabilizing compensators, that is, 


X 




X 


' B, 

0 

z \ 


’ A 

0 


* . . 




z 2 

— 




Z\ 

z 2 

V 

+ L 

0 

0 

; 


o ; 

0 



i=1 

• 

I 

- . 




Zv 


0 

0 


u, 

v, 


(7.9) 


y. 


C, o ... o 
o o Ci o 


*2 


(7.10) 


where Vi = and C x is provided by the user. The cost may be further modified to 
J = \ r (x t Qx + 2J2 + £ njRju, + £ ifi?, 2 i,) dt . (7.11) 

^ Jo <=i i=i »=i 

7.2.2 Design of an Optimal Decentralized Controller 

The modes and mode shapes used in this study are listed below: 


Mode 

Freq. (Hz) 

Bay 28 

Tip 

(1) 1st x-z 

0.187 

.0149 

.0486 

(4) 2nd x-z 

1.36 

.0647 

-.0243 

(7) 3rd x-z 

3.93 

.0111 

.0167 

(10) 4th x-z 

6.84 

-.0567 

-.0107 


Each of the actuators with their associated (inner-loop) compensators (already 
designed) can be described by a fourth order model: 


»rc A c x c ~b B C1 6 C 4" 

/ = C c x c + D Ci 6 c + D CJ p 


(7.12) 

(7.13) 


where, 


I 
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(7.14) 


A c = 


-T 0 0 0 

0 -T (1 - g 2 )r o 

0 0 0 1 

Ka _ Ksi _ 31 Kg 0 


B et = 


(1 - 3i) r 
0 
0 

gi Ko 


Be, = 


0 

-(1 ~ 9i) T 

0 

92 Kq 


C c = | —h 0 ho g 2 h o 0 j 


(7.15) 


(7.16) 


(7.17) 


D CI = -giK 0 (7-18) 

D CJ = -g 2 K 0 . ( 719 ) 


In the above, 8 C is the control input, p is the displacement of the structure (in the x 
direction) at the location of the actuator, and / is the force output of the actuator 
as applied to the structure. The parameters used in our study are, t = 9.9298s , 

gi = 3, g 2 = 6, A' 0 = 83.333 and m = 11.4 kg. 

A two-mode model is used to describe the x-z bending plane dynamics of the 
structure: 


X. = A.x. + B\r + B 2 ,*f 2S 
V' = I>. 

y* = p i = C i ,x, + D'; t f t + D', 2 »f 26 


(7.20) 

(7.21) 

(7.22) 


where i takes on the values w t” and “28” and these denote the quantities at the tip 
and bay 28 locations, respectively. 

The equations given above can now be combined to give a state-space model with 
12 states. In doing this 0.2 and 0.3% damping was assumed for the structure. 
A quadratic regulator was designed with second order dynamics admitted for the 
compensators in the feedback controllers. 
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7.3 Decentralized Design using Overlapping Decompositions and Loop- 
Shaping 

7.3.1 Background Theory 


Consider a system with two decentralized control agents described by 


x = Ax -|- Bi ui + B 2 u 2 , 

(7.23) 

Vi - 

(7.24) 

V 2 = C 2 x 

(7.25) 

where x E R n is the state, U{ E R m ' is the input, and yi E R m ' 
i th control agent (i = 1,2). Let the state x be partitioned as 

is the output of the 


*1 


x = 


*0 


*2 


(7.26) 


where x< £ R n ' (i = 1,0,2). The partitioning is usually done in such a way that r, 
corresponds to the part of the state space which is strongly observable and control¬ 
lable only by the i th control agent (i = 1,2), and x 0 corresponds to the part which 
is strongly observable and controllable by the both agents [15]. 

Consider the transformation: 


x — Tx = 


Xi 

x 0 


x 0 

x 2 


<E R ", 


n = n + n 0 . 


(7.27) 


The expansion of the original system (7.23) with respect to the transformation (7.27) 
is given by (see [15] for developments) 


x = Ax + Bu, i(0) = Tx 0 E R n 


(7.28) 


y = Cx 


(7.29) 


where 


A = 


A i 
A 21 


An 

A2 


= TAT 1 + M a , 


S2Q 


(7.30) 




B± 

*i 

B 2 \ 

B\ 2 

b 2 

= TB + A/b , 

(7.31) 

C = 

' Cl 

. Cai 

Cl 2 
C 2 

= CT 1 + A/c , 

(7.32) 


T 1 is a generalized inverse of T satisfying T J T = /, and A/ 4 , A/b, and Me are com¬ 
plementary matrices satisfying appropriate inclusion conditions. It can be shown 
that the response x(f) of (7.23) is related to the response i(<) of (7.27) by 

x(0=T 7 x(<) V<> 0 . (7.33) 

Furthermore (7.23) and (7.27) have equivalent input/output descriptions 

G{s) = C(sl - A)-'B = G{s) = C(sl - A)~'B . 

Now consider the uncoupled expanded model described by 

x = Ax + Bu 
y = Cx 

where 

A = 0 

0 a 2 

5 _ \B X 0 

[0 B 2 

r _ \ Ci 0 

6 0 c 2 . 

The transfer function matrix (TFM) for the uncoupled expanded model 

G{s) = C{sl - AY'B = block diag (Gi(a),G 2 («)) (7.40) 

where G,(s) is the TFM of the local model 

= A{£i + B,Ui ( 7 - 41 ) 

Vi = Cii> , ( 7 - 42 ) 


(7-34) 

(7.35) 

(7.36) 

(7-37) 

(7.38) 

(7.39) 
is given by 
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for agent i (i = 1,2). Note that this model differs from one that is obtained I>y 
simply ignoring the interactions in the original system as it retains portions of the 
state space which affects each channel. 

Let G(s) and G(s) be related by 


G(s) = G(s)(I + E(s)) 


(7.43) 


where E(s) is the so-called multiplicative error matrix between the true TFM Ct(<s) 
(or equivalently G(s)) and the uncoupled expanded TFM G(s). An upper bound on 
the norm of E(ju>) can be found as 


<r{E{ju>)) < 


- 6W) 



(7.44) 


where <r(-) and <r(*) indicate, respectively, the maximum and the minimum singular 
values of the indicated matrices. 

In practice, model (7.23)-(7.25) may not represent the physical system exactly, that 
is there may be some modeling uncertainties present. In such a case, a total error 
function, e m (u>), must be defined, for example as 


e m (u>) = e(u>)e c (u>) + c 0 (w) 


(7.45) 


to represent the error between the uncoupled expended system model and the true 
system. Here, e(u>) represents an upper bound on 5(E(ju? )), such as the one given in 
(7.44), e c (u>) > 1 accounts for modeling uncertainties in interactions, and e 0 (u>) > 0 
accounts for uncertainties in subsystem models. 

Once the total error function is chosen as in (7.45), our approach proceeds with 
designing decentralized controllers based on local models (7.41)-(7.42). These con¬ 
trollers are designed such that the local stability robustness requirement [15], 

a(T,(ju)) < —| r , Vu; (7.46) 

and acceptable performance requirements are satisfied at each local station. Here 

TM = (I + GMKM)-'GMK,(s) , 
is the closed-loop design TFM for the i th agent. 
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i = 1,2 


(7.47) 
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7,3.2 Design of the Decentralized Controllers 


An initial try at decomposition was made with the dynamics of the individual 
actuator/sensors assigned to separate subsystems and the two-mode model being 
kept in the overlapping portion. Thus the subsystem states are described by: 



(7.48) 


This is a design choice, although there indeed are guidelines in the selection process. 
In fact, with a larger model for the structure, portions of the structure model (for 
example, the “flexible modes”) can also be assigned to separate subsystems as based 
on a measure of effectiveness. In a structure like a beam this would be directly 
related to the mode-shape values at the individual actuators. 

An upper bound on the error introduced due to the interactions between the sub¬ 
systems is determined as 

e(u>) = ---r . (7.49) 

Furthermore, it is assumed that the modeling uncertainties can be represented by: 


e c (u>) = 



and 


eo(ui) 


“((%)* + »)* 

((I0u>) 2 + 1)» 


(7.50) 


(7.51) 


These functions, together with the final e m (u>) given by (7.45), are plotted in Fig¬ 
ure 7.1. 

With the present decomposition, Linear Quadratic Gaussian controllers were de¬ 
signed for each subsystem to obtain good performance while satisfying the stability- 
robustness requirements. Gaussian white noise processes of intensity 1 and 100 
are assumed to be present respectively at the input and the output of each sub¬ 
system. Local quadratic performance indices are chosen with control weightings 
= p 2 = 100, and state weightings Qi — LjL, (i = 1,2), where 


S"3i 



Lx = 


10- 3 C c 0 


0 


ioon 


(7.52) 


and 


L 2 = 


0 10" 3 C c 

ioor 38 o 


(7.53) 


The resulting control and filtering gains are: 


Ki = [.0297 -.0297 

-. 0.888 

.2329 -.2171 .3369 .1076 . 0103 ] 

( 7 . 54 ) 

K 2 = [-.0625 .0736 

-.4039 

-.0384 .0147 -.0147 -.0439 . 0409 ] 

( 7 . 55 ) 


H, = 


' -.0764 ' 


' -.0033 ' 

-.0272 


-.0315 

.0104 


-.1281 

.0917 

, h 2 = 

-.0862 

-.0233 

-.1445 

-.0509 


-.0135 

.0361 


.0101 

.0101 


.01817 


( 7 . 56 ) 


Thus the individual controllers are of the form 


Uj = — A, z, 


(7.57) 


where 


Zi = (A* - B t Ki - H,Ci)zi + Hi(y, - D { Ui) . (7.58) 

Magnitudes of the resulting closed-loop TFM’s are plotted together with the sta¬ 
bility bounds in Figures 7.2 and 7.3. It is observed that the stability robustness 
requirement (7.46) is indeed satisfied. 

As a final check, the controllers were applied to a 4-mode “truth model". The re¬ 
sults are summarized in Table 7.1. It is observed that 13% damping in the first 
and 5% damping in the second mode were achieved. Furthermore this was accom¬ 
plished with relatively small control and filtering gains. The design is also robust 
to modeling errors and plant variations. The closed-loop system is guaranteed to 
remain stable as long as the perturbations are bounded by e m (u>). 
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Table 7.1: Decentralized Control Design Results 


Mode 

Structural Eigenvalues (damping, %) 

Freq. 

(Hz) 

OPEN-LOOP 

CLOSED-LOOP 


With actuators 

Design Model 

“Truth” Model 

0.187 


-0.0038 ± j'1.1556 
(0.33) 

-0.1511 ± jT.1724 
(12.78) 

-0.1514 ± j'l. 1727 
(12.80) 

1.36 


—0.1342 ± . 78 .6144 
(1.56) 

-0.4862 ± j'8.5773 
(5.66) 

-0.4396 ± j'8.6359 
(5.08) 

3.93 

-0.1235 ± j'24.693 
(0.5) 

-0.1248 ±;24.696 
(0.5) 


-0.1236 ±;24.6972 
(0.50) 

6.84 

-0.2149 ± j'42.977 
(0.5) 

-0.2186 ±;42.995 
(0.51) 


-0.1954 ±;43.0042 
(0.45) 


S3G 







7.4 Direct Low—Order Controller Design 

In this section we discuss application of the direct (fixed-order) design approach 
to reduced order controller design for the COFS mast. The details of the optimal 
projection design methodology, with the maximum entropy stochastic modeling of 
the uncertainties (MEOP approach), were discussed in earlier in this report. We 
thus omit further detail here, and refer the reader to Equation (3.1), (3.2) for the 
general system and controller structure, Equation (3.3) for the cost criterion, and 
Equations (3.4)-(3.7) for the coupled Riccati and Lyapunov equations to be solved 
in the MEOP approach. 

In the design results to follow, the model employed is as given above, for a full 
10 mode representation. We note, however, that the designs of this section do not 
include the effects of actuator dynamics. 


7.4.1 LQG Parameter Selection 

For comparison purposes we first present the results of an LQG design (we omit 
details here) involving the full 20 state model. The design objective is to design 
compensation to dampen oscillation of the COFS mast, structure as much as possible 
while avoiding high gains. The four design parameters for this exercise are R i, R 2 
(see Equation (3.3)) in the controller Riccati equation and Vj, V 2 in the filter Riccati 
equation. 

After trying severed different combinations of the design parameters and conducting 
simulations, the following parameters were chosen: 


Rj = 1(T 2 


0 0 
0 R 


R 2 = 10- 3 / 4 x4 


(7.59) 


where 


R = diag { 20, 20, 10, 10, 5, 5, 5, 1, 1, 1 } , 


(7.60) 


and 


Vi = 


0 0 

0 120x20 


, F 2 = 10- 2 / 4 x4 . 


(7.61 


The results of this design (referred to later as the “20th order controller”) are 
presented in comparative plots of the next section. 
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Table 7.2: LQG & OP Controllers 


Order 

j 

20 

16 

14 

12 

10 

8 

6 

Cost j 

2.848 

2.848 

2.849 

2.849 

2.867 

2.895 

2.932 


Table 7.3: LQG & MEOP Controllers 


Order 

20 

16 

14 

12 

10 

8 

6 

Cost 

3.681 

3.561 

3.561 

3.562 | 

3.565 

3.568 

3.804 


7.4.2 Reduced-Order Design Results 

For the optimal projection design, Table 7.2 lists the cost for each design, begin¬ 
ning with the full 20th order controller, including OP designs of order 16, 14, 12, 
10, 8, and 6. Taking into account the vibrational energy of the mast flight system, 
controller designs with orders down to about 12 provide acceptable damping rela¬ 
tive to the full LQG design; the energy profiles of each controller design, and the 
open loop case, are given in Figure 7.4, where the energy expression is given by 
Equation (3.11). For the MEOP designs, where robustification (using the above- 
mentioned noise intensities) was carried out with respect to the modal frequencies 
of the first bending modes in all three planes, the results are given in Table 7.4.2 
and in Figure 7.5. Again, a 12th order compensator is acceptable considering both 
the cost and the ability to dissipate energy. 
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8. FUTURE DIRECTIONS 


We have identified several areas to pursue for continuation in this research, not 
only for the SCOLE configuration but for the MiniMast apparatus as well. These 
directions for future investigation are summarized below. 

Nonadaptive Control Approaches 

The issue of spillover or the robustness with respect to unmodeled dynamics is 
another important subject in large space structures control. Large flexible structures 
are inherently infinite dimensional, but we can only work on a finite dimensional 
model due to the constraints in computer capability. In this regard, the system 
response to unmodeled dynamics should also be investigated in the future. We 
are, therefore, focusing much of our attention in current and future design efforts 
towards robust controller designs. 

To this end, we intend to investigate asymptotic LQG synthesis procedures, such as 
LQG/LTR, as well as the MEOP design approach. We continue to believe that the 
MEOP design approach is a very useful tool for control design in the application area 
of large space structures. This is due to the fact that the compensator order can be 
reduced significantly from the full order LQG design, allowing easy implementation. 
A comparison between the various techniques in terms of stability robustness to 
variations in all system parameters (for SCOLE) should be carried out. 

Adaptive Control Approaches 

Preliminary results presented in this report have shown that the indirect adaptive 
LQ controller is initially very effective in vibration damping but effects of unmod¬ 
eled modes can eventually cause the control input and output to grow without 
bound. The variable dead zone method and the empirical transfer function estima¬ 
tion scheme are being pursued as viable options for improving the robustness of the 
indirect adaptive LQ controller. 

Another method for improving robustness is to introduce an adaptive filter in order 
to reduce the effects of control and observation spillover. Methods for designing 
this filter on-line are currently under investigation. Preliminary results on simple 
models indicate that the addition of filtering improves both the transient response 
and the identification of system parameters. 

S ystem Identification for Contr ol 

The empirical transfer function estimation (ETFE) method is under current, inves¬ 
tigation for real-time identification of general ARMA model formulations. Using fi¬ 
nite length discrete Fourier transforms (DFTs) of the input/output data, the ETFE 
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method provides both an estimate of the plant transfer function and a frequency- 
dependent upper bound on the modeling error between the estimated transfer func¬ 
tion and the actual plant transfer function. We have obtained new results along 
these lines for on-line estimation which enables the design of a frequency weighted 
optimal LQ-adaptive controller. We believe that the ETFE approach is particu¬ 
larly well suited to flexible structure identification and control due to the frequency 
response profiles typical of such systems. 

As for off-line identification techniques, we are investigating experiment design and 
input synthesis algorithms for flexible structure identification in a related project 
within the Control Research Laboratory at Ohio State. We therefore hope to in¬ 
corporate many of those ideas in future work with modeling and identification of 
SCOLE and MiniMast. 

Variable Stru ctur e Control 

Due to the dynamic coupling of rigid body motion and flexible modes during slew 
maneuvers, as well as to the inherent complexities resulting from nonlinear dy¬ 
namics, we believe the application of Variable Structure Control to be extremely 
promising for control of slewing. The SCOLE apparatus and accompanying model 
offers a suitable testbed for such ideas. We have examined the feasibility of VSC 
for vibration suppression in our previous work on the NASA Grid and SCOLE, and 
based on the our experience we propose to investigate such techniques for control of 
slewing on the SCOLE model. Indeed, experiments to date have been performed at 
NASA/LaRC [16] employing “on-off”, bang-bang type controls on the apparatus, 
and several works have recently appeared advocating such ideas for general rapid 
slewing, reorientation maneuvers of large space structures [17]. 

Time Optimal Control 

Time Optimal slewing of the SCOLE structure presents some extremely interesting 
problems. Initial experimental results have been reported by Montgomery et al 
[18] where the attitude control problem for SCOLE was treated as time—optimal 
uncoupled pitch and roll motions. The model that the basic design was based on 
included only the rigid body modes; the effects of the vibrational modes were dealt 
with by an ad-hoc scheme. Recently Barbieri and Ozgtiner [19] have been able 
to derive the time-optimal feedback control for systems with both rigid body and 
purely vibrational modes. The SCOLE configuration will be utilized as an example 
for evaluating these control strategies. 

Experimentation 

From our most recent visit to LaRC and subsequent evaluation, we have determined 
appropriate actions to take for improvement. Since a large portion of the difficulties 
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experienced centered on the computing power available, particularly evident in the 
observer-based designs, we will be investigating the following solutions. After noise 
problem analysis, we propose to install a Butlerworth filter to narrow the effec¬ 
tive sensor bandwidth and eliminate excess sensor noise and unmodeled dynamics 
(modal frequencies). In order to speed up computation time, we plan to add state 
space discretization routines to discretize the continuous time observer equations 
at runtime. This will essentially result in double to triple the current maximum 
sample rates allowing us to try the MRAC and MEOP control laws again with a 
much faster sampling rate. Moreover, we need to delineate the problem with the 
A^-axis rate gyro and take corrective action. 

The experimental software program is installed on the SCOLE computer. With the 
cooperation of a technician at Langley, we are now prepared to run experiments over 
the phoneline. In the coming months we intend to take advantage of this capability, 
provided that help is available at LaRC. With regard to the computer setup in the 
laboratory, we feel that the overall procedure could be improved with the availability 
of some sort of local plotting routines, or graphics capabilities, to evaluate control 
law performance at runtime. Also, some sort of diagnostic programs to evaluate 
system behaviour before any control laws are tried would be a tremendous help. 
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ABSTRACT 

This paper documents the linear least-square 
identification procedure used to obtain an 
empirical model of the vibrational dynamics of 
SCOLE, a laboratory apparatus used to test 
parameter identification techniques and control 
laws for large, flexible space structures. 

Testing is done autonomously by exciting the 
structure from a quiescent state with torque 
wheels and recording the time history data of rate 
gyro sensors. The torque wheels are then shut 
down and free-decay data is recorded. The 
free-decay portion of the data is analyzed using 
the Fast Fourier transform to determine the best 
model order to use in modelling the response. 
Linear least-square analysis is then used to 
select the parameters that best fit the output of 
an Autoregressive CAR) model to the data. The 
control effectiveness of the torque wheels is then 
determined using the excitation portion of the 
test data, again using linear least squares. This 
report describes the system model assumed and the 
experimental apparatus and procedures used. Also, 
typical experimental data are presented that 
reflect the performance of the identification 
algorithms. 


INTRODUCTION -. 

Control of distributed-parameter systems is a 
technology that will be required for future 
spacecraft, but there are both theoretical and 
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t VPI&SU Aerospace Engineering Co-op, Spacecraft 
Control Branch. 

** Computer System Manager, PRC Kentron, Inc. 


practical problems which must be overcome to 
create an effective design capability within the 
control systems community. Recognizing this, NASA 
has built a national research facility, the 
Spacecraft Control Laboratory Experiment (SCOLE) 
(1,2], to support this research effort. This 
facility provides researchers with a highly 
flexible test article, sensors, actuators, and 
digital control processing capability. The test 
article mimics the Space Shuttle with a large, 
flexible, offset-feed antenna attached to the 
payload bay (figures 1 and 2). 

Much interest has been expressed by the 
research community concerning SCOLE. This Is 
reflected in the technical output of two workshops 
held at the NASA Langley Research Center [3,41. 
Using SCOLE, control laws for a multi-input/output 
structural dynamics system can be implemented in 
real time from any remote site that has a computer 
terminal and modem communications capability. In 
order to design effective control systems, the 
technical community needs a good model. 

Previously, work has been reported on the 
empirical determination of the rigid body 
rotational dynamics of SCOLE (reference 5). To 
this time, however, there has been no 
comprehensive testing of its vibration dynamics. 

The technique used to identify the dynamics 
of SCOLE is an extension of the one developed in 
[6). The model used to represent the dynamics is 
an Autoregressive Moving Average (ARMA) model. In 
[6] the test data was assumed to contain only a 
single mode and the test input was modified to 
produce only that mode. That technique is valid 
for structures wherein that can be achieved. 
Unfortunately, for SCOLE, generating a response 
that contains only one mode is very difficult. 

Thus, a new procedure had to be evolved. The 
object of this paper is to develop and tutorially 
describe a method applicable to SCOLE that Is 
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capable of handling multiple inodes in the test 
data. 


The report is organized as follows: first, 
the experimental apparatus and the models used to 
represent the vibration dynamics are described. 
Then, the parameter identification procedure and 
computational algorithms are presented. Next, 
typical examples of the experimental work are 
presented. Presentation and discussion of these 
experimental results complete the report. 


SCOLE APPARATUS 

The SCOLE hardware and support software is 
described in detail in 12]. In this report, only 
those elements of SCOLE used in this research are 
described. Referring to the schematic diagram, 
figure 2, SCOLE contains two major structural 
elements of interest: a planar, hexagonal, 
tubular structure representing an antenna 
reflector; and a single tubular flexible mast 
connecting the antenna to the platform. For this 
experiment, the platform is fixed to ground; only 
the mast and antenna portions are dealt with in 
this paper. The system actuators consist of three 
torque wheels that produce torque in three 
mutually orthogonal directions (see figures 1 and 
2 ). 


Experiments are run on SCOLE using a digital 
computer that has a UNIX-like operating system 
called UNOS. Programming is accomplished In a 
combination of C and FORTRAN 77 programs. The 
computer has analog-to-digital converters used for 
sampling the rate gyro data, digital-to-analog 
converters used to command the torque wheels, and 
a process timer which achieves precise internal 
timing of the data sampling process. This 
equipment has been added to the original system 
along with software drivers which can be evoked 
from either C or FORTRAN 77 programs. 


The output is of the form 
y = Hx 

For a rate sensor, H ■ [0 cl, c Is the mode 
slope at the sensor location. 

For digital control implementation, the 
control input Is assumed constant over the uniform 
sample time interval of T seconds, and the 
continuous-tlme model is converted to its 
discrete-time equivalent by integration of 
equation (2) over the Interval. Thus, the 
difference equation describing the motion appears 
as 


Vi " *V ru k (4) 

where * » e AT , T - X T e At dt B = (*-I)A _1 B since 
A is nonsingular. The matrix I is a 2x2 identity 
matrix. 

An autoregressive form of the discrete time 
model may be accomplished by taking the 
z-transform of the equation (4) and solving for 
the sampled sensor output, y^, in terms of the 

input actuator 


*k * Vk-l* Vk-2 + (b lVl + Vk-2 

where 

a l = *11 + *22 

®2 * * 12*21 " * 11*22 


♦ b^iL 0 )b (5a) 


b l “ r 22 


b 2 = *2i r i2 ' *li r 22 


SYSTEM MODELS 


The model we seek for SCOLE should 
incorporate the actual natural frequencies, 
damping factors, and control effectiveness 
coefficients of the system. To this end, each 
mode of the vibrational dynamics of SCOLE Is 
modeled as a single-input, single-output system 
described by 

x = Ax + Bu (1) 


wherein 


A = 


'0 1 
_ -u 2 -2<u> 


and 



( 2 ) 


x - modal state vector 
u - control input vector 
b - modal parameter of actuator location 
u - natural frequency of mode 
£ - damping ratio of mode 


b * c • b 

Taking the z transform of equation (5a) and 
sir ranging like terms (and letting the a 
coefficients absorb their respective signs) 
results in the output equation form 

(z 2 + az ♦ a z )Y(z) = (b^ + b 2 )U(z)b (5b) 

The above derivations are for the single mode 
case. A more general form of the sampled output 
equation, which Includes multiple modes, is 
presented. This case still deals with single 
Input, single output only. It can be shown that a 
decoupled multiple mode state space system can be 
written in block diagonal form, using the A 
matrix of equation (2) as the block elements for 
each mode. Because of the diagonal property of 
this system and the fact that the corresponding 
output matrix has the H matrix of equation (3) as 
elements, the z transformed multiple mode output 
equation can be written in the following compact 
form 


5+9 



N 

d(z) Y(z) = l b 
1 = 1 1 


, i 

(rijZ 


i " 

V n 

j=i * 


(z) U 


( 6 ) 


2 i I 

where d{z) = f[ d (z) , d (z) = z + a z+ a 

j = i J ^ c 

and N is the number of modes. The super- and 
subscripted quantities refer to the 1th and jth 
modes and the b^ terms are the modal control 

effecjiveness coefficients of the torque wheel. 

The a 1 coefficients are the AR coefficients of 

1,2 

the 1th mode and n can be directly computed 
1 1,2 

given a i . Thus, the model satisfies 


, 2 K 2N-1 o. 

< z + V + ... + * 2H z ) Y 


E C3, jZ 2 "' 1 * ft 
1=1 ' 


,. 2 zZ,, ' a+ ••• *A.»*° )b . U C7) 


where: , .... - coefficients of d(z) 


i» .... a - coefficients of the 
’ product of righthand 

side of equation (6) 

This is the model form used to represent the 
dynamics for each test conducted herein. 


PARAMETER IDENTIFICATION PROCEDURES 

Linear least-square estimation (LSE) is used 
to identify the difference equation (7) 
parameters. This method was selected because of 
its computation efficiency and implementation 
simplicity [7,8], The identification process is 
carried out for each mode and each actuator. A 
structural excitation test is conducted wherein 
the structure, 'initially at rest, is sinusoidally 
forced by a single actuator at the predicted 
natural frequency of the mode of interest. The 
actuator Is turned off when the output of the 
sensor of interest reaches a predetermined 
magnitude. The data recording Is continued to 
obtain free-decay data. 

The test data is processed in a two-step 
operation. For the first step, the AR 

coefficients (a ,a , ...) of the model are 
1 2 

identified using the free-decay portion of the 
data. For the second step, the identified AR 
coefficients are used to obtain the control 

effectiveness coefficients (b ,b , ...) with data 

1 2 

taken from the excitation portion. 


Step 1 - Identification of the Free Decay 
Coefficients 

First, the spectral content of the 
free-decay portion of the data is examined to 
determine the number of modes present in the data. 
The order of the model is taken to be twice the 
number of modes determined. Then, the free-decay 


data Is digitally filtered to suppress noise an 
the signal due to any modes not wanted in the 
model. Finally, the filtered data is processed 
using the standard least-squares estimation to 
identify the AR coefficients for the number of 
modes selected for modelling. 


The Identification of the AR (a 1 ) 

1,2 

coefficients for each mode generally depends on 
the data base used in the estimation. As more 
significant data is added, the estimates should 
converge to a value and the varisince of the 
estimates will improve to a limit, which depends 
on the measurement noise and the model. After 
this "convergence" occurs the mean and variance of 
the estimates should remain constant. Therefore, 
the variation of the estimates is examined as data 
is added to the data base and the mean of the 
estimates is taken over the last several data base 
additions. To ascertain confidence in the 
estimates, the variance of the estimates is also 
checked. 

The natural frequency and damping factor are 
obtained from a 1 and a^ by finding the roots 

z of Its char^cteristlc equation and using the 
relation z=e In the primary strip, where T 
is the uniform sampling period. The following 
equations for the natural frequency and damping 
factor are obtained for both the overdamped and 
underdamped cases: 


Overdamped 

2 

a^/4 + a^ > 0 

(8a) 

*1,2 " 

+ in [* \ 

(8b) 


“ * / s l s 2 

(8c) 


C « -(s 1 +s 2 )/2w 

(8d) 

Underdamped 

— aj/4 + a 2 < 0 

(9a) 


= 2T ln(_a 2 ) 

(9b) 

“d 

, . -1 / , 4a_ / 2 

= 1 tan / -1 - 2/a, 

(9c) 


u = y ♦ <r 2 

a 

(9d) 


C * - <r/tj) 

(9e) 
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Step 2. - Identification at £ Cbfttrgl. 

Effectiveness Parameters 

Once the values of the AR coefficients for 
the node of Interest are determined, a similar 
linear least-squares scheme used In [6] is 
employed to obtain the control effectiveness of 
the torque wheel actuator with respect to the mode 
of interest. 


of the testing for the mode 1 tests. Figure 4a 
shows the excitation portion of the mode 1 test 
with the x-axis (the 1-axls of figure 3) torque 
wheel. The data was sampled at 50 samples/second 
and the actuator commands were updated at this 
same frequency. The free-decay portion is shown 
in figure 4b. Although the time scales both start 
at 0, the free-decay portion Is, In fact, a direct 
extension of the excitation portion. 


From equation (7), an error equation can be 
defined as 


2M n 

e * * E * y k ~ E W. .u + ■ • 

k " i-i 

+ B u )b 

,2N k-2M 1 


( 10 ) 


To calculate the b ’s by least squares, the 

parameter J , defined in equation (11), is 
minimized with respect to the b^s. The J 

parameter represents the total error for all the 
samples taken. 

J- 1/2 JeJ HD 

k-1 * 


s - total number of k samples 

The resulting least square equation is In matrix 
form, shown in equation (12). 

V = P* l 0 (12) 


Step 1 — The free-decay portion of the data run 
Is first analyzed to determine the number of modes 
present in the data. Figure 5a Is a Fast Fourier 
transform of the free-decay data. The transform 
is only plotted to the Nyquist frequency. 
Experience has shown that the accommodation of 
noise in a finite data base does not allow one to 
estimate to the Nyquist frequency but to about one 
decade down. In this case about 5 Hz. Since the 
data includes a significant peak at around 10 Hz. 
It should be filtered to remove this mode and the 
noise in the signal. The sensor data was 
digitally filtered using an 8th order Butterworth 
filter with the cutoff frequency at 3 Hz. This 
was taken to provide significant attenuation of 
the 10 Hz peak In the unfiltered spectrum. 

Figure 5b shows the spectral content of the 
filtered data. Indeed, the 10 Hz. mode has been 
removed and three modes remain. The separation of 
the magnitudes of these remaining modes is a 
factor of 200. This means that extraction of the 
coefficients of the AR model for all three modes 
will be difficult with a limited data base. 

Hence, we select a two mode model for the purpose 
of estimating the dynamics of the free-decay data. 
The ARMA model Is thus of order 4. 


V - Nxl vector, with b^s as elements 

p - NxN matrix, with element P 

s r 2M 2M -j 

-E I E p ,..V. 0 E.VA-. 

k-1 L ■-! 


0 - Nxl vector, with element Q 


S r 2H 

■ E K ♦ E « 

k=l L »-l 


y ) 

wT k-m 


l 

2K 1 

\m 1 J 


For the case of a single mode, the above equation 
should be the same as the least square equation in 
[61. 


RESULTS AND DISCUSSIONS 

A finite-element analysis of SC0LE has been 
conducted to calculate the modes and frequencies 
of Interest. Figure 3 shows the first vibration 
mode shape and frequency and that of the third 
mode. The first mode Is essentially bending. The 
second mode is also a bending mode and looks 
similar to mode 1 except that the bending is 
orthogonal to that of mode 1. The third mode Is a 
torsional rotation of the reflector about the 
mast. The remaining modes are higher In frequency 
and appear similar to the classical Euler bending 
modes of a flexible beam. The predicted natural 
frequency of the first mode Is .43 Hz. This is 
the frequency used during the excitation portion 


Least square estimation is used to determine 
the AR parameters from the data. The variation of 
the coefficients as data Is added to the data 

base for estimation Is shown In figure 6. The 
abscissa is interpreted as follows: the last point 
(100) corresponds to the estimate with the entire 
free-decay data base (this Includes 1465 data 
points sampled at 50 samples/sec.) and the first 
point (0) corresponds to a data base of 1365 
points, approximately 27 seconds. This variation 
is typical and the other two a parameters have a 
similar variation. The mean and variance of the 
parameters are: 


*1 = 

3.8403 

var (aij) 

°2 * 

-5.6792 

var(a 2 ) 

*3 _ 

3.8363 

var() 

®4 “ 

-0.9980 

var(a^) 


8.567(10) 
6.787(10) 
6 . 220 ( 10 ) 
6.650(10) 


The characteristic polynomial corresponding to 
these coefficients has factors which correspond to 
mode natural frequencies and damping ratios of 
(see equations 8 and 9): 

» .4610 Hz ^ = .0011 

« - 3.1479 Hz C = .0024 

2 2 




Step 2 — The data for the frequency and damping 
ratio obtained in step 1 has been used in the 
algorithm of [63. Examination of the excitation 
portion of the test data (figure 4b), reveals the 
presence of a disturbance induced when the torque 
wheel is started. It is believed that this is the 
result of start-up friction in the torque wheel 
assembly. This disturbance disappears in about 7 
seconds. Hence, this part of the data was not 
included in the data base for estimation. The 
variation of the least square control 
effectiveness estimate as the last 400 samples are 
added to the data base is shown in figure 7. For 
the lowest number of points in the data base (0 on 
the abscissa) the estimate is seen to not have 
converged and more data is needed. As points are 
added the estimate statistics converge. It Is 
important to use the mean taken over a converged 
portion of figure 7 and not a single point. The 
value of the control effectiveness mean taken over 
the 400 data points of figure 7 is .0171 and the 
variance is 1.0576(10) . 
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This paper describes the linear least-square 
Identification procedure used to obtain an 
empirical model of the vibrational dynamics of 
SCOLE. Empirical data is obtained autonomously by 
exciting the structure from a quiescent state with 
torque wheels and observing the time history data 
of rate gyro sensors. The torque wheels are then 
shut down and free-decay data is recorded. The 
free-decay portion of the data is analyzed using 
the digital Fourier transform to determine the 
best model order to use in modelling the response. 
Linear least-square analysis is then used to 
select the parameters that best fit the output of 
an autoregressive (AR) model to the data. The 
control effectiveness of the torque wheels is then 
determined using the excitation portion of the 
test data, again using linear least squares. 

Typical experimental data are presented that 
reflect the performance of the identification 
algorithms. 
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Figure 2.- Schematic of 
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Figure 3.- First bending and torsional node 
shapes and frequencies. 




Figure 4.- Typical torque wheel test data 
including the excitation and 
free-decay portions. 
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Figure 5.* Fast Fourier transform magnitude plots 
of the rate gyro output during the 
free-decay portion. 
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Figure 7.- The control effectiveness parameter 

for the first mode using an average of 
the AR parameters of a single mode 
model. 
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ABSTRACT 

Future NASA missions will depend on 
significantly large, flexible, and complicated struc¬ 
tures in outer space. These structures will require 
very stringent pointing and vibration suppression 
requirements. The design of the active controllers to 
achieve these objectives requires accurate knowledge 
of the dynamic behavior of these structures. The 
structures are lightly damped distributed parameter 
systems and have an infinite number of frequencies of 
vibration. This paper presents a system identifcation 
technique that represents the structural dynamics in 
terms of a reduced-order modal model. The finite- 
element method is used to reduce the order of the 
physical model, and then the model is transformed 
from a physical to a modal form. The maximum 
likelihood method is used to parameterize the modal 
model, based on experimental observations of the 
structural dynamics. 


NOMENCLATURE 

gSt ) Generalized Coordinate Vector 
u_(t ) Input Vector 

x(0 Predicted Output 

M Mass Matrix 

K Stiffness Matrix 

B p Input Influence Matrix 

C p Output Influence Matrix 

O Eigenvalue Matrix 
<J> Eigenfunction Matrix 

co Modal Frequency 

3 i(r) Modal Displacement Vector 

iL(f ) Modal Acceleration Vector 

0(f ) Vector of System Parameters 

Damping Ratio 
7(0) Cost Function 

T Period of Oscillation 

£(* ) Predicted Acceleration Vector 

zSs ) Observed Acceleration Vector 

G Covariance Matrix of Noise 

p£z ) Position Vector of Instrument on Mast 


L Length of the Mast 

p Density of the Mast 

A m Cross-sectional Area of the Mast 
E Modulus of Elasticity of the Mast 
I Area Moment of Inertia of the Mast 

G ^ Shear Modulus of the Mast 

M x Mass of the Shuttle End of the Mast 

M 2 Mass of the Reflector End of the Mast 

J i Mass Moment of Inertia of the Shuttle Mass 

J 2 Mass Moment of Inertia of the Reflector Mass 

l ^ Polar Moment of Inertia of the Mast 

W (©) Jacobian Matrix 
H (0) Hessian Matrix 

INTRODUCTION 

Many future NASA missions would utilize 
significantly large, flexible and very complicated 
structures in outer space* These large space structures 
would require very stringent pointing and vibration 
suppression requirements. The active controller that 
can acheive these objectives will have to be designed 
with very accurate knowledge of the dynamic 
behavior of the structure to ensure performance 
robustness to a variety of disturbances and uncertain¬ 
ties. It is recognized by control'engineers that there 
are certain inherent problems in the design of active 
controllers for this class of large flexible space struc¬ 
tures. A space structure is a very lightly damped 
distributed parameter system and hence has an 
infinite number of frequencies of vibration. A large 
number of these frequencies have extremely low 
values and they are very closely spaced together. 
Additionally, numerical solutions of the higher fre¬ 
quencies tend to yield large errors in their calcula¬ 
tions. Thus, it is a challenge to design a state-of-the- 
art multivariable control system for a lower order 
model of a large flexible space structure. 

Several techniques for designing control systems 
for vibration suppresion for flexible spacecraft have 
been proposed, and various computer and experimen¬ 
tal studies have been conducted to validate these con¬ 
trol schemes. An experimental test article was 
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designed under the cognizance of the Spacecraft Con¬ 
trol Branch at NASA Langley Research Center for 
testing additional methods and for directly comparing 
the various competing techniques on a common exper¬ 
imental model (SCOLE). 

This test article consists of a softly supported 
dynamic model of a flexible beam and reflector gril¬ 
lage attached to the space shuttle. The control objec¬ 
tive includes the task, of directing the line-of-sight 
(LOS) of the antenna-like configuration toward a 
fixed target, under conditons of limited control 
authority and random disturbances. 

IDENTIFICATION OF SCOLE ASSEMBLY 

The first task in the control system design is that 
of system identification. The dynamic behavior of the 
experimental model must be described analytically. 
The finite-element method can be used as a structural 
analysis method to define the reduced-order model 
for the infinite-dimensional SCOLE test assembly. 
The finite- element method represents the continuous 
body as a finite number of discrete basic elements 
that each possess mass and stiffness properties that 
can be assembled to form mass and stiffness matrices. 

System identification combines experimental 
observations and theoretical predictions to form an 
accurate model of an unknown system. The observa¬ 
tions are used to quantify the unknown model 
parameters and to minimize the error between the 
observed system dynamics and the model-predicted 
dynamics. 

The SCOLE test assembly can be modeled as a 
free-free beam with rigid end masses. The objective 
of the identification process is to describe the dynam¬ 
ics of the antenna and mast with respect to the shut¬ 
tle body. The observations that will be used to 
parameterize the system are the linear and angular 
acceleration rates of the mast and antenna. The 
available data consists of two sets of measurements 
of mast accelerations in the x and y planes, one set of 
measurements of the x and y accelerations of the 
reflector, and the angular acceleration about x, y, and 
z measured at the junction of the mast and antenna. 
Only the motion of the antenna and mast with 
respect to the shuttle body is of interest so that the 
pitch and roll of the antenna adequately describes the 
motion in the z-direction. 

The data set thus describes the system as one 
that can be represented by ten degrees-of-freedom. 
The frame of reference is chosen so that the origin is 
at the universal joint suspension point of the test 
article with the z-axis running along the center of the 
mast. The desired form of the model is one that 
transforms the physical properties of mass and 
stiffness into the modal properties of natural fre¬ 
quency and damping ratio. The ten degrees-of- 
freedom of the physical model results in a truncated 
modal model with five modes. 


BASIC ANALYTICS 

The structural damping of the assembly is very 
light and is neglected in the physical model. The 
model can thus be written as 

MUt ) + KgJj ) = B p u(t ) (1) 

x(t) = C p aSf ) (2) 

where gS. 0) = i_(0) = g# are the initial conditions. 

Let 0= diagonal jo^ 2 <o 2 2 ...o> n 2 j be the 

eigenvalue matrix corresponding to the physical 
model and be the corresponding eigenfunction 
matrix of the eigenvalue problem. 

It can be shown that 

t>. (3) 

Now normalize $ so that 

& M <*> = / (4) 

and let 

gSf ) = <*a(f )• (5) 

This results in the modal transformation 

£.(* ) + ) (6) 

X(0 = 0*21(0 (7) 

where B m = <$> T B p and C m — C p and the initial 

conditions are = Qr'g#, = 3r l g. The modal 
equation is equivalent to the physical equation in the 
sense that it produces the same input-output descrip¬ 
tion of the structure. The viscous damping can be 
modeled in terms of , the damping ratio of the ith 
mode. Including the damping and reordering the 
state variables, the dynamic model can be written as 
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where fj_= o> = natural frequency vector and 
7j,= damping ratio vector. 
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where the subscript j indicates the sensor locations 
along the z axis. 

The physical properties of the test assembly are 
used to evaluate the output vector £(f ) as follows 
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where ) is the vector of unknown parameters to 
be adjusted. 


ALGORITHM DEVELOPMENT 
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The coefficients A, B, C, D, E, and F are evaluated 
from the boundary conditions and C x and C 2 are 
evaluated from the initial conditions. The first 
coefficient of the spatial parameters (A, E) is set equal 
to one and the resulting system is reduced so that the 
following system of equations is obtained 


The system identification algorithm uses the 
maximum likelihood method to improve the accuracy 
of the finite-element model. The model is assummed 
to be of the same form as the true model and con¬ 
verges to the true value as the parameters converge to 
their actual values. The cost functional , 7(0), is 
minimized so that the error betwen the model- 
predicted output, y_( f )» anc * observed output, 
z£t ), approaches zero. The algorithm uses the 
modified Newton-Raphson method to minimize the 
error. 

The cost functional, 7(0) can be expressed as 
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The first measurements in each data set , z_(t x ), are 
used to evaluate the initiaJL modal conditions and the 
coefficients of the modal parameter as follows 
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The cost function, 7(0), is minimized by adjust¬ 
ing the vector of system parameters, Q, by an itera¬ 
tive method in which the new parameter estimate, 
+ lf is computed on the basis of the present esti¬ 
mate Qn • 
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APPLICATION TO EXPERIMENTAL DATA 

The algorithm is used to evaluate eight sets of data 
and each set represents different initial conditions in 
order to excite vibrations corresponding to roll, pitch, 
and yaw motions. The system identfication is carried 
out for each data set to yeld the modal model. Each 
data set is a sequence of 1200 observations that were 
sampled at a 0.025 s sampling rate. The parameters 
are evaluated for each observation and then the 
resulting set of parameters are evaluated using sta¬ 
tistical analysis to define the dominant natural 
modes. 

CONCLUSIONS 

The identification technique presented in this paper 
has been found to be very effective in modeling the 
SCOLE test article and for active controller designs. 
The method has some convergence problems but has 
the potential to be extremely useful for solving the 
system identification problem of other large flexible 
spacecrafts. 
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SUMMARY 


The Spacecraft Control Laboratory Experiment Is a facility for the 
investigation of control techniques for large flexible spacecraft. The control 
problems to be studied are slewing maneuvers and pointing operations. The facility 
implements the salient characteristics of a flexible satellite with distributed 
sensors and actuators. 

The flexible satellite is represented by a continuous structure consisting of a 
large mass and inertia connected to a small mass and inertia by a slender, flexible 
beam. The structure is suspended by a single cable mounted to a universal joint at 
the system C. G. The sensors for the experiment consist of aircraft quality rate 
sensors and servo-accelerometers. The shuttle attitude will be determined through a 
combination of inertial measurements and optical sensing techniques. Actuators for 
the experiment consist of Control Moment Gyros, reaction wheels, and cold gas 
thrusters. Computational facilities consist of micro-computer-based central 
processing units with appropriate analog interfaces for implementation of the 
primary control system, the attitude estimation algorithm and the CMG steering law. 
Details of the experimental apparatus and the system software are presented in this 
paper. 


ABSTRACT 

A laboratory facility for the study of control laws for large flexible 
spacecraft is described in the following paper. The facility fulfills the 
requirements of the Spacecraft Control Laboratory Experiment (SCOLE) design 
challenge for a laboratory experiment, which will allow slew maneuvers and pointing 
operations. The structural apparatus is described in detail sufficient for 
modelling purposes. The sensor and actuator types and characteristics are described 
so that identification and control algorithms may be designed. The control 
implementation computer and real-time subroutines are also described. 

INTRODUCTION 

A modelling and control design challenge for flexible space structures has been 
presented to the technical community by the NASA and IEEE (ref. 1). The Spacecraft 
Control Laboratory Experiment (SCOLE) was constructed to provide a physical test bed 
for the investigation and validation techniques developed in response to the design 
challenge. The control problems to be studied are slewing maneuvers and pointing 
operations. The slew is defined as minimum time maneuver to bring the antenna 
line-of-sight (LOS) pointing to within an error limit of the pointing target. The 
second control objective is to rotate about the line of sight and stabilize about 
the new attitude while keeping the LOS error within the bound 5. The SCOLE problem 
is defined as two^design challenges. The first challenge is to design control laws, 
using a given set of sensors and actuators, for a mathematical model of a large 
antenna attached to the space shuttle by a long flexible mast. The second challenge 
is to design and Implement the control laws on a structural model of the system in a 
laboratory environment. This report gives preliminary specifications of the 
laboratory apparatus so that Interested investigators may begin design and 
simulation for the laboratory experiment. 

The laboratory experiment shown in figure 1 attempts to Implement the 
definition of the modelling and control design challenge within reasonable limits of 
the 1-g atmospheric environment. The experimental facility exhibits the essential 
SCOLE characteristics of a large mass/inertia (space shuttle model) connected to a 




small mass/inertia (antennae reflector) by a flexible beam. Control sensors and 
actuators are typical of those which the control designer would have to deal with on 
an actual spacecraft. Some trades are made in terms of structure, sensors, 
actuators, and computational capability in order to develop the experiment in a 
timely and cost-effective manner. To this end, the basic structure is made of 
homogeneous, continuous elements. It is suspended from a steel cable with the 
positive z-axis of the shuttle pointing up, thus minimizing the static bending of 
the antenna mast. The suspension point Is a two-degree-of-freedora gimbal for pitch 
and roll with yaw freedom supplied by the suspension cable. The sensors are 
aircraft quality rate sensors and servo-accelerometers. The shuttle attitude will 
be determined through a combination of inertial measurements and optical sensing 
techniques. 

The shuttle control moments are provided by a pair of two-axis control moment 
gyros (CMG's). Mast-mounted control torques can be applied by a pair of two-axis 
reaction wheels. The reflector-based forces are provided by solenoid-actuated cold- 
air thrusters. Reflector mounted torque devices are a trio of high-authority 
reaction wheels. Computational facilities consist of micro-computer-based central 
processing units with appropriate analog interfaces for Implementation of the 
primary control system, the attitude estimation algorithm, and the CMG steering 
law. All of the elements which make up the SCOLE experiment are described in detail 
In the following text. 

The description of the apparatus covers five major groups: The basic 
structural elements are described and pertinent dimensions and structural properties 
are provided. The sensor locations and their dynamic properties are presented. The 
actuator locations and estimated dynamic properties are also given. The mass prop¬ 
erties of the combined structure, sensor and actuator system are given. Finally, 
the computing system and analog Interfaces are described. 

The contents of this report are considered accurate at the time of 
publication. All of the planned SCOLE components are Implemented and are available 
to the user at a raw signal level. However, due to continued refinement of some of 
the components, specific details of the system may change over the life-time of the 
experimental apparatus. 


STRUCTURES 

The SCOLE is comprised of three basic structures, the shuttle, the mast, and 
the reflector panel. The assembly of these individual components and the global 
reference frame are shown In figure 2. 

The shuttle planform is made from a 13/16-Inch steel plate and has overall 
dimensions of 83.8 by 54.0 inches. Its total weight Is 501.7 pounds. The shuttle's 
center-of-mass is located 3.4 inches below the experiment's point of suspension, and 
26.8 inches forward of the tail edge (fig. 3). 

The mast is 120 Inches long. It is made from stainless steel tubing and weighs 
4.48 pounds. One-inch thick manifolds are mounted to the mast at each end. The 
assembly of these parts and their dimensions are shown in figure 4. 

The reflector panel is hexagonal In shape, made from welded aluminum tubing, 
and weighs 4.76 pounds (fig. 5). It is located 126.6 inches below the SCOLE's point 
of suspension. The center of the reflector is located at 12.0 inches in the x 
direction and 20.8 inches in the y direction from the end of the mast. 
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The complete system Is suspended from an 11-foot cable attached at the system 
center-of-gravity via a universal joint. Roll and pitch rotational freedom Is 
provided by pillow-block ball bearings which have an estimated break-out torque of 
0.1 ft-lb. The universal joint is shown in figure 6. It is fixed to the shuttle 
plate, and the system center-of-gravity is made to coincide with the center-of- 
rotation by means of an adjustable counter balance system. 


SENSORS 

The sensors for the experiment consists of nine servo-accelerometers and two, 
3-axis rotational rate sensing units. An optical sensor will provide yaw attitude 
of the shuttle. The power supplies for these sensors are mounted on the shuttle 
plate to minimize the number of large gauge wires which must cross the universal 
joint suspension point. Only a single 115 VAC cable and 33 signal wires cross the 
universal joint. The wires for the sensors are routed on the shuttle and along the 
mast • 


Accelerometers 

All nine accelerometers have a frequency response which is nearly flat up to 
350 Hz. Linearity is within 0.17 percent of the full-scale output. A typical 
calibration is presented in figure 7. Individual calibrations are available on 
request* 

The shuttle-mounted accelerometers shown in figure 8a sense the x, y, and z 
accelerations. These sensors are distributed away from the suspension point to aid 
inertial attitude estimation. The locations and sensitive axis are shown In 
figure 8b. 

The mast-mounted accelerometers shown in figure 9a sense x and y acceleration 
at locations about one-third of the mast length from each end. The positions and 
sensing axis of the devices are shown In figure 9b. 

The reflector-mounted accelerometers are shown in figure 10a. They are 
positioned in the center of the reflector below the thrusters and sense the x and y 
accelerations. The coordinates and sensing axis of the devices are shown In 
figure 10b. 


Rate Sensors 

The rotational rate sensors are three-axis, aircraft-quality Instruments. The 
frequency response is approximately flat to 1 Hz and -6 db at 10 Hz. Linearity Is 
about 0.6 percent full scale. A typical calibration is shown in figure 11. The 
range Is 60 deg/sec for the yaw and pitch axis and 360 deg/sec. for roll. The 
threshold is 0.01 deg/sec. 


The shuttle-mounted rate sensor package, shown In figure 12a, senses 
three-axis, rigid body angular rates of the shuttle plate. Its coordinates and 
sensing axis are presented In figure 12b. 

The mast-mounted rate sensor package, shown in figure 13a, senses three-axis 
angular rates at the reflector end of the mast. Its coordinates and sensing axis 
are presented in figure 13b. 
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The sensor information required for control system design is summarized in 
Table I. The sensor type is listed in column 2, and its sensed variable is listed 
in column 3. The analog interface channel is listed next. The coordinates of the 
parts with respect to the universal joint are listed in the next three columns. The 
sensitivities in terras of analog-to-digital converter units are listed next. The 
error list shows the RMS deviation of the rate sensors or the percentage of full- 
scale linearity error for the accelerometers. The linear range of the instruments 
is listed in the next-to-last column. 


OPTICAL SENSOR 

An optical sensor will be provided to determine yaw attitude of the shuttle. 
The optical sensor is a planar photo-diode with appropriate optics mounted on the 
ground. The outputs of the sensor are proportional to the position of an infared 
light source on the shuttle. The sensor data is processed by a dedicated micro¬ 
controller and then sent to the main CPU over a serial data link for transformation 
to attitude angles. No photographs or calibration data are available for this 
device. 


ACTUATORS 

The actuators consist of both proportional and on-off controllers. Shuttle 
attitude control is provided by a pair of two-axis control moment gyros (CMG's). 

Mast vibration suppression can be achieved with a pair of orthogonally mounted reac¬ 
tion wheel actuators positioned at two stations on the mast. Reflector forces are 
provided by four cold gas jets. Reflector torques are provided by three ortho¬ 
gonally mounted reaction wheels at the end of the mast. As with the sensors, all 
devices are inertial, and the power supplies and amplifiers are mounted on the 
shuttle. Fifteen command signal wires cross the universal joint. All actuators 
were manufactured in house. 


Control Moment Gyros 

The CMG's each have two gimbals which are equipped with individual direct drive 
DC torque motors. The momentum wheel is mounted in the inner gimbal and driven by 
two permanent magnet DC motors. The nominal operational momentum is about 2.5 ft— 
lb-sec. The gimbal torque motors are driven by current amplifiers so the output 
torque will be proportional to the command voltage sent to the amplifier. The 
gimbal torquers will produce +/— 1.5 ft—lbs at frequencies up to 1kHz. The gimbals 
are instrumented with tachometers and sine—cosine potenlometers to facilitate 
decoupled control of the shuttle attitude angles. A dedicated computer will be used 
to control the CMG gimbals. Routines will be provided so that users may command 
decoupled shuttle torques or gimbal torque commands. 

The sensitivity calibration curve of a typical gimbal motor is shown in 
figure 14. No other calibration data are available for the CMG's. 

The forward CMG is shown in figure 15a. Note that the outer gimbal is fixed 
and parallel to the pitch axis of the shuttle. The inner gimbal is nominally orien¬ 
ted so that the spin axis of the momentum wheel is parallel to the shuttle z-axis. 
The second CMG is mounted at the rear of the shuttle so that the outer gimbal is 
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parallel to the z-axis. The inner gimbal is nominally oriented so that the rotor 
spin axis is parallel to the shuttle x-axis. The coordinates and nominal axis of 
actuation of the CMG's are shown in figure 15b. 


Reaction Wheels 

The mast-mounted reaction wheels consist of aluminum disks with inertia of 

about 0.00027 lb-ft-sec 2 mounted directly on the drive shaft of a 20 oz-in 
permanent magnet DC motor. The motors are powered by high bandwidth current 
amplifiers. A torque sensitivity plot is presented in figure 16. No other 
calibration data are available. A typical reaction wheel assembly is shown in 
figure 17a. The two actuator locations and their axis of actuation are shown in 
figure 17b# 

The mast end mounted reaction wheels consist of D C permanent magnet pancake 
motors which are mounted with the armature fixed to the structure. The stator and 
case of the motor are allowed to rotate via a slip ring assembly, thus providing 
high inertia mass to fixed mass efficiency. 

The motors are powered by high bandwidth current amplifiers* The torque 
capability of these devices is estimated to be about 50 oz-in. No sensity plot or 
other calibration data is presently available. The three-axis reaction wheel 
assembly is shown in figure 18a. The actuator locations and their axis of actuation 
are shown in figure 18a. 


Thrusters 

The control forces on the reflector are provided by solenoid actuated cold gas 
jets. The thrusters are mounted in the center of the reflector and act in the x-y 
plane. The jets are supplied by a compressed air tank mounted on the shuttle. The 
pressurized air travels through the mast to the solenoid manifold, which gates the 
air flow between the regulated supply tank and the thrusters as shown in figure 19. 
Thrust is Initiated by opening the solenoid with a discrete command. The rise time 
and transient oscillation of thrust is shown in figure 20. The magnitude and dura¬ 
tion of the thrust before the air supply is depleted at 60-psi nozzle pressure is 
shown in figure 21. The pertinent data from figures 20 and 21 are tabulated in 
Table II. 

The thrusters are shown in figure 22a. Their location and axis of actuation 
are shown in figure 22b. 

The actuator information required for control system design is summarized in 
Table III. The actuator type and direction of action are listed In Column 2. The 
analog interface channels are listed in Column 3. The coordinates of the devices 
are listed in the next three columns. The sensitivities of the actuators in terms 
of digital-to-analog converter units are shown in Column 7. An estimate of the 
thruster RMS deviations exhibited in figure 19 is presented as error data. No other 
error data are available. The maximum range of the system actuators is shown in the 
next-to-last column. 
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MASS PROPERTIES 


The position and weight of the various pieces of equipment, which collectively 
form the SCOLE apparatus, are cataloged in Table IV. Distances are measured from 
the point of suspension to the approximate center-of-mass of each component. Each 
major component is listed in the second column of the table. The x, y, z 
coordinates are listed next. The weight of each component is listed in Column 6. 
The remaining columns are the mass moments and moments of inertia. The totals for 
the complete system are presented on the bottom row. 


COMPUTER SYSTEM 

The main computer for control law implementation will be a micro - computer based 
on the Motorola M68000 microprocessor. The computer has 2.0 M-byte of random access 
memory and a 40 M-byte hard disk. The operating system is based on UNIX with C, 
Fortran and Pascal compilers available for applications programming. The computer 
has 12 serial ports and 1 parallel port. Terminals are connected on four of the 
ports and an answer-only modem Is attached to another. One port Is used for an 
originate-only modem. A line printer is attached to another port. The optical 
sensor is connected to a serial port. The IBM PC, which Is used to drive the CMC s. 
Is also connected to a serial port. 

Analog interfaces consist of a four-bit, output-only discrete channel: an 8 
bit discrete output port, an 8-bit discrete Input port, 8 digital-to-analog 
converters, and 64 analog-to-digital converters. All converters are 12-bit devices 
with a range of +/-10v. These Interfaces are shown schematically in figure 23. The 
CMG control software required for the PC should be relatively transparent to the 
controls designer who will be operating on the CRDS computer. 

Subroutines for accessing the analog interfaces and setting the digital 
sampling interval are described in Appendix A. The most commonly used routines are 
listed below. 

For accessing the analog devices: 

getadc - read the analog-to-digital converters 
setdac - set the digital-to-analog converters 
thrust - set the cold-gas thrusters. 

To control the sampling interval: 

rtime - sets the sample period marks the beginning 
of a real-time loop. 

A time-line of the synchronization of the sample Interval using the routine 
rtime and the analog interface routine usage is shown in figure 24. The basic 
operation is as follows: 

The user first calls rtime with flag-.true., and a valid sample Interval. 
After setting the timer period, the routine starts the user's real-time 
routine and the interval clock. The user routine will use some or all of 
the subroutine calls shown. When the user computations and actuator 
commands are complete, the routine must return to the top of the real-time 
loop and once again call rtime. If this occurs before the end of the 
sample Interval, the time-out condition will be inhibited and rtime will 
wait for the next rising edge of the sample interval clock and then return 
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to the calling program. If the user computations take longer than the 
sample interval, a time-out condition will be signaled to the operator 
when rtirae is called. The user may choose to ignore the condition and 
continue or may take specific steps to alleviate the condition. 

The procedure for logging on to the computer Is as follows: 

Set communication parameters to 1200 baud, 7 bit, even parity. 

Dial in to the Langley data communication switching system at 
804-865-4037. When connected, type a carriage return. 

To the system prompt "ENTER RESOURCE CODE" type "acrl." 

Wait until "GO" and the "name:" prompt appear on the screen. 

Type In your log-in information. All investigators will be given a three- 
letter log-in name (usually the university affiliation). 

Some useful system commands are listed below: 

To transfer a file to the experiment computer from a smart terminal, type 
cat) filename <cr> and then enable the upload function of the local 
terminal. When the upload is complete, type <cntl>d. 

To transfer a file from the experiment computer to a smart terminal, type 
cat filename then enable the download function of the local terminal and 
type a <cr>. 

To list the contents of a directory, type 1 <cr>. 

To look at a file, type p filename* 

To compile a FORTRAN program and link with real-time system commands and 
TCS graphics, type frt filename* 

Note: filename must have the extension .for. The executable code will be 

under filename without the .for extension. To run, simply type filename 
without any extension. 

To list the system commands, type 1 /bin. to get a description of any 
command, type describe cownd. 

System user guides will be available upon request from the Spacecraft Control 
Branch, M/S 161, NASA Langley Research Center, Hampton, VA, 23665-5225. Other 
details for operating in a real-time mode will be provided at the time of 
implementation. 


CONCLUDING REMARKS 


The SCOLE laboratory facility is an experimental apparatus which permits 
ground-based Investigation of identification and control algorithms for large space 
structures. The facility exhibits structural dynamics similar to those expected on 
the large satellites. The sensors and actuators are typical of those, which may be 
used on an operational satellite. The computational system Is reasonably sized with 
current technology processors and permits ready access to the facility for 
interested Investigators. 
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The description of the structural assembly, the sensor and actuator 
configuration, and software provided in this paper should be sufficient for SCOLE 
Investigators to begin designing identification and control algorithms for the SCOLE 
facility. 
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APPENDIX 


REAL-TIME SYSTEM SUBROUTINES 


Analog I/O system command. 


NAME: 

getadc Samples the analog-to-digital converters. 


COMMAND: 

getadc [-s] 


DESCRIPTION: 

This command is used to sample the analog-to-digital converters and display selected 

channels at the terminal. The +/- 10.0 volt input range is scaled to +/- 1.0 units 

so a single bit is worth .00049 units. The channels to be displayed are selected 
with the -s option. This option displays a menu which allows the user to set the 
print flags for the Individual ADC channels. Specific choices are: 

1) turn on all print flags, 

2) turn off all print flags, 

3) turn on a range of print flags, 

4) turn off a range of print flags, 

5) display current print flags, 

6) save current print flags. 

If the command Is executed without the -s option, the last set of print flags Is 
used to selectively display the ADC channels. 


USES: 

getadc.flags 


DIRECTORY: 

/bin 


SOURCE: 


/usr/csc/ele/getadc.c 
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Analog I/O system command. 


NAME: 

getadc Samples the analog-to-dlgital converters. 


COMMAND: 

getadc [-s] 


DESCRIPTION: 

This command is used to sample the analog-to-digital converters and display selected 

channels at the terminal. The +/- 10.0 volt input range Is scaled to +/- 1.0 units 

so a single bit is worth .00049 units. The channels to be displayed are selected 
with the -s option. This option displays a menu which allows the user to set the 
print flags for the individual ADC channels. Specific choices are: 

1) turn on all print flags, 

2) turn off all print flags, 

3) turn on a range of print flags, 

4) turn off a range of print flags, 

5) display current print flags, 

6) save current print flags. 

If the command is executed without the -8 option, the last set of print flags is 
used to selectively display the ADC channels. 


USES: 

getadc.flags 


DIRECTORY: 

/bin 


SOURCE: 

/usr/csc/ele/getadc.c 
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Analog I/O system command 


NAME: 

setdac Sets the digital-to-analog converters. 


COMMAND: 

setdac [—s] [-0] 


DESCRIPTION: 


This command is used to set the voltage on a range of digital-to-analog output 
channels. The +/- 1.0 unit output range is scaled to +/- 10.0 volts so a single 
unit is worth .0049 volts. The channels to be set are selected by executing the 
command with the -s option. This option displays a request for the range of 
channels to be set, and then queries for individual channel values in terms of 
units. All DAC channels may be set to zero by executing the command with the ”-0” 
option. 


USES: 

Nothing. 


DIRECTORY: 

/bin 


SOURCE: 

/usr/csc/ele/setdac.c 
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UNOS system command 


NAME: 


lterm Terminal emulator. 


COMMAND: 

term [—s] [-S] 


DESCRIPTION: 

This command connects the user terminal to the Langley central data communication 
switch at 1200 baud. This is a dumb terminal emulator which provides rudimentary 
file transfer capabilities. No attempt is made to emulate control codes of any 
particular terminal for editing purposes. 

The emulator commands are as follow: 

@ Return to UNOS (operating system.) 

! To download a file. 

“ To upload a file. 

To excute a system command. 

? For help. 

The options are: 

-8 300 baud 

-S 900 baud 

Upload means to transfer a file from the Charles River computer to the remote 
computer. The remote computer must have some mechanism for receiving the text. 

Download means to transfer a file from the remote computer to the Charles River 
computer. No attempt is made to check for existence of the receiving file name 
before saving the downloaded file. 


USES: 

/doc/cmds/1term.help 


DIRECTORY: 

/bin 


SOURCE: 

/jpw/lterm.c 
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Analog I/O system command 


NAME: 

8751_test Test the serial communication link to the 8751 boards. 


COMMAND: 

8751 test 


DESCRIPTION: j 

This command facilitates verification and calibration of the 8751 micro-controller ' 

interface over the RS—232 serial ports. The command queries for voltages to be 
output by the digital-to-analog converters on the 8751 boards. The data input is in 
terms of units with 2047 equal to 9.9951 volts and -2048 equal to -10.0000 volts. 

If the Input line contains only one value, all active boards are sent that value. 
Otherwise, individual values are sent. 

I 


USES: 

raotintQ 
motsub() 


DIRECTORY: 

/bin 


SOURCE: 

/u8r/rdb/8751 test.for 
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Analog I/O system subroutine 


NAME: 

getadcO ( C callable ) Sample a range of analog-to-dlgital converters. 


CALL: 

Int error, first_adc, last_adc; 
float adc_data_polnter; 

Int getadc( flrst_adc, last_adc, &adc_data_polnter) 
error * getadc(flrst_adc, last_adc, &adc_data_polnter) 


DESCRIPTION: 

This subroutine samples a range of analog-to-digltal converters. The +/- 10.0 volt 
Input range is scaled to +/- 1.0 units so a single bit Is worth .00049 units. The 
arguments are: 

first adc (int) First converter to be sampled (numbering starts from 

~ zero.) 

last adc (Int) Last converter to be sampled (maximum is 63.) 

&adc_data_pointer (*) Starting location for storing sample data. Data are 

floating point values with a range of +/- 1.0. 


RETURNS: 

error - 0 indicates valid transfer, 
error » Indicates bad range. 


USES: 

Nothing 


LIBRARY: None 


SOURCE: 
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/u8r/csc/ele/getadc_c.c 



Analog I/O system subroutine 


NAME: 

setdacQ ( C callable ) Set a range of digital-to-analog conveters. 


CALL: 

Int error, first dac, last dac; 
float dac data pointer; 

int setdac( first dac, last dac, & dac data pointer ) 


DESCRIPTION: 

This subroutine sets a range of digital-to-analog converters. The +/- 1.0 unit 
output range is scaled to +/- 10.0 volts so a single unit is worth .0049 volts. The 
arguments are: 

first dac (int) First converter to be set (numbering starts from 

zero.) 

last_dac (int) Last converter to be set (maximum is 7.) 

&dac_data_polnter (*) Starting location DAC data. Data are floating point 

values with a range of +/- 1.0. 


RETURNS: 

error ■ 0 Indicates valid transfer, 
error *-1 Indicates bad range. 

error > 0 Indicates "error” number of data words out of range. 


USES: 

Nothing. 


LIBRARY: 

None. 


SOURCE: 

/usr/csc/ele/setdac c.c 
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Analog I/O system subroutine. 


NAME: 

getadcO ( Fortran callable ) Sample a range of analog-to-digital 

converters. 


CALL: 

integer error, getadc 

error ■ getadc( f irst_adc, last_adc, adc_data_array) 


DESCRIPTION: 

This subroutine samples a range of analog-to-digital converters. The +/- 10.0 volt 
Input range Is scaled to +/- 1.0 units so a single bit is worth .00049 units. The 
arguments are: 

first adc (integer) First converter to be sampled (numbering starts from 
” zero.) 

last adc (integer) Last converter to be sampled (maximum is 63.) 

adc_data_array (real) Starting location for storing sample data. Data are 

floating point values with a range of +/- 1.0. 


RETURNS: 

error * 0 indicates valid transfer, 
error *-l Indicates bad range. 


USES: 

getadc_w.j . 


LIBRARY: 

/lib/acrl rt lib f.j 


SOURCE: 

/usr/csc/ele/getadc_f .c 
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Analog I/O system subroutine 


setdacO ( Fortran callable ) Set a range of digital-to-analog 

converters. 


CALL: 

integer error, setdac 

error * setdac( first dac, last dac, dac_data_arrray ) 


DESCRIPTION: 

This subroutine sets a range of digital-to-analog converters. The +/- 1.0 unit 
output range is scaled to +/- 10.0 volts so a single unit is worth .0049 volts. The 
arguments are: 

first dac (Integer) First converter to be set (numbering starts 
~ from zero.) 

last dac (integer) Last converter to be set (maximum is 7.) 

dac_data_array (real) Starting location DAC data. Data are floating point 

values with a range of +/- 1.0. 


RETURNS: 

error » 0 Indicates valid transfer, 
error »—1 Indicates bad range. 

error > 0 Indicates "error" number of data words out of range. 


USES: 

setdac w.j 


LIBRARY: 

/lib/acrl rt lib f.j 


SOURCE: 

/usr/csc/ele/setdac f.c 
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Analog I/O system subroutine 


NAME: 

■ot8ub() ( Fortran callable ) Send scaled voltages to the 8751 

micro-controllers which in turn set 
individual DACs. 


CALL: 

integer torque 

call raotsub ( torque ) 


DESCRIPTION: 

This subroutine sends the motor torque command data to the 8751 micro-controller 
boards which in turn load the data into the digital-to-analog converters* 

torque(6) (integer) Array of dimension 6 which contains the scaled data 

to be output on the 8751 DAC's. The range of the data 
is +2047 for +9*9951 volt output to -2048 for -10.0000 
volt output. 

REQUIRES: 

Call to motlnt. 


RETURNS: 

Nothing. 


USES: 

motsub w.j 


LIBRARY: 

/lib/acrl rt lib f.j 


criMPf* 1 • 

/Jpw/8751C0M/torsub.c 
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Analog I/O system subroutine 


NAME: 

fInIshO ( Fortran callable ) Close serial ports to 8751s. 

CALL: 

finishO 

DESCRIPTION: 

This subroutine closes the serial communication lines to the 8751s. There are no 
arguments. 

RETURNS: 

Nothing. 

USES: 

finish_w.j. 

LIBRARY: 

/lib/acrl_rt_lib_f. j 

SOURCE: 

/jpw/8751COM/torsub.c 
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Analog I/O system subroutines. 


NAME: 

thrust ( Fortran callable ) Set the discrete ports to activate the 

thrusters. 


CALL: 


integer*2 thrust - thrust ( x,y ) 


DESCRIPTION: 

This subroutine sets the states of the four discrete outputs on the Parallel 
Interface/Timer. It was designed for the thrusters on the SCOLE facility. The 
arguments can have one of three values: 1, 0, or -1 corresponding to positive, 
none, and negative thrust respectively. 

The arguments are: 

x (integer) State of x thruster, 

y (integer) State of y thruster. 


REQUIRES: 

Nothing. 


RETURNS: 

Nothing. 


USES: 

thrust w.j 


LIBRARY: 

/lib/acrl rt lib_f. 


SOURCE: 
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/j pw/TIMER/pitdsc.j 


Analog I/O system subroutine 


NAME: 

rtiae() (fortran callable) Mark the start of the real-time loop and 

change the sample time interval if required* 

CALL: 

call rtime(tau,flag,k) 

DESCRIPTION: 

An internal memory mapped timer is used to control the timing of real-time 
operations in the Charles River Computer. The programmable clock is required to 
generate a start pulse and a stop pulse for each sampling Interval of the control 
process. The maximum interval is eighty-five seconds and the minimum interval Is 5 
micro-seconds. The timer is a Motorola M68230 Parallel Interfac Timer (Pl/T) which 
provides versitile double buffered parallel Interfaces and 24-bit programmable timer 
for M68000 systems. 

Note: The call to rtime() should be made just Inside the real time loop. The 
arguments are: 

tau (real) is the sample period, 

flag (logical) is the Indicator to either maintain the same value of tau or 
pass in a new value, 

k (integer) is the timout parameter. If k Is returned from rtlme 

containing a 1, this indicates a normal return. If k is returned a 0, 
a timout has occurred, and appropriate action should be taken. The 
user must supply his/her own timout procedure. 

REQUIRES: 

Nothing. 

RETURNS: 

k 

LIBRARY: 

/lib/acrl_rtJLib_f.j 

Add " /jpw/rtirae.obj " to FORTRAN compile command. 

USES: 

trasetQ 
InitO 
rtwate() 
cktim() 
pint^w.j 

SOURCE: 

/jpw/TIMER/plnt.c 
/jpw/TIMER/plnt_w.m 
/jpw/rtime.for 
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TABLE II. TYPICAL THRDSTKR DATA FOR SIXTY PSI NOZZLE PRESSURE 


Peak Thrust 
Steady State Thrust 
Rise Time 
Thrust Duration 


0.641 lb 
0.32 lb 
0.032 seconds 
24 seconds 
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TABLE III. ACTUATOR PARAMETERS FOR CONTROL SYSTEM DESIGN 
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TABLE IV. SCOLE APPARATUS MASS PROPERTIES AND COMPONENT LOCATIONS 
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TABLE IV. SCOLE APPARATUS MASS PROPERTIES AMD COMPONENT LOCATIONS (CONT'D) 
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Figure 1. The SCOLE experiment apparatus. 
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Figure 2. Basic SCOLE structural assembly. 







Bored cylinder 
4.00 outer diameter 
t = 0.25 depth = 1.50 


Top surface t = 


Antenna 
Weight 4.76 lb 



O.750 outer diameter 


All units are in inches 
♦ Length does not include length of insert. 


Figure 5. Reflector assembly• 
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Figure 6. Universal joint suspension point. 
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Figure 7. Typical accelerometer calibration. 
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Figure 8a. Shuttle-mounted accelerometers 
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Figure 9a. Mast-mounted accelerometers. 
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Figure 8b. Coordinates and sensing axis of 
shuttle-mounted accelerometers. 
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Figure 9b. Coordinates and sensing axis 



of mast-mounted accelerometers. 
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Figure 10b. Coordinates and sensing axis of 

Figure 10a, Reflector-mounted accelerometers• 

reflector mounted accelerometers* 
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Figure 11. Typical rate sensor calibration 
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Figure 12a. Shuttle-mounted, three-axis 


rate sensor• 



Figure 13a. Mast-end-mounted, three-axLs 


rate sensor. 
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Figure 12b. Coordinates and sensing axis of 
8huttie-mounted rate senspr. 
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Figure 13b. Coordinates and sensing axis of 
mast—end—mounted rate sensor. 
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Figure 15b. Coordinates and actuation ITxls 


of shuttle-mounted CMG's. 
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Figure 16. Typical reaction wheel torque sensitivity curve. 



Figure 17a. Mast-mounted reaction wheels. 


Figure 17b. Coordinates of mast-mounted 


reaction wheels. 
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Figure 18b. Coordinate and axis of actuation 


for mast-end mounted reaction 
wheels. 


Figure 18a. Mast-end mounted reaction wheels. 
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Figure 19. Thruster air supply schematic. 
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Figure 20. Thruster startup and transients. 



Figure 21. Thrust magnitude and duration. 
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Figure 22a. Reflector-mounted thrusters. 
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Figure 22b* Coordinates of reflector-mounted thrusters* 
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OPTICAL ATTITUDE SENSOR 
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CMG CONTROLLER 


CMG COMMANDS 


Figure 23. SCOLE computer interfaces * 
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SUMMARIZE THE GRID SCHEME OF WILLIAMS AND 
MONTGOMERY 

DESCRIPTION OF THE LQG DESIGN FOR THE 
SCHEME 

PLANS FOR COMPLETION OF THE WORK 
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OBJECTIVE OF WORK 


DEVELOP k TEST A PRACTICAL SENSOR ARM 
SCHEME USING SCOLE 


APPROACH 

USE SCHEME PREVIOUSLY DEVELOPED FOR THE 
GRID DY WILLIAMS AND MONTGOMERY 
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GRID ARM SCHEME - SUMMARY 


USE SINGLE, ON-LINE, KALMAN FILTER 

USE SPRT TO CHECK THE ZERO MEAN CHARACTER 
OF THE ESTIMATED MEASUREMENT ERROR 

IF FAILURE IS DETECTED, ISOLATE USING 
FAILURE SIGNATURE IN THE ESTIMATED 
MEASUREMENT ERRORS 
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Logic for Failure Detection 
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Behavior of the 
SPRT Decision Variable 
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Failure threshold 






Effect of Bias Failure on Estimation 
Langley Grid Experiment 



613 














Failure Detection and Reconfiguration 

Langley Grid Experiment 
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TORQUE WHEELS USED FOR THE REGULATOR 




SCOLE CONFIGURATION 
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ABSTRACT 


The computer software available to those implementing identification 

and control algorithms on the Spacecraft Control Labratory Experiment 

(SCOLE) is described in the following paper. This paper is intended to 

be used as a guide and reference material. j 
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The Charles River Data System, commonly referred to as Chuck or the 
CRDS, is currently being used as the Aerospace Control Research Lab 
(ACRL) micro-computer for realtime programming for the SCOLE 
(Spacecraft Controls Laboratory Experiment). The CRDS is ba^ed on the 
Motorola 68020 microprocessor. The operating system is UNIX V with a 
Korn Shell. The computer has twelve serial ports and one parallel 
port. Terminals, five answer only communication lines, an originate 
only communication line, and a line printer are connected to the 
serial ports. 


Accessing the CRDS : 

The Charles River Data System (CRDS) can be accessed through LaTS 
(group name ACRL, group members 59382, 59383, 59384, 59385) or the 
Communication Machinery Corporation (CMC) Internet hardware and 
software. (Internet address 128.155.23.20) The following Internet 
applications are available for the UNIX-UNIX system communication: 


1. Remote copy (rep) 

-transfers file and directory copies 

2. Remote login (rlogin) 

-establishes remote UNIX system login 

3. Remote shell (rsh) 

-executes one UNIX command on a remote system 

The following Internet applications are available for the UNIX-UNIX 
and/or UNIX-non-UNIX system communications: 

1. File transfer (ftp) 

-transfers file and directory copies 

2. Network virtual terminal (telnet) 

-establishes remote system login 

Upon establishing communications with the CRDS, the user will be 
prompted for a login name and password. The CRDS system administrator 
will provide the user with a login name; the user Is encouraged to add 
a password. 

Editors : 

The Charles River Data System supports three editors-ved, vi, and ed. 
On-line help i»s iavailable for each editor by typing "help ved "y "man 
vi I more", Yman ed I more". Ved and vi are full screen editors. Ed is, 
a line editor. 

FIle Naming Conventions: 

File naming conventions exist on the CRDS. A table of appropriate file 
name extensions follows. These conventions must be followed to ensure 
proper functioning of the system procedures. 
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PROGRAM TYPE 


EXTENSION 


£ source programs 
assembly language 
object flies 
FORTRAN source programs 
Pascal source programs 
editor backup files 
listing files 

pre-processed assembly language 


. c 
. m 
. o 

. for 
. pas 

.bak (ved editor only) 
. 1st 

.s (FORTRAN 77) 


Control Characters : 

Control characters influence terminal processes. To use a control 
character, strike the CTRL key followed by the letter corresponding to 
the desired command. A table of control characters commonly used on the 
CRDS follows. 


CONTROL CHARACTER COMMAND EXECUTED 


~c 

A d 

A r 

~s 

~u 

n 


stops command execution 
returns to prompt level 
end of file 
user logged out 

continues display halted with ~s 
redisplays current input line 
suspends display such as listing 
contents of file to the screen 
erases current line 
suspends command execution 


Commands : 

On-line help is available for the following commands by typing "help" 
or "describe'* followed by the command name: 

addname, ar, archive, build, cat, cc, ccom68, cd, change, chmod, 
chown, command, common, compare, contig, copy, count, cpw, cut, date, 
debe, debug, delete, describe, difference, diskusage, dump,echo, 
email, false, find, gcc, hexodec, hist, id, kill. Id, link, list, 
lpd, lpmalnt, ‘l^q, lpr, makedir, match, message, mount, move, Hlce, 
p, paste, pmisk, #, prfmt, ps, pwd, repeat, resume, setdate, 
setmask, setprl, setuser, size, sleep, space, su, suname, 
suspend, sysversion, tall, ted, tee, term, time, translit, true, 
ttymodes, unmount, uname, unique, uptime, version, wall, wcc, when, 
whenq, whenrun, which, who, write, yes. 

The following summary presents the most commonly used commands on the 
CRDS: 
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Command 


Description 


cat "flle(s)" 

Displays contents of flle(s) or standard input 
if not supplied 

cd "dir" 

Change working directory to ’dir’ 

cp "f 1 lei" ' "f*l le2" 

Copy "fllel" to "flle2" 

cp "f1le(s)" "dir" 

Copy "file(s)" into "dir" 

date 

Displays the date and time 

echo "args" 

Displays "args" 

In "fllel" "file2" 

Link "filel" to "file2" 

In "flle(s)" "dir" 

Link "file(s)" into "dir" 

Is "f11e(s)" 

List "file(s)" 

Is "dir(s)" 

List files in "dir(s)" or in current directory if 
"dlr(s)" is not specified 

mkdir ”dir(s)" 

Create directory "dir(s)" 

mv "fllel" "f1le2" 

Hove "fllel" to "flle2" (simply rename it if both 
reference the same directory) 

mv "flle(s)" "dir" 

Move "file(s)" into directory "dir" 

ps 

List information about active processes 

pwd 

Display current working directory path 

rm "file(s)" 

Remove "flle(s)" 

rmdlr "dlr(s)" 

Remove empty directory "dir(s)" 

sort "flles(s)" 

Sort lines of "file(s)" or standard input if not 
supplied 

wc "file(s)" 

Count the number of lines, words, and characters 
"flle(s)" or standard input if not supplied 

who 

Display who’s logged in 

[Source: Charles 

River Data System’s Commands Manual] 


Subroutines : 

On-line help is available for the following commands by typing "help" 

or "describe" followed by the subroutine names: 

abs, acos, alarm, asctime, asln, astoc, astod, astof, astol, 
astol, atan, atan2, atod, atof, atoi, atol, atoq, atos, bsearch, 
capitalize, ceil, clock, close, conv, cos, cosh, creat, ctime, 
cvtlme, dayofyear, delay, drand, drand48, dup, dup2, cprlntf, 
execl, execle, execlp, execv, execve, execvp, exit, exp, fabs, 
fbufcnt, fclose, fexecl, fexecle, fexecv, fexecve, fflush, 
fgetc, fgetllne, fileopen, fllepos, flleread, flleseek, fllesize, 
filewrlte, find, findline, floor, flush, fmod, fopen, format, 
fprlntf, fputc, fputs, frac, fread, frexp, fscanf, fseek, fshow, 
fsstat, Ffcftat, ftell, ftime, ftoes, ftofs, ftok, fwrite, 'gamma, 
gcvt, gjetallargs, getargs, getbroken, getc, getchar, getcmask, 
getcwd, gethz, getllne, gettlme, gmtlme, gtty, hypot, 
index, loctl, ls_addressable, ltoq, JO, kill, ldexp, In, 
localtime, log, loglO, longjmp, lrand, lsearch, lseek, ltoq, 
map_lo_page, modf, open, parseranges, pause, peekc, pipe, 
poly, pow, proctlme, putc, putchar, qadd, qadd3, qcmp, 
qdiv, qdlv3, qdlvi, qmod, qmod3, qmul, qmul3, qneg, qneg2, 
qsort, qsub, qsub3, qtos, qtol, qtol, qtoq, qtos, quicksort, 
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rand, random, randomize, read, rewind, scanf, setjmp, setwd, 
signal, sin, sinh, sleep, spfun, sprlntf, srand, sscanf, 
ssignal, stat, status, stoq, strcat, strcatl, strchr, strcmp, 
strcpy, streql, strlndex, strlen, strncat, strncmp, strncpy, 
strrchr, system, tan, tanh, tell, time, times, tlmestr, 


tlmetoSitr, ttyname, wait, yO. 

.i 11*. 


Compilers: 


FORTRAN 77, C, and Pascal compilers are available on the CRDS. On-line 
help is available for each compiler by typing "help" followed by 
f77 (FORTRAN 77), cc (C), and pascal (Pascal). 


Many macros exist to compile, link and assemble code. The Following Is a 
brief description of each macro. On-line help Is available for cm, af, 
and as by typing "help" followed by the macro name. 


cm : 

file: FORTRAN 77 main program source code "main, for" 

command: cm main (prefix only) 
use: produce object code 


as: 

file: FORTRAN 77 subroutine source code "sub.for" 

command: as sub (prefix only) 
use: produce object code 

af: 

file: FORTRAN 77 main program source code "maln2.for" 

command: af main2 (prefix only) 

use: produce executable code by linking object code to all 

available system libraries (SCOLE, math, realtime, strip 
chart Included) 

executable code Is stored In the file "maln2" 

Note : This macro is not to be used with main programs 
calling subroutines not found In the system 
1lbrarles. 


C: 

file: C main program or subroutine source code "maln.c" or 

"sub.c" 

command: C main or C subc (prefix only) 

use: produce C object code main.o or subc.o 

S: 

file: FORTRAN pre-processed assembly code 

"file.s“ produced with f77 command 
command: S file (prefix only) 

use: produce assembly language code flle.m 

A: 

file: FORTRAN assembly code "flle.m" produced with S macro 

command: A file (prefix only) 

use: produce object code flle.o ready to be linked. 
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Note : macros cm and as perform the same functions as the f77 command 
used in conjunction with the S and A macros. 


Linker : 

The CRDS linker command is "id". On-line help is available for the linker 
by typing "help Id". The standard link line is as follows: 

Id FORTRAN maln.o C/FORTRAN subrout 1 ne(s).o -Klibrary names) o=fn 


Id C maln.o C/FORTRAN subroutine(s) .o -Klibrary names) o=fn 
(fn = executable file name) 

One or more of the following libraries may be linked to the object code: 
LIBRARY COMMAND 


scole : 

SCOLE FORTRAN library 

(-lscole) 

real : 

C realtime library 

(-lreal) 

mt8500 : 

C MT8500 strip chart 



recorder library 

(-lmt8500) 

matx : 

FORTRAN matrix library 

(-lmatx) 

unos : 

UNOS system library 

Note: This library must 
always be linked and must 
be linked as the last 

1lbrary. 

(-lunos) 


On-Line Help : 

Every effort has been made to make the CRDS a user-friendly machine. 
Help and system updates are updated as required. The 

following menu presents all of the on-line help and update reports of 
the CRDS. 

Should any questions arise pertaining to software and/or hardware 
problems with the CRDS, please feel free to contact Danette B. Lenox at 
X46619 or x44473 ((804)864-6619 or (804)864-4473). 


On-line help is available on the following topics: 

1. SCOLE 'Help: 

A. ' type "h_make" for help with the SCOLE program build 

macro "make" 

B. type "flow" for a copy of the SCOLE program flowchart 

2. System bug report: 

A. UNOS VII bugs type "bJJNVII" 

B. UNIX bugs type "bJJNIX" 

C. C compiler bugs type "b_UNCII" 

D. FORTRAN compiler bugs type "b_F0RTRAN" 
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3. System updates: 

A. type "updates" 

4. Basic system subroutine and command help: 

A. type "llst_all" to see a listing of the names of the 
basic commands for which help Is available 

5. Complied help: 

A. type "absoft" for help with the absoft F0RTRAN77 compiler 

B. type "comp_pascal" for help with the PASCAL compiler 
G. Library help: 

A. type "real" for help with the real time library 

B. type "math" for help with the matrix math library 

C. type "LINPACK" for help with the LINPACK libraries 

D. type "mt8500" for help with the MT8500 strip chart 
recorder library 
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SCOLE SOFTWARE OVERVIEW 
SECTION II 
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The SCOLE model Is operated by a set of software generlcally referred 
to as the SCOLE program. Control laws can be easily Implemented and 
tested on the model through minor SCOLE software alterations. 

1. SCOLE Files 

Each user must create a directory called "SCOLE" on his account. 

Subroutine and main program source and object code will be located on 
directory SCOLE. Flies main, for, cpmenu.for, preexp.for, 
posexp.for, rtloop.for, makefile, gains.dat and select.set 
on directory /I1b/SCOLE/TEMPLATE should be copied onto the user’s 
SCOLE directory. 

These files perform the following functions: 

A. main.for : 

.. source code main program 

B. cpmenu.for: 

.. source code for controls parameters menu 

C. posexp.for : 

.. source code for post-experiment menu 

D. preexp.for : 

.. source code for pre-experiment menu 

E. rtloop.for : 

.. source code for real time loop 

F. makefile : 

..makefile used to create SCOLE program (refer to topic #2) 

G. gains.dat : 

..data file used by SCOLE library routines altered by 
variable setting In menus 

H. select.set : 

..data file used by SCOLE library routines altered by 
variable setting In menus 

Each of these files should be altered according to the user’s require¬ 
ments. The codes have been "commented" to provide guidelines for 
alterations. All additional source code for subroutines written by 
the user for Implementation of the SCOLE contol law should be 
located on SCOLE and must have the appropriate extension. 

2. Building SCOLE Program 

The command "make" builds the SCOLE program with the code located on 
the user’s SCOLE directory using the directions found In makefile, the 
math library (^ib/llbmatx. a), the strip chart library (/llb/ll|bmt8500. a), 
the real time; library (/llb/1lbreal. a), the UNOS library (/llb/libunos. ^), 
and the SCOLE library (/lib/1ibscole.a). 

The command Is invoked by typing: 

make [exec=name]. 

The SCOLE program executables are by default written to file ’scole’. Should 
the user desire to save the executables In another file, simply append 
"exec = name" where name Is an arbitrary file name to which the 
executables are written. 
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The formation of the program Is based on relationships expressed In 
’makefile’ as follows: 

The modification dates of source code on the user’s SCOLE directory 
are compared to that of the corresponding object code on SCOLE. If 
the date of the source code Is more recent, the subroutine Is 
recompiled. The modification dates of the two SCOLE Include flies 
(/HOLD/SCOLE/varlbs.comm and /HOLD/SCOLE/ blk.comm) and any user 
supplied Include files are also checked. If these dates are more 
recent than the object codes Including the files, the corresponding 
source Is recompiled. If errors occur In compilation, formation of 
the SCOLE program will not be successful. FORTRAN syntax errors can 
be noted by studying the ".1st" file related to the subroutine In 
question. If no errors occur during compilation, the object code Is 
linked to the libraries listed above, forming a set of executable code 
stored in the file scole or "name". Watch out for "unresolved 
externals". When user adds files to the SCOLE directory, the source 
name must be added to the FSUBS list In makefile and the object name 
must be added to the FOBJ list in makefile. When the user adds 
Include files to the SCOLE directory the file name must be added to 
the INCVC list. 

3. Running the SCOLE Program 

The SCOLE program Is run by simply typing scole or "name". 

If the user encounters timeouts during the real time execution, Increase 
the sample period and please contact the system administrator for informat 1 
on running the program with priority privileges. 

4. Scole Software Variables 

The following variables are commonly used In the SCOLE software 
environment and are passed by way of common blocks In Include file 
/HOLD/SCOLE/var1bs.comm: 


VARIABLE 

DESCRIPTION 



shurgr 

shuttle rate gyro roll 



shurgp 

shuttle rate gyro pitch 


shurgy 

shuttle rate gyro yaw 



refrgr 

reflector rate gyro roll 


refrgp 

reflector rate gyro pitch 


refrgy 

reflector rate gyro yaw 


shaccx 

shuttle acceleration x 

direction 


shaccy 

shuttle acceleration y 

direction 


shaccz 

shuttle acceleration z 

direction 


mlaccx 

mast station 1 -upper- 

acceleration 

X 

mlaccy 

mast station 1 -upper- 

acceleration 

y 

m2accx 

mast station 2 -lower- 

acceleration 

X 

m2accx 

mast station 2 -lower- 

acceleration 

y 

refacx 

reflector acceleration 

x direction 


refacy 

reflector acceleration 

y direction 


shtorr 

shuttle torque roll 



shtorp 

shuttle torque pitch 




direction 

direction 

direction 

direction 






shtory 

refwhx 

refwhy 

refwhz 

mlwhx 

ml why 

m2whx 

m2why 

tx 

ty 

phi 

theta 

psl 

P 

q 

r 

fx 

fy 

ctheta 

stheta 

cpsl 

spsl 

cphl 

sphl 

bis 


shuttle torque yaw 

reflector wheel x direction 

reflector wheel y direction 

reflector wheel z direction 

mast station 1 -upper- wheel x direction 

mast station 1 -upper- wheel y direction 

mast station 2 -lower- wheel x direction 

mast station 2 -lower- wheel y direction 

thrust x direction 

thrust y direction 

bank angle 

pitch angle 

yaw angle 

shuttle rotational velocity about x 

shuttle rotational velocity about y 

shuttle rotational velocity about z 

thruster gains x direction 

thruster gains y direction 

cosine theta 

sine theta 

cosine psl 

sine psl 

cosine phi 

sine phi 

thruster least square coefficient 


5. SCOLE FORTRAN Library Routines 

The following subroutines are found in the SCOLE library. Users are not 
given privileges to alter the subroutines. 


MT8500 STRIP CHART RECORDER SUBROUTINES: 


chan 

cscale 

lbmenu 

scree 

spdset 

stmenu 

strcht 

strset 


change number of strip chart channels 

accept new strip chart scale factors 

set strip chart recorder labels 

replays saved data to strip chart 

set strip chart recorder speed 

strip chart recorder menu 

send data to the strip chart recorder 

select program variable to be output to the strip 

chart recorder 


PARAMETER DEFINITION AND FILE I/O SUBROUTINES: 


collnt 

dltlme 

ertrap 

lselek 

lhelp 

readf 

skipln 

tscole 


gain and dead band levels for thruster slew menu 

duration and Interval time change 

error trapping In keyboard input 

user menu selection option 

on-line help 

reads standard SCOLE data files 

skips line 

get current date 
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CMC SUBROUTINES: (CMC software In process of being updated) 


cmgtor send 3 axis shuttle torque commands to cmg controller 

cmgtst cmg test 

gettor get shuttle torque commands from Joystick 

glmbal *cmg gimbals - four torque commands to gimbals' 

sensor prints scaled sensor measurements from cmgs 

MEASUREMENT/DATA ACQUISITION SUBROUTINES: 


attlt 

damlnt 

getac 

getyz 

lsqang 

sblas 

slwlnt 

toggle 


determine roll and pitch attitude of shuttle from 

accelerometers and rate sensor 

thruster damping menu 

get accelerometer data 

get y and z data from the op_eye 

least square shuttle attitude estimation 

determines sensor bias by averaging over twenty seconds 

accepts gain and dead band levels 

toggle electromagnet (not In use) 


OPTICAL SENSOR SUBROUTINES: 

(Optical Sensor software In process of being 
updated) 

opacal calibrate optical sensor 

opang returns yaw and pitch angles as determined from optical 

sensor outputs 

opashw prints yaw and pitch angles as determined from optical 


opatst 

opmenu 

qqstop 

yz 


sensor outputs 
test optical sensor 

sets up optical sensor calibration and test 
stops program 

get y and z data from the op_eye 


DATA OUTPUT SUBROUTINES: 


stof11 
stoset 
uwrlte 


store experimental data In arrays to be output to file 
select program variable to be saved on file 
output data to file 


6. SCOLE Realtime Library Calls 

The realtime library consists primarily of the following FORTRAN 
callable C routines. Help Is available for each routine by typing 
"describe" followed by the routine name In question. 


getadc - 

This 'subroutine samples a range of analog-to-dlgital converters. 
+/- 10.0 volt Input range Is scaled to +/- 1.0 units so a single 
unit Is worth .0049 volts. 


lswtch - 

This subroutine checks the state of a switch as low or high. 

oswtch - 

This subroutine sets the state of a switch as low or high. 


The 
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rtime - 


An internal memory mapped timer Is used to control the timing of 
real-time operations In the Charles River Computer. The programmable 
clock Is required to generate a start pulse and a stop pulse for each 
sampllhg Interval of the control process. The maximum lhterval.is 
eighty-five seconds and the minimum Interval Is 5 micro-seconds. The 
timer 'lsl-a Motorola M68230 Parallel Interface Timer '(PI/T) which 
provides versatile double buffered parallel Interfaces and 24-bit 
programmable timer. 


setdac - 

This subroutine samples a range of dlgltal-to-analog converters. The 
+/- 1.0 unit output range Is scaled to +/- 10.0 volts so a single 
unit Is worth .0049 volts. 


strip - 

This function outputs data to the analog strip chart recorder. 


thrust - 

This function sets the states of the four discrete outputs on the 
Parallel Interface/Timer. It was designed for the thrusters on the 
SC0LE facility. The arguments can have one of three values: 1, 0, or 
-1 corresponding to positive, none, and negative thrust respectively. 

rnum - 

This function returns the current date since Jan 1, 1980 In seconds. 

mprep - 

This function clears the terminal screen (Wyse50 terminal). 

mtfeed- 

Thls function advances the paper in the MT8500 chart recorder 
to the next top-of-form mark (form feed). 

mtrls- 

Thls function releases the MT8500 chart recorder from the real¬ 
time mode and leaves It in a ready state. 

mtstar- 

This function starts the MT8500 chart recorder’s waveform 
printing - printing of scale grid can be turned off and on via 
subroutine grid argument. 

mtstop- 

Thls function stops the waveform printing from the MT8500 chart 
recorder. 

mtspee- 

Thls function sets the waveform printing speed of the MT8500 
recorder. 

mtchan- 

Thls function sets the number of channels to be printed from the 
MT8500 recorder. 

cmgtor- (software in process of being updated) 

This function commands the torque wheels. 
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A table of 
follows. 

the C subroutines 

commonly called from the 

SCOLE program 

SUBROUTINE/ LIBRARY 

FUNCTION 

C DECLARATIONS 

FORTRAN DECLARATIONS 

thrust 

real 

lnt thrust(x.y) 
lnt *x, *y 

integer thrust 

Integer x, y 

Iswtch 

real 

short lswtch(n) 
lnt *n 

loglcal'2 Iswtch 
Integer n 

oswtch 

real 

short oswtch(state,n) 
lnt 'state, 'n 

integer*2 oswtch 
integer n, state 

rtlme 

real 

lnt rtlme(flag,tau) 
float *tau 
lnt 'flag 

loglcalM rtlme 
real tau 
logical flag 

setdac 

real 

long setdac( 
dac_dat a_pointer, 
last_dac,first_dac) 

integer*4 setdac 
real dac_data_pointer, 
last_dac,first_dac 

getadc 

real 

long getadc( 
adc_data_polnter, 
last_adc,flrst_adc) 

lnteger'4 getadc 
real adc_data_poInter, 
last_adc,flrst_adc 

mtfeed 

mt8500 

long mtfeedU 

integer*4 mtfeed 

mtrls 

rat8500 

long mtrls() 

lnteger'4 mtrls 

mtstar 

mt8500 

long Btstar(grld) 
lnt 'grid 

lnteger'4 mtstar 
Integer grid 

mtstop 

rat8500 

long ntstop() 

lnteger'4 mtstop 

mtspee 

nt8500 

long atspeeC 
unlts_arg,speed_arg) 
char •unlts_arg 
float 'speedjarg 

integer*4 mtspee 
character*10 units_arg 
real speed_arg 

mtchan 

mt8500 

long mtchan(nchan) 
long *nchan 

lnteger*4 mtchan 
Integer nchan 

cmgtor 

real torz, 
(software 

cmg long cmgtor(torz,tory 

torz) 

tory, torx 

In process of being updated) 

, integer*4 cmgtor 

strip 

real 

long strlp(data) 
lnt 'data 

lnteger'4 strip 
real data 

mprep 

real 

long mprep() 

integer*4 mprep 

rnum 

real 

long rnura(ans) 
long *ans 

integer*4 rnum 
lnteger*4 ans 
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7. SOFTWARE FLOWCHART 




1 - 

Iselek 



| - 

readf 



I - 

sblas 



1 ~ 

cmgraen 



1 - 

stmenu-1- lbmenu 1 


1 - (1) preexp 

-1 “ 

stoset 


i 

| - 

opmenu 


i 

1 - 

dltlme 1- (2) slwlnt 


i 

1- (1) 

cpmenu-1- (2) damlnt 


1 

1 " 

tstjet 1- (1) colint 


i 

1 " 

ertrap 

(1) 

main 1 




1 

I - 

strcht 


1 

l-(3) 

ntspeed 


1 - (1) rtloop 

- l-(3) 

mtstart 


1 

l-(2) 

rt lme 


1 

l-(2) 

getadc 


1 

1 - 

attlt 


1 

l-(2) 

lswtch 


! 

1 " 

damp-1- (2) lswtch 


1 

1 * 

lsqang- - —1- getac 


1 

I - 

slew 1- getyz 


1 

l-(2) 

thrust 


1 

1 -(2) 

setdac 


1 

1 - 

stofll 


1 

1 -(2) 

oswtch 


1 

1 - 

cagtor 


1 

l-(3) 

atstop 


1 

i 

l-(3) 

atfeed 


1 

1 

( — 

iselek 


1 - (1) posexp 

- 1 - 

uwrlte 


1 

1 

1 — 

scree 


1 

1 - (3) mtrls 



(1) 

user file 



(2) 

C realtime library routine (/llb/llbreal. a) 

(3) 

C strip chart library routine 

(/lib/llbmt8500. a) 


all other routines 

from SCOLE library (/llb/llbscole.a) 
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system 


I- ertrap 











SCOLE TEMPLATE FILES SOURCE CODE(/Iib/SCOLE/TEMPLATE) 

SECTION III 


I 
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program main 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C DOCUMENTATION c 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c This main program calls the 3 main modules - preexperiment, 

c real time loop and postexperiment. This is the main calling module, 

c 
c 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c DATA/VARIABLE DICTIONARY c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

c 


include /HOLD/SCOLE/varibs. comm 


+ 

+ 


Integer iabrt 
common /store/ bias(25) 
common /scale/ sfsrr, sfsrp, 
sfsax, sfsay, 
sfrax, sfray, 

common /dimen/ xd, yd, zd 
common /cdebug/dbug 
dimension stodatC10,10000) 


sfsry, 

sfsaz, 

sfstr, 


sfrrr, 

sfmlx, 

sfstp, 


sfrrp, 

sfmly, 

sfsty 


sfrry, 

sfm2x, 


sfm2y, 


logical dbug 
character ans*3 


c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c PROCEDURE SECTION c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
dbug=.false, 
iabrt * 0 


USERS CAN ADD CODE TO READ FILES, PERFORM INITIALIZATION OPERATIONS, 
OR WHATEVER ONE DESIRES IN THIS CODE 


10 


c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


if(dbug)print*, ’DEBUG MODE ON' 

THE MAIN PROGRAM FIRST CALLS SUBROUTINE PREEXP WHICH PRESENTS 
THE PRE-EXPERIMENT MENU BY WHICH USERS CAN INTERACTIVELY SET 
UP PARAMETERS 

call preexp (iabrt) 

IABRT IS A VARIABLE SET TO 1 IN "PREEXP.FOR" WHEN THE USER CHOOSES 
TO EXIT THE SCOLE PROGRAM 


if(dbug)print *,'call preexp from main* 
if (iabrt.eq.l) goto 10000 
if(dbug)print *,'call rtloop from main' 


NEXT THE MAIN PROGRAM CALLS SUBROUTINE RTLOOP WHICH CONTAINS THE 
REAL TIME LOOP OR THE ACTUAL EXPERIMENT 


E 

I 
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c 


call rtloop(stodat) 


i 


g 


c STODAT IS THE ARRAY IN WHICH DATA ARE SAVED IN THE REAL TIME LOOP 

c VARIABLES TO SAVE ARE CHOSEN BY SELECTING THE MAIN MENU OPTION TO 

c SET-UP THE DATA FILE 

c 

if(dbug)print *, p call posexp from main’ 
c 

c AFTER THE EXPERIMENT IS COMPLETE, THE MAIN PROGRAM CALLS POSEXP - 

c THE POST-EXPERIMENT MENU FOR DATA PROCESSING 

c 

call posexp (iabrt,stodat) 
if (iabrt.eq. 1) goto 10000 
goto 10 
10000 continue 


c RELEASE MT8500 STRIP CHART RECORDER 
c 

im = mtrls() 

stop 

end 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*******«•******************•*********«•*•********•»*********************+** 

* 

Block data subprogram * 

* 

t**«t**tt*»f******»**«**«***§*«**§«•§§***#**§*»****************♦*********** 


USERS SHOULD NOT DELETE ANY OF THE INCLUDE FILES BUT ARE FREE 
TO ADD ANY OF THEIR OWN FILES OR DATA 


blockdata const 


inc1ude /HOLD/SCOLE/var1bs. comm 
inc1ude /HOLD/SCOLE/b1k.c omm 


c 


end 
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subroutine cpraenu 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c This is the menu for ’SET UP CONTROL PARAMETERS’ module, 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c DATA/VARIABLE DICTIONARY c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
integer ians, mode, lhno, idef, iopt 

character nans*l, procse*20, hlpfil*20 

common /info/ procse, hlpfil 
common/cdebug/dbug 
logical dbug 
integer*4 iprep, mprep 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c PROCEDURE SECTION c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
if(dbug)print*,’Subroutine cpmenu: * 
c 

c USERS CAN ADD AND/OR DELETE CONTROL PARAMETER OPTIONS IN THIS MENU 
c 

c THIS MENU IS PRESENTED AS AN ILLUSTRATION 
c 

10 iprep = mprepO 

write!*, * )’-’ 

write!*, * )’l ** CONTROL PARAMETER MENU ** I* 

write!*, * )’-’ 

write!*, * )’l 1) SET THRUSTER DAMPING PARAMETER. I’ 

write!*, * )’l 2) SET SLEW PARAMETER. I’ 

write!*, * )*l 3) SET COLLOCATED MAST DAMPER PARAMETERS!’ 

write!*, * )’l 4) RETURN TO PREVIOUS MENU. [DEFAULT] I’ 

write!*, * )’-’ 

imaxp = 4 
lhno = 10 
idef = 4 
imode = 7 

call iselek (imaxp,imode,ians,nans,lhno,idef) 
if (nans. eq. ’M’) then 
goto 10 
else 

iopt = ians 
end if 
c 

c USERS SHOULD INCLUDE/DELETE OPTION CODE IN THIS IF. . . THEN. ..ELSE BLOCK 
c TO CORRESPOND TO HIS/HER OPTIONS 
if (iopt.eq.1] then 
call damint 

else if (iopt.eq.2) then 
call slwint 

else if (iopt.eq.3) then 
call colint 

else if (iopt.eq.4) then 
goto 10000 
end if 
goto 10 
10000 return 
end 
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subroutine posexp (iabrt,stodat) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c DOCUMENTATION c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

c This is the menu for the ’POST-EXPERIMENT MENU’ module 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c DATA/VARIABLE DICTIONARY c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
save dbug 

dimension stodat(10,10000) 

integer iabrt, ians, mode, lhno, maxp, idef, runs 

character nans*l, procse*20, hlpfil*20 

common /time/ period, ttime, runs 

common /info/procse,hlpfi1 

common/cdebug/dbug 

logical dbug 

integer*4 iprep, mprep 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c PROCEDURE SECTION c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
if(dbug)print*,’Subroutine posexp: ’ 
c 

c OUTPUT FILE FORMATS INCLUDE MATRIX-X, MATLAB, FORMATTED, UNFORMATTED 

c FILES OPTIONS PRESENTED IN LIBRARY ROUTINE UWRITE 

c 

10 iprep = mprep() 

c 

c USER IS GIVEN THE CHOICE AT THE END OF THE EXPERIMENT TO 

c SAVE DATA IN A FILE, REPLAY IT TO THE STRIP CHART, OR 

c END THE EXPERIMENT 

c 

write!*,*) 
write!*,*) ’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
write!*,*)’ 
c 

lhno = 10 
mode = 6 
idef = 0 
maxp = 4 

call iselek (maxp,mode,ians,nans, lhno, idef) 
if (nans.eq.’M’) then 
goto 10 

else if (ians.eq.1) then 
goto 10000 

else if (ians.eq.2) then 

call uwrite(runs,period, stodat) 
else if (ians.eq.3) then 
call scree 


I ** POST-EXPERIMENT MENU ** I 


1) PREPARE FOR NEXT RUN 

2) SAVE DATA TO OUTPUT FILE 

3) REPLAY SAVED DATA TO STRIP 

CHART 

4) STOP EXPERIMENT 
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else 

write(*,124) 

124 formatC Are you sure ? (y/n) - 

maxp = 0 
mode = -13 

call iselek (maxp,mode,ians,nans,lhno,idef) 
If (nans.eq.’Y*) then 
iabrt = 1 
goto 10000 
end If 
end if 
goto 10 
10000 continue 
c 

c USERS CAN ADD POST-PROCESSINF CODE AT THIS POINT 
c 

end 
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subroutine preexp (iabrt) 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c This is the menu for the ’PRE-EXPERIMENT MODULE' cc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccc 
c DOCUMENTATION c 

ccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccc 
c PROCEDURE SECTION c 

ccccccccccccccccccccccccccccccccccccccccc 

c 

include /HOLD/SCOLE/varibs. comm 
common /info/ procse, hlpfil 
common /cdebug/dbug 
logical dbug 

integer c, rret, iabrt, runs, ians, idef, mode, lhno, maxp 
integer system 
Integer*4 mprep, iprep 

character schar*l, nans*!, procse*20, hlpfil*20 
if(dbug)print*,'Subroutine preexp:' 
c 

c Set default values 
c 

iabrt = 0 
c 

c SCOLE LIBRARY ROUTINE TO READ FILES GAINS. DAT AND SELECT. SET ON SCOLE 

c 

call readf 

if (dbug)print*, ’ Returned from readf:' 

1000 continue 

900 iprep=mprep() 

c 

c TTIME AND PERIOD ARE SUPPLIED BY THE USER AS REQUESTED IN THE MAIN 

c MENU 

c 

runs=ttime/period 
c 

c USERS SHOULD READ ALL OF THEIR DATA FILES AT THIS POINT 

c ALL VARIABLES READ SHOULD BE PASSED AMONG THE SUBROUTINES 

c ON HIS ACCOUNT BY MEANS OF INCLUDE FILES 

c 

c USERS ARE FREE TO ADD/DELETE MENU OPTIONS AS LONG AS THEY TAKE 

c INTO ACCOUNT HOW THE OPTIONS ARE USED AND WHAT SUBROUTINES 

c THEY CALL 

c 

write(*, * )'- ’ 

write(* )’l ** PRE-EXPERIMENT MAIN MENU ** I’ 

wrlte(*, • )’- ’ 

write(\ * )’l 1) OBTAIN RATE SENSOR AND ACCELEROMETER I' 

write(*, * ) 9 I BIAS. I 1 

write(*, * )’l 2) SET UP AND TEST CMGS. I' 

write(*, * )'l 3) SET UP STRIP CHART. I' 

write(*, * )'l 4) SET UP DATA FILE. I ' 

wrlte(*, * )'l 5) SET UP OPTICAL SENSOR L’ 
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wrlte(“, 

* 

)’ i 

6) 

SET UP DURATION AND INTERVAL TIME. 

1 ’ 

write(* f 

* 

)■ i 

7) 

SET UP CONTROL PARAMETERS. 

1 ’ 

write(*, 

* 

)■ i 

8) 

TEST JET THRUSTERS. (NOT OPERATIONAL) 

1 ’ 

write(*, 

* 

n 

9) 

START EXPERIMENT 

1 ’ 

write(* f 

• 

)• i 

10) 

TOGGLE DEBUG MODE 

1 ’ 

write(*, 

* 

)’ i 

11) 

ENTER OPERATING SYSTEM 

1 ’ 

write(*, 

write(*, 

* 

* 

}’ i 
)*- 

12) 

QUIT. 

I ’ 


lhno = 10 
mode = 6 
idef = 0 
maxp * 12 

call iselek (maxp,mode,Ians,nans,lhno,idef) 
iprep = mprepO 
if (nans.eq.’M*) then 
goto 900 
else 

iopt » Ians 
end if 
c 

if (iopt.eq. 1) then 
call sbias(period) 
else if (iopt.eq.2) then 

print*, > **** # **THIS SUBROUTINE NOT IMPLEMENTED***** 
c call cmgmen 

else if (iopt.eq.3) then 
call stmenu 

else if (iopt.eq.4) then 
call stoset 

else if (iopt.eq.5) then 
call opmenu 

else if (iopt.eq.6) then 
call ditime(period,ttime) 
else if (iopt.eq.7) then 
call cpmenu 

else if (iopt.eq.8) then 
c call tstjet 

c this subroutine is not in use at this time 1-20-87 

print*,'*******THIS SUBROUTINE NOT IN USE AT THIS TIME******* 1 
go to 900 

else if (iopt.eq.9) then 
goto "10000 

else if (iopt.eq.10) then 
dbug * .not. dbug 
else if (iopt.eq.il) then 

write(*,*)* Use A D to return to this program' 
is » system("command") 
else if (iopt.eq.12) then 
writeT*,124) 

124 formate Are you sure ? (y/n) - ’ ,\) 

maxp = 0 
mode = -13 

call Iselek (maxp,mode, ians,nans, lhno, idef) 
if (nans.eq.'Y') then 
iabrt = 1 
goto 10000 
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else 

goto 900 
end if 
endif 
goto 900 
c 

100 format (a, f7. 2, a, f6. 4, a) 
c 

9996 call ertrap(ios,1,0) 
goto 1000 

9997 call ertrapdos, 1,0) 
goto 900 

10000 return 
end 
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This is the main subroutine for the real-time loop - all 
estimation, control, and data storage functions should be 
accomplished in this subroutine. 


c 

c 

c 

c 

c 

c 

c 

c 


* 

* 

* 

* 

* 


subroutine rtloop(stodat) 
c 

c USE COMMON BLOCKS TO PASS DATA BETWEN REAL TIME ROUTINES TO 

c TO IMPROVE PROCESSING SPEED, 

c 

include /HOLD/SCOLE/varibs. comm 
save 

dimension dac0(8) 
dimension stodatC10,10000) 

dimension adc(18), dac(8), sfsens(18), adcs(18) 

common/cdebug/dbug 

common/store/bias(25) 

common /scale/ sfsrr, sfsrp, sfsry, sfrrr, sfrrp, sfrry, 

+ sfsax, sfsay, sfsaz, sfmlx, sfmly, sfm2x, sfm2y, 

+ sfrax, sfray, sfstr, sfstp, sfsty 

common /dimen/ xd,yd,zd 
c 

c The Absoft compiler accepts only the first array argument and first 

c element of the common block to which the first array argument is 

c equivalenced. All subsequent equivalences are performed in the order 
c of array elements and common block variables 



equivalence 

( 

adcs(1), 

shurgr) 

_ 

c 

equivalence 

( 

adcs(1), 

shurgr) 


c 

+ 

.( 

adcs(2), 

shurgp) 


c 

+ 


adcs(3), 

shurgy) 


c 

+ 

.( 

adcs(4), 

refrgr) 

e: 

c 

+ 

.( 

adcs(5), 

refrgp) 

= 

c 

+ 

.( 

adcs(6), 

refrgy) 

i 

c 

+ 

.( 

adcs(7), 

shaccx) 

i 

c 

+ 

,c 

adcs(8), 

shaccy) 

= 

c 

+ 

,( 

adcs(9), 

shaccz) 


c 

equivalence 

( 

adcs(10), 

mlaccx) 


c 

+ 

,( 

adcs(11), 

ralaccy) 


c 

+ 

,( 

adcs(12), 

m2accx) 


c 

+ 

.( 

adcs(13), 

m2accy) 


c 


,( 

adcs(14), 

refacx) 

fi 

c 

+ 

,( 

adcs(15), 

refacy) 


c 

+ 

A 

adcs(16), 

shtorr) 

= 

c 

+ 

A 

adcs(17), 

shtorp) 

■ 

c 

+ 

A 

adcs(18), 

shtory) 

_ 


c 

equivalence (dac(l), refwhx) 
c 



equivalence (sfsens(1), 

sfsrr) 

_ 

c 

equivalence (sfsens(1), 

sfsrr) 

_ 

c 

+ ,(sfsens(2). 

sfsrp) 

: . 

c 

+ ,(sfsens(3), 

sfsry) 

= 

c 

+ ,(sfsens(4), 

sfrrr) 
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c 

+ 

,(sfsens(5), 

sfrrp) 

c 

+ 

,(sfsens(6), 

sfrry) 

c 

+ 

,(sfsens(7), 

sfsax) 

c 

+ 

,(sfsens(8), 

sfsay) 

c 

+ 

,(sfsens(9), 

sfsaz) 

c 

+ 

,(sfsens(10), 

sfmlx) 

c 

+ 

,(sfsens(11), 

sfmly) 

c 

+ 

,Csfsens(12), 

sfm2x) 

c 

+ 

,(sfsens(13), 

sfm2y) 

c 

+ 

,(sfsensC14), 

sfrax) 

c 

+ 

,(sfsens(15), 

sfray) 

c 

+ 

,(sfsens(16), 

sfstr) 

c 

+ 

,(sfsens(17), 

sfstp) 

c 

+ 

,(sfsens(18), 

sfsty) 


c 

c Function declarations 
c 

integer*2 oswtch, io 
logical*2 iswtch, il, ils 
logical*4 rtime 
integer*4 cmgtor 
integer*4 iprep, mprep 
integer*4 setdac , getadc , thrust 
integer*4 mtspeed, mtstart, mtstop, mtrls, mtfeed 
c 

character*! schar 
c 

logical flag, reset, cmgclr,dbug 
logical*4 insw(8) 

Integer c, rret 
integer runs 

integer reply, x, y, dampx, dampy, slewx, slewy 
c 

data insw / 8*.false. / 

data dacO / 8*0.0 / 

if(dbug)print*,* Subroutine rtloop: ’ 

c 

c 

c DELETE THIS CALL TO STRCHT IF THE USER DOES NOT WISH TO 

c USE THE MT8500 STRIP CHART RECORDER 

c 

call strcht 
iprep = mprep() 
c 

write(*,*) 
write(*,124) 

124 format(’ EXPERIMENT RUNNING. 

1/,’ WAIT FOR MENU 1 ) 

c 

reset = .true, 
flag * .true, 
c 

refrgr=0.00 
refrgp^O.00 
refrgy=0.00 
c 
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shtorr =0.0 
shtorp =0.0 
shtory =0.0 
c 

c START STRIP CHART RECORDER - KT850Q 

c OPTIONAL CODE - STRIP CHART NOT NECESSARY 

c DELETE LINES BETWEEN ASTERISKS IF STRIP CHART NOT DESIRED 

£*»»*****««****•**••**** #•**•«**•«***«************************** 

im = mtspeedtstrspd, trim(spdunt)//char(0)) 
im = mtstart(l) 

c 

c DO LOOP FOR REAL TIME EXECUTION 

c RUNS = DURATION/SAMPLE RATE AS SET IN MAIN MENU 

c 

do 1 1=1,runs 
c 

if(dbug)print*.’.,......call rtime’ 

c 

c ’ IF(RTIME( PERIOD, FLAG))’ THEN TIMEOUT HAS OCCURRED 

c "RTIME" PERFORMS THE DISCRETE SAMPLE INTERVAL SYNCHRONIZATION 

c 

If (rtime(period, flag)) then 
910 write(*,125) 

125 fornatC Do you wish to continue ? (y/n) - * ,\) 

read(*, *(a* ,err=910) schar 

if ( schar .eq. ’ n* .or. schar .eq. * N* ) goto 10000 
Ir = rtimeCperiod,.true.) 

end If 
c 

c Acquire sensor data, remove bias, and scale, 
c 

c GET MEASUREMENT DATA AND SCALE 

c . ...... 

c SEE NASA TMB9057, "DESCRIPTION OF THE SPACECRAFT CONTROL LABORATORY 

c EXPERIMENT (SCOLE) FACILITY" FOR INFORMATION ON THE SENSOR CHANNELS 

c 

c SEE EQUIVALENCE AT THE TOP OF THIS SUBROUTINE FOR SCOLE VARIABLE 

c ASSOCIATION WITH "ADC" 

c : -- v ■ T ' : ' • • 

c BIAS IS DETERMINED BY SUBROUTINE "SBIAS. FOR" 

c 

ig ■ getadc(0,17,adc) 
c 

do 10 kk = 1,18 

10 adcs(kk) = (adc(kk) - bias(kk)) / sfsens(kk) 

c 

c OPTIONAL - KALMAN FILTER-BASED SHUTTLE EULER ANGLE ESTIMATE 

c 

call attitf reset ) 
c 

c THE CODE BETWEEN THE ASTERIKS IS AN EXAMPLE OF A SET OF CONTROL 

c LAWS INCORPORATING THE SWITCH BOX FOR OPERATOR INTERVENTION 

c 

c USERS SHOULD PUT THEIR OWN CONTROL LAW OR ESTIMATION ALGORITHM HERE 

c 

c 

CAjb 

0“v 


\ 










c 


x * 0 
y = 0 


c Determine vibration suppression commands only if slew is not activated 
c 

dampx = 0 
dampy = 0 

if( iswtch(l) ) then 

call damp( dampx, dampy ) 
x = dampx 
y = dampy 

end If 
c 

c Determine slew commands only if vibration suppresion is not activated 
c 

slewx=0 

slewy=0 

ifC iswtch(3) ) then 

call slew( slewx, slewy ) 
x - slewx 
y * slewy 

end if 

if( iswtch(l) .and. iswtch(3) ) then 
if ( slewx+dampx .eq. 0)then 
x = 0 

else 

x = isign(l, slewx + dampx ) 

end if 

if ( slewy+dampy .eq. 0)then 
y = 0 

else 

y = isignd, slewy + dampy ) 

end if 

endif 

if( iswtch(6) )then 
x = 0 

endif 
tx = x 
ty = y 

it = thrust(x,y) 
c 

c USER SETS DACS TO ARRAY DAC VALUES 

c 

is = setdac(0,6,dac) 
c 

c DATA SENT TO STRIP CHART ( DELETE LINE IF NOT USING STRIP CHART) 

c 

c 

c USER SHOULD COMMENT OUT THIS CALL TO STRCHT IF HE/SHE DOES 

c NOT WISH TO USE THE MT8500 STRIP CHART RECORDER 

c 

call strcht 
c 

c DATA SAVED IN ARRAY FOR POST-EXPERIMENT PROCESSING 

c 

call stofild,stodat) 
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cccccccccccccccccccccccccccccccccccccc 

c 

c*******«****«•««****•***•»************************************* 

c 

1 continue 

ccccccccccc END OF REAL TIME LOOP cccc 
333 continue 
10000 continue 
c 

c USER MUST SHUT DOWN ACTUATORS FOR SAFETY REASONS 

c 

c SHUTDOWN ACTUATORS! 
c 

is = oswtch( 1,0) 
it = thrust( 0, 0 ) 
is - setdac(0,7,dacO) 
ic = cmgtor(0.,0.,0. ) 
c 

shtorr =0.0 
shtorp =0.0 
shtory =0.0 
c 
c 

c END STRIP CHART SESSION 

c (OPTIONAL - DELETE NEXT TWO LINES IF NOT USING STRIP CHART) 

im = mtstop() 
im = mtfeedO 
return 
end 
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SCOLE Makefile 


.SUFFIXES: .for .o 

FSUBS = cpmenu.for posexp.for preexp.for rtloop.for 
FMAIN = main.for 

OSUBS = cpmenu.o posexp.o preexp.o rtloop.o 
OMAIN = main.o 

INCVC=/HOLD/SCOLE/varibs.comm 
INCBC=/HOLD/SCOLE/b1k. comm 

LIBS » /lib/1ibscole.a /Iib/1lbreal.a /lib/1ibmt8500.a /lib/libmatx.a 
exec * scole 

SCOLE: $(OMAIN) $(OSUBS) $(LIBS) 

Id main.o $ (OSUBS) -lscole -lreal -lmt8500 -lmatx -lunos o=$(exec) 
delete -s *bak *lst 
echo $(exec) re-created 
.for.o : 

f77 -lrujkp -z 240 $< 

/sys/pp $*.s>$*.m 
aa -argi2s $*.m o=$*. j 
reltocoff $*.J $*.o 
delete -s $*.s S*.m $*. j 
echo $*.for compiled 
$(OMAIN): $(FMAIN) $(INCVC) $(INCBC) 
f77 -lujkmp -z 240 *V for 
/sys/pp $*.s>$*.m 
aa -argi2s $*.m o=$V j 
reltocoff $*.j S*.o 
delete -s $*.s $*.m $VJ 
echo $*.for compiled 
$(OSUBS): $( INCVC) *(INCBC) 

f77 -lrujkp -z 240 **.for 
/sys/pp $*.s>$*. m 
aa -argi2s $*.m o=$*.j 
reltocoff $*.j S*. o 
delete -s $*.s $*.m $*. J 
echo $*.for compiled 


.IGNORE : 




SILENT : 

PRECIOUS: $(FMAIN) *(FSUBS) S(OSUBS) $(OMAIN) 
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SCOLE INCLUDE FILES 


SECTION IV 
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/HOLD/SCOLE/varibs. comm 


The file /HOLD/SCOLE/varibs. comm is included in the majority of the SCOLE 
library subroutines and should be included in user supplied subroutines. 

The file contains the following common blocks and declarations. 


cccccccccccccccccccccc Sinclude file VARIBS.COMM cccccccccccccccc 


implicit real (m) 


common 

/varibs/ shurgr, 

shurgp, 

+ 

shaccx. 

shaccy, 

+ 

m2accy, 

refacx, 

+ 

refwhx. 

refwhy, 

+ 

m2why, 

tx, 

+ 

P. 

q. 

+ 

ctheta, 

stheta, 

+ 

bis 



shurgy, 

refrgr, 

refrgp, 

refrgy, 

shaccz. 

mlaccx, 

mlaccy, 

m2accx. 

refacy. 

shtorr, 

shtorp, 

shtory, 

ref whz, 

mlwhx, 

mlwhy, 

m2whx, 

ty, 

phi, 

theta, 

psi, 

r. 

fx, 

fy. 


cpsi. 

spsi, 

cphi, 

sphi, 

dbly, 

sgainy. 

dbphi, 

dbthe, 

dbpsi, 


common /gains / dgain, dblx, 

+ sgainr, sgainp, 

+ clx, cly, clz 


character*6 vdat,sdat,varib, cnum 
character*20 unit,sunit,vunit 

common /strdat/ sf (8), varib(42), unit (42), cnum(42), fsr (8) , sfmt, 
+ istosel (15), vdat (15), vunit (15), istrsel (8), sdat (8), sunit (8) 


c 

c 

c 

c 


common /consts/ pi 

common /time / period, ttime, runs 

character*10 spdunt 

common /mt8500/ strspd, spdunt 


common /chng_var/change(8) 
logical change 

cccccccccccccccccccccc END OF Sinclude file VARIBS.COMM ccccccccccc 


A description of the elements of each common block follows, 
COMMON VARIBS: 

VARIABLE DESCRIPTION 


shurgr 

shu rgp 

shurgy 

refrgr 

refrgp 

refrgy 

shaccx 

shaccy 


shuttle rate gyro roll 
shuttle rate gyro pitch 
shuttle rate gyro yaw 
reflector rate gyro roll 
reflector rate gyro pitch 
reflector rate gyro yaw 
shuttle acceleration x direction 
shuttle acceleration y direction 




!! !n!IJI!! , !!" F '!!! I 


shaccz 

shuttle acceleration z direction 



mlaccx 

mast station 1 -upper- acceleration 

X 

direct ion 

mlaccy 

mast station 1 -upper- acceleration 

y 

direction 

m2accx 

mast station 2 -lower- acceleration 

X 

direction 

m2accx 

mast station 2 -lower- acceleration 

y 

direct ion 

refacx 

reflector acceleration x direction 



refacy 

reflector acceleration y direction 



shtorr 

shuttle torque roll 



shtorp 

shuttle torque pitch 



shtory 

shuttle torque yaw 



refwhx 

reflector wheel x direction 



refwhy 

reflector wheel y direction 



refwhz 

reflector wheel z direction 



mlwhx 

mast station 1 -upper- wheel x direction 

ml why 

mast station 1 -upper- wheel y direction 

m2whx 

mast station 2 -lower- wheel x direction 

m2why 

mast station 2 -lower- wheel y direction 

tx 

thrust x direction 



ty 

thrust y direction 



phi 

bank angle 



theta 

pitch angle 



psi 

yaw angle 



P 

shuttle rotational velocity about x 



q 

shuttle rotational velocity about y 



r 

shuttle rotational velocity about z 



fx 

thruster gains x direction 



fy 

thruster gains y direction 



ctheta 

cosine theta 



stheta 

sine theta 



cpsi 

cosine psi 



spsi 

sine psi 



cphi 

cosine phi 



sphi 

sine phi 



bis 

thruster least square coefficient 




COMMON CAINS: 


VARIABLE 

DESCRIPTION 


dgaln 

dblx 

dbly 

sgainr 

sgalnp 

sgalny 

dbphi 

dbthe 

dbpsl 

clx 

cly 

clz 


thruster damping routine gain controller 

thruster damping routine dead band level x velocity 

thruster deunping routine dead band level y velocity 

thruster slew routine gain roll 

thruster slew routine gain pitch 

thruster slew routine gain yaw 

thruster slew routine roll dead band level 

thruster slew routine pitch dead band level 

thruster slew routine yaw dead band level 

collocated mast damper thruster slew routine gain roll 

collocated mast damper thruster slew routine 

gain pitch 

collocated mast damper thruster slew routine gain yav 
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COMMON STRDAT: 


VARIABLE DESCRIPTION 


sf, fsr 
varib 

unit 

cnum 

istosel 

vdat 
vuni t 
istrsel 

sdat 

vdat 


scale factors on eight dac channels sent to MT5800 
strip chart recorder 

forty-two SCOLE supplied variables from which user 
chooses fifteen to save to data file and eight to 
record on the MT8500 strip chart recorder 
forty-two units corresponding to the forty-two 
SCOLE supplied variable choices above (varib) 
diameter string nrrny consisting of elements from 
' V to 9 42 1 

array of indices of fifteen elements chosen 
from forty-two element varib array to be saved 
in a data file 

variable name (character) of fifteen varib array 
elements chosen to be saved in a data file 
units of fifteen varib array elements chosen to be 
saved in a data file 

array of indices of eight elements chosen from 
forty-two element varib array to be recorded on the 
MT8500 strip chart recorder 

variable name (character) of eight varib array 
elements chosen to be recorded on the MT8500 strip 
chart recorder 

units of eight varib array elements chosen to be 
recorded on the MT8500 strip chart recorder 


COMMON TIME: 


VARIABLE 

DESCRIPTION 

period 
tt ime 

runs 

real time loop sampling period (sec) 
real time loop duration (sec) 

number of passes in real time loop (ttime/period) 

COMMON MT8500: 


VARIABLE 

DESCRIPTION 

strspd 

spdunt 

MT8500 strip chart recorder speed 

MT8500 strip chart recorder speed units 

COMMON CHNG..VAR: 


VARIABLE 

DESCRIPTION 

change 

array to correlate variables sent to MT8500 strip 
and apropriate units 
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/HOLD/SCOLE/b1k.comm 


The file /HOLD/SCOLE/blk.comm is included in the SCOLE template main 
program located on /lib/SCOLE/TEMPLATE. 

The file contains the following common blocks and declarations. 
COMMON SCALE: 

VARIABLE DESCRIPTION 


sfsrr 

scale 

factor 

sfsrp 

scale 

factor 

sfsry 

scale 

factor 

sf rrr 

scale 

factor 

sfrrp 

scale 

factor 

sfrry 

scale 

factor 

sfsax 

scale 

factor 

sfsay 

scale 

factor 

sfsaz 

scale 

factor 

sfmlx 

scale 

factor 

sfmly 

scale 

factor 

sf m2x 

scale 

factor 

sf m2y 

scale 

factor 

sfrax 

scale 

factor 

sfray 

scale 

factor 

sfstr 

scale 

factor 

sfstp 

scale 

f actor 

sfsty 

scale 

factor 


The values of the above listed variat 
elements of array sfsens in blk.comm. 


shuttle rate gyro roll 

shuttle rate gyro pitch 

shuttle rate gyro yaw 

reflector rate gyro roll 

reflector rate gyro pitch 

reflector rate gyro yaw 

shuttle acceleration x 

shuttle acceleration y 

shuttle acceleration z 

mast station 1 -upper- acceleration x 

mast station 1 -upper- acceleration y 

mast station 2 -lower- acceleration x 

mast station 2 -lower- acceleration y 

reflector acceleration x 

reflector acceleration y 

shuttle torque roll 

shuttle torque pitch 

shuttle torque yaw 

►les correspond to the eighteen 


Scale factors have units of unity/engineering unit. 

The units for the rate sensors are units/radian/second. 
The units for the accelerometers are units/inch/second^2. 


COMMON CMGSET: (CMC software In process of being updated) 
VARIABLE DESCRIPTION 


rpml cmg rpm - default set to 3000 

rpm3 cmg rpm - default set to 3000 


sfl 

scale 

sf2 

scale 

sf3 

scale 

pd 

servo 

tl 

servo 

td 

servo 

b2 

servo 


factor cmg - default set to 10. 

factor cmg - default set to 10. 

factor cmg - default set to 10. 

gain prop. - default set to 60. 

gain integ - default set to .0005 

gain deriv - default set to 9000. 

gain damp - default set to ,9925 
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The values of array sfsens are set In data statement sfsens. 

data sfsens / 7.860, 11.2351, 4.7080, .1975, 1.15, 1.15, 
+ .00518, .00518, .00259, 

+ 4*2.591e-3, 2*5.18e-4, 3*1. / 


The forty-two character elements of array varlb are defined In data 
statement varlb. 


data 

varlb / 





’refrgy’ 

+ 

’shurgr’, 

’shurgp’, 

’shurgy’, 

’refrgr’, 

’refrgp’, 

+ 

’shaccx’, 

’shaccy’, 

’shaccz’, 

’mlaccx’, 

’mlaccy’, 

’m2accx’ 

+ 

’m2accy’, 

’refacx’, 

’refacy’, 

’shtorr’, 

’shtorp’, 

’shtory’ 

+ 

’refwhx’, 

’refwhy’, 

’refwhz*, 

’mlwhx’ , 

’mlwhy’, 

’ m2whx’, 

+ 

’m2why’, 

’tx’ , 

’ty\ 

’phi’. 

’theta’ , 

’psl’, 

+ 

* « » 

P * 

» _f 

Q * 

•r\ 

’fx’ , 

’fy\ 

’sphi’, 

+ 

+ 

’ctheta’, 
’bis’ / 

* stheta p , 

'cps1’, 

’sps1’ , 

’cphi’ , 


The forty-two character units of array unit are defined in data statement 
unit. 


data unit / 

6* * rad/s’,9*’ln/s**2’,10*’ ft-lbs’ , 2*’ lbs’ , 
3*’rad’,3*’rad/s’,2*’N',6*’units’,’units’/ 


The forty-two character elements of array cnum are defined in data statement 
cnum. 


+ 

+ 

+ 


data cnum /’1’,’2’,’3’,’4’,’5’, 
’ll’,’12’,’13’,’14’,’15’,’16’ , ’ 
’21’,’22’,’23’,*24’,’25’,’26’,’ 
’31’,’32’,’33’,’34’,’35’, ’36’ , ’ 


’6’,’7’,*8’,’9’,’10’ 
17’,’18’,’19’,’20’, 
27’,’28’,’29’, ’30’ , 
37’,*38’,’39’ , ’40’ , ’ 


41’,’42’/ 
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SCOLE LIBRARY FILES 


SECTION V 
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SUBROUTINE ATTIT 


CALL: call attit(reset) 

logical reset 

Kinematic observer used to determine the roll and pitch attitude 
of the shuttle from the shuttle accelerometers and rate sensor. 
Reset is input as true for the first call to attit. 


SUBROUTINE CHAN 


CALL: call chan 

User prompted for number of channels to be printed on the MT8500 
strip chart recorder (1-8 channels). 


SUBROUTINE CMGMEN 

(CMG software In process of being updated) 
CALL: call cmgmen 

User presented a cmg set-up menu as follows: 

CMG SETUP MENU 

1. INPUT RPM 

2. TAKE ADC TARE 

3. TEST STEERING LAW CONTROL 

4. TEST GIMBAL COMMAND CONTROL 

5. CHANGE SERVO GAINS 

6. CHANGE INPUT SCALE FACTOR 

7. CHECK CMG SENSOR OUTPUTS 

8. RETURN TO PREVIOUS MENU [DEFAULT] 

sets rpml, rpm3 
sets array adct(12) 
calls subroutine cmgtst 
calls subroutine gimbal 
sets pk, ti, td, b2 
sets sfl, sf2, sf3 
calls subroutine sensor 


SUBROUTINE CMGTOR 

(CMG software in process of being updated) 

CALL: img = cmgtor(xcmd,ycmd,zcmd) 

FORTRAN callable C routine 
integer*4 cmgtor, img 

Three axis shuttle torque commands (x, y, z) sent to cmg controller. 


upnon 

Option 

Option 

Option 

Option 

Option 

ODtion 
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SUBROUTINE CMGTST 

(CMC software in process of being updated) 
CALL: call cmgtst 

Cmgs tested. 


SUBROUTINE COLINT 

CALL: call colint 

User presented the following menu to accept the gain and dead-band levels 
for the thruster slew routine as shown in the documentation for the real 
time loop - rtloop. 


I •• COLOCATED MAST DAMPER MENU ** I 


I 1) SET THE GAIN FOR ROLL I, clx 
I 2) SET THE GAIN FOR PITCH I, cly 
I 3) SET THE GAIN FOR YAW I, clz 
I 4) RETURN TO MAIN MENU [DEFAULT] I 


Values of clx, cly, clz passed in common block gains in include file 
/HOLD/SCOLE/varibs.comm. 


SUBROUTINE CSCALE 


CALL: call cscale 

User prompted for scale factors by which to scale data being sent to 
the MT8500 strip chart recorder. 


SUBROUTINE DAMINT 

CALL: CALLED FROM SUBROUTINE CPMENU (call damint) 

User presented following menu to accept gain and dead band levels for 
thruster damping routine. 

Variables dgain, dblx, dbly, set and passed in common block gains in' 
include file /HOLD/SCOLE/varibs.comm 


i ** THRUSTER DAMPING MENU ** I 


I 1) SET THE GAIN FOR THE CONTROLLER I, dgain 

I 2) SET DEAD BAND LEVEL OF X-VELOCITY I, dblx 

I 3) SET DEAD BAND LEVEL OF Y-VELOCITY I, dbly 

I 4) RETURN TO MAIN MENU [DEFAULT] I 
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SUBROUTINE DITIME 


CALL: call ditIme(period, ttiroe) 

period: sample period - set in subroutine ditime 
ttime : run duration - set in subroutine ditime 

User inputs test duration time (ttime) and sample period (period). 


SUBROUTINE ERTRAP 

CALL: call ertrap(ercode,ertype,mode) 

ercode - passed in value of ios in unsuccessful file open 
(ie. open(..,iostat=ios,..)) 
ertype - passed in as 1 to execute subroutine code 
passed in as 0 to abort program 
mode - passed In as 0 to continue program 
passed In as 1 to abort program 

Subroutine called after a file open as follows: 
cal1 ertrap(Ios,1,0). 


SUBROUTINE GET AC 

CALL: CALLED FROM SUBROUTINE LSQANG (call getac(ax, ay, az)) 

Subroutine gets data from accelerometers, 
ax set based on adc6 
ay set based on adc7 
az set based on adc8 


SUBROUTINE GETTOR 

CALL: call gettor 

Shuttle torque commands gotten from joystick. 


SUBROUTINE GIMBAL 

(CMG software In process of being updated) 

CALL: call gimbal(test) 

logical test 

Four torque commands sent to gimbals. 

If switch 2 is set on the switch box the subroutine code is not executed. 
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SUBROUTINE ISELEK(MAXP,MODE,IANS,NANS,LHNO,IDEF) 


CALL: call iselek(maxp,mode,ians,nans,lhno,idef) 

maxp, mode, Ians, lhno and ldef sire integers 

nans should be declared character*l In the calling subroutine 

The user chooses a menu selection line based on the mode chosen. 

The following modes are available. 

MODE DESCRIPTION 


1 1-maxp, H, M, Q or * 

2 1-maxp, H, Q or * 

3 1-maxp, H, M or * 

4 1-maxp, H or * 

5 1-maxp, Q or * 

6 1-maxp, H, M, Q 

7 1-maxp, H, Q 

8 1-maxp, H, M 

9 1-maxp, H 

10 1-maxp, Q 

11 1-maxp, M 

12 1-maxp 

13 Y or N 

14 [Y] or N 

15 [Y] or [N] 

16 H, Y or N 

17 H, [Y] or N 

18 H, Y or [N] 

19 [MORE] enter <cr> to continue . . . 

20 Enter <cr> to continue . . . 

Where, 

maxp = maximum range of selection. 

[xx] = default value enclosed in brackets. 

H = display help message; returns ’M’ in nans and 0 in ians. 

M = returns’M’ in nans. 

Q = gives user option to abort program; if not confirmed 
an ’M’ will be return in nans and 0 in ians. 

* = return ’A’ in nans and 0 in ians. 

<cr> = return ’C’ in nans and 0 in ians. 

IANS paissed out of iselek equal to zero if user inputs invalid selection 
number entered else equal to valid selection number input. 

NANS passed out of iselek as follows: 

C modes 19, 20 

Y modes 14, 17; all yes choices with modes 13, 16, 18 

N modes 15, 18; all no choices with modes 13, 16, 17 

D choices 1 - maxp 

M choices H, M, Q 

A choice * 
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If nans equals M the user should direct program flow back to the menu In 
question. (ie.. If(nans. eq.’M’)go to xxx). 

IDEF is Input equal to MAXP. 

LHNO is input with modes incorporating the H option for on-line 
help as follows: 

0 user does not know what standard SCOLE subroutines 

are documented on-line; menu presented subroutines 
for which on-line help Is available 

1 user knows which subroutine for which help is 
desired; user inputs name In argument hlpfil 

2 user desires help for a routine that is not 
one of the SCOLE standard routines; a help file 
for the subroutine Is located on the SCOLE working 
directory; user inputs file name prefix in argument 
hlpfil (file extension Is ".doc"). 

In the calling subroutine, LHNO should be set to 1 or 2 and hlpfil should be 
declared character*20 and be set to the name of the calling subroutine. If 
lhno Is 2 then the user must have created a help file on the current working 
directory with prefix equal to the subroutine name and suffix ’.doc’. 

The common block info In subroutine iselek and lhelp passes parameters 
procse and hlpfil. 


SUBROUTINE LBMENU 

CALL: CALLED FROM SUBROUTINE STMENU (call lbmenu) 

USERS SHOULD NOT NEED TO CALL THIS SUBROUTINE 


User presented a menu els follows with options to choose the MT8500 strip 
chart recorder labels. 


I ** CHOOSE LABELS ** I 

I 1) USE PROGRAM VARIABLE NAMES AS DESCRIPTORS. I 
I 2) USE FULL DESCRIPTION OF PROGRAM VARIABLES. I 
I 3) CHOOSE FILE TO USE FOR LABEL DATA. I 
I 4) CREATE CUSTOM LABELS FILE. I 
I 5) RETURN TO PREVIOUS MENU. (DEFAULT) I 


SUBROUTINE LHELP(HLPFIL,LHNO) 

CALL: call lhelp(hlpfi1,lhno) 

This subroutine will allow the user to obtain on-line help for 
the standard SCOLE subroutines smd/or user supplied subroutines. 
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LHNO Is integer input bis follows: 


0 user does not know what standard SCOLE subroutines 

sire documented on-line; menu presented of help 
available for subroutines 

1 user knows which subroutine for which help is 
desired; user Inputs name In argument hlpfil 

2 user desires help for a routine that is not 
one of the SCOLE standard routines; a help file 
for the subroutine is located on the SCOLE working 
directory; user Inputs file name prefix in argument 
hlpfil (file extension Is ".doc"). 

Hlpfil must be declared chanacter*20 in the calling subroutine and set to the 
name of the subroutine for which help is desired if lhno is set to 1 and set 
to the prefix of the name of the help file if lhno is set to 2. 

The common block info as found in subroutines Iselek and lhelp passes 
parameters procse and hlpfil. 

common /Info/ procse, hlpfil 

This subroutine is called from subroutine iselek when modes including the H 
menu option are chosen. 

Users can include a call to lhelp from a menu if iselek is not used. 


SUBROUTINE LSQANG 

CALL: call lsqang(ctheta,stheta,cpsi,spsi,cphi,sphi) 


Subroutine determines leaist square angle estimation. 

Arguments set in subroutine lsqang an follows 
ctheta - cosine theta 
stheta - sine theta 

cpsi - cosine psi 

spsi - sine psi 

cphi - cosine phi 

sphi - sine phi. 


SUROUTINE READF 


CALL: call readf 

SCOLE formatted direct access files gains.dat and select.set, copied onto 
the user’s account from /lib/SCOLE/TEMPLATE, are read. 

The section of the code that reads the files follows. 
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c Read from "gains.dat" 
c 

open(5,err=9996, iostat=ios,form=’formatted’,access=’direct’ , 

+ recl=80,file=’gains.dat’,status=’old’) 

if(dbug)print*,’Subroutine readf -file gains.dat opened’ 
read(5,25,rec=2)dgain,dblx,dbly 

25 format(3(el2.5,lx)) 
read(5,26,rec=l)ttime,period 

26 format(2(el2. 5,lx)) 

read(5,27, rec=3)sgainr,sgainp,sgainy,dbphi,dbthe,dbpsi 

27 format(6(el2. 5, lx)) 

read(5,39,rec=4)clx, cly,clz 

39 format(3(el2. 5, lx)) 

close(5) 
c 

c Read from "select.set" 
c 

open(6,err=9996, iostat=ios,form=’formatted’,access=’direct’ 

+ recl=l12,file=’ select. set’,status=’ old’ ) 
read(6,28, rec=l)(fsr(i),i=l,4) 
read(6,28,rec=2)(fsr(i),i=5,8) 

28 format(4(el2. 5, lx)) 

read(6,2000,rec=3)strspd,spdunt 
2000 format(f10. 2,lx,alO) 

read(6,29,rec=4)(istosel(i),i=l,15) 

29 format(15(i3, lx)) 

read(6,30, rec=5)(vdat(i),1=1,15) 

301 format(3(a20,lx)) 

30 format(15(a6,lx)) 

read(6,301,rec=6)(vunit(i),i=l,3) 
read(6,301, rec=7)(vunit(i),i=4,6) 
read(6,301,rec=8)(vunit(i),i=7,9) 
read(6,301,rec=9)(vunit(i),i=10,12) 
read(6,301,rec=10)(vunit(i),i=13,15) 
read(6,31,rec=11)(istrsel(i),i=l,8) 
read(6,32,rec=12)(sdat(l),1=1,8) 
read(6,301,rec=13)(sunit(i),1=1,3) 
read(6,301, rec=14)(sunit(i),i=4,6) 
read(6,304, rec=15)(sunit(i),i=7,8) 

304 format(2(a20,lx)) 

close(6) 
c 

31 format(8(13,lx)) 

32 format(8(a6,lx)) 

321 format(8(a20,lx)) 

To alter values in these files, files should be opened as shown and records 
read and written with the indicated formats. 
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VARIABLE 


USE 


SET IN SUBROUTINE COMMON BLOCK IN 

/HOLD/SCOLE/VARIBS.COMM 


dgain 

controller gain 

damint.for 

gains 

dblx 

dead band level x 
velocity 

damint. for 

gains 

dbly 

dead band level y 
velocity 

damint.for 

gains 

ttime 

run duration 

ditime.for 

parameter 

period 

sample period 

ditime.for 

parameter 

sgainr 

gain roll thruster 
slew 

slwint.for 

gains 

sgainp 

gain pitch thruster 
slew 

slwint.for 

gains 

sgainy 

gain yaw thruster 
slew 

slwint. for 

gains 

dbphi 

roll dead band level 
thruster slew 

slwint.for 

gains 

dbthe 

pitch dead band level 
thruster slew 

slwint.for 

gains 

dbpsi 

yaw dead band level 
thruster slew 

slwint.for 

gains 

clx 

gain roll colocated 
mast damper 

colint.for 

gains 

cly 

gain pitch colocated 
mast damper 

colint.for 

gains 

clz 

gain yaw colocated 
mast damper 

colint.for 

gains 

fsr 

MT8500 scale factors 

cscale.for 

strdat 

strspd 

MT8500 speed 

spdset.for 

strdat 

spdunt 

KT8500 speed units 

spdset.for 

strdat 

Istosel 

array of indices of 

15 variables of 42 
SCOLE variables to be 
saved to a data file 

stoset.for 

strdat 

vdat 

array of 15 variable 
names selected to be 
saved to data file 

stoset.for 

strdat 

vunit 

array of 15 variable 
units selected to be 
saved to data file 

stoset.for 

strdat 

istrsel 

array of indices of 

8 variables of 42 

SCOLE variables to be 
sent to the MTT8500 

strset.for 

strdat 

sdat 

array of 8 variable 
names selected to be 
sent to the MT8500 

strset.for 

strdat 

sunit 

array of 8 variable 
units selected to be 
sent to the MT8500 

strset.for 

strdat 
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SUBROUTINE SBIAS 


CALL: call sbias(tau) 

tau - sample period for 20 second averaging 

Subroutine determines sensor bias by averaging over 20 seconds. 
Bias values passed in array bias in common store and written to 
file calibration.dat. Bias(l) bias for adc(l), bias(2) bias for 
adc(2), etc. 


SUBROUTINE SLWINT 

CALL: CALLED FROM SUBROUTINE CPMENU (call slwint) 

User presented the following menu to set the gain and dead-band levels 
for the thruster slew routine. 


I •• THRUSTER SLEW MENU ** I 


1) SET THE GAIN FOR ROLL i, sgainr 

2) SET THE GAIN FOR PITCH I, sgainp 

3) SET THE GAIN FOR YAW I, sgainy 

4) SET ROLL DEAD BAND LEVEL I. dbphi 

5) SET PITCH DEAD BAND LEVEL I, dbthe 

6) SET YAW DEAD BAND LEVEL I. dbpsi 


I 7) RETURN TO MAIN MENU [DEFAULT] 


Variables passed in common gains in include file /HOLD/SCOLE/varibs.comm 


SUBROUTINE SPDSET 


CALL: call spdset 

User given option to change MTS500 strip chart recorder speed. 


SUBROUTINE STOFIL 

CALL: call stofi1(n,stodat) 

n - pass number in real time loop 

stodat - array of data points to be stored; dimensioned 
15 by runs where 15 is the number of variables 
to be saved and runs is the number of passes 
in the real time loop calculated to be duration/period. 

Subroutine called in real time loop for data point storage. 
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SUBROUTINE STOSET 


CALL: call stoset 


User selects program variables (15) to be output to data file 
from the 42 SCOLE variables saved In common block varibs in 
1nc1ude /HOLD/SCOLE/varibs.comm 


common /varibs/ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


shurgr, shurgp 
shaccx, shaccy 
m2accy, refacx 
refwhx, refwhy 
m2why, tx, 

P> q» 
ctheta, stheta 
bis 


shurgy, refrgr 
shaccz, mlaccx 
refacy, shtorr 
refwhz, mlwhx, 
ty, phi, 

r, fx, 

cpsi, spsi. 


refrgp. 

refrgy, 

mlaccy, 

m2accx, 

shtorp, 

shtory, 

mlwhy, 

m2whx, 

theta. 

psi, 

fy. 


cphi, 

sphi, 


SUBROUTINE STRCHT 


CALL: call strcht 

The values of the eight SCOLE variables selected in subroutine strset sent 
to the KT8500 strip chart recorder. These values are scaled by the scale 
factors set in subroutine cscale before being sent to the strip chart 
recorder. Subroutine called in a real time loop. 


SUBROUTINE STRSET 


CALL: call strset 


User selects program variables (8) to be output to the MT8500 strip chart 
recorder from the 42 SCOLE variables saved in common block varibs in 
include /HOLD/SCOLE/varibs. comm. 


common /varibs/ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


shurgr, shurgp 
shaccx, shaccy 
m2accy, refacx 
refwhx, refwhy 
m2why, tx, 
p. 

ctheta, stheta, 
bis 


shurgy, refrgr 
shaccz, mlaccx 
refacy, shtorr 
refwhz, mlwhx, 
ty, phi, 

r, fx, 

cpsi, spsi. 


refrgp, 

mlaccy, 

shtorp, 

mlwhy, 

theta, 

fy. 

cphi, 


refrgy, 
m2accx, 
shtory, 
m2whx, 
psi, * 

sphi, 


SUBROUTINE TSCOLE 
CALL: call tscoleCnm, nd,ny) 

The current date is output in integer form: month - run, day - nd, and 
year - ny. 
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SUBROUTINE UWRITE 


CALL: call uwrite(nsamp,tau,stodat) 

nsamp - number of samples, precalculated to be duration / period 

tau - sample period 

stodat - array of data points stored 

User presented the following menu of options to save run data to a file 
after experiment complete. 

*»*#**••****#*•*•••»**»•#*•••»*»••••*••* 

DATA OUTPUT MENU 

•**•**•*•»•*******•*******••••*•**••*••» 

1. CRDS DATA FILE FOR ANALYSIS 

2. MATRIX-X DATA FILE 

3. MATLAB DATA FILE 
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MT8500 STRIP CHART RECORDER LIBRARY ROUTINES CALLED BY USERS 


SECTION VI 
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SUBROUTINE MTCHAN 


NAME: mtchan(nchan) 

DESCRIPTION: This subroutine sets the number of channels to be 
printed from the MTB500. 

CALL: 

FORTRAN: 

lnteger*4 mtchan, kepone 
Integer nchan 
kepone = mtchan(6) 

C : 

integer nchan; 
mtchan(nchan); 

INPUT ARGUMENTS: 

nchan Integer 1-8 = number of channels to be printed. 

OUTPUT ARGUMENTS: 

return: 1 = successful call 

0 = unsuccessful call 

REQUIRES: 

Nothing 

LIBRARY: /I lb/1lbmt8500.a 


SUBROUTINE MTDATE 


NAME: mtdateO 

DESCRIPTION: This subroutine sets the time and date on the MT8G00 

chart recorder. 

(This subroutine Is not normally called by applications 
date is automatically set by other routines.) 


CALL: FORTRAN: 


c 

IntegerM mtdate, im 
lm = mtdateO 

INPUT ARGUMENTS: 

None 

OUTPUT ARGUMENTS: 

mtdatef); 


return: 1 = successful call 

0 = unsuccessful call 
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REQUIRES: 


Nothing 


USES: 

LIBRARY: 


cvtlme 

/llb/1ibmt8500.a 


SUBROUTINE MTFEED 


NAME: mtfeedO 

DESCRIPTION: This subroutine advances the paper In the MT8500 chart 
recorder to the next top-of-form mark (form feed). 


CALL: 

FORTRAN: 

lnteger*4 mtfeed, J 
J=mtfeed() 

C : 

mtfeed(); 

INPUT ARGUMENTS: 

None 

OUTPUT ARGUMENTS: 

return : 1 (no check for a successful call) 

REQUIRES: 

Nothing 

LIBRARY: /llb/1ibmt8500.a 


SUBROUTINE MTFILE 


NAME: mtf1le(’fname’) 


DESCRIPTION: This subroutine sets the channel labels on the MT8500 

chart recorder according to the filename supplied In fname. 


CALL: 

FORTRAN: 

lntegerM mtf lie, Igor 

character*(*) fname 

Igor = mtf1le("llabel.doc") 

C : 

chart] fname; 
mtflie (fname); 

INPUT ARGUMENTS: 

fname string pointing to sequential data file of labels 
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OUTPUT ARGUMENTS: 

return: 1 = successful call 

0 = unsuccessful call 


REQUIRES: 

All labels must end in a <CR>, be written sequentially, 
and be < 100 characters in length. 

LIBRARY: /Iib/1ibmt8500.a 

SUBROUTINE MTINIT 


NAME: mtinit(nchan) 

DESCRIPTION: This subroutine initializes the MT8500 chart recorder for 

use from its "Ready" state. Communication lines are 
opened between the chart recorder and the CRDS on /dev/tty4. 
The MT8500 is initialized to the "Realtime" state and 
passed the number of channels to be printed. Finally, date 
and time are passed to the recorder via the mtdate subroutine. 

This subroutine is NOT called by user programs! 

CALL: 

FORTRAN: 

integer*4 mtinit, i 
integer nchan 
i=mtinit(nchan) 

C : 

int nchan; 
mtinit(nchan); 

INPUT ARGUMENTS: 

nchan {integer 1-8} = number of channels to be printed. 

OUTPUT ARGUMENTS: 

return: 1 = successful call 

0 = unsuccessful call 

REQUIRES: 


Nothing 


LIBRARY: 


/lib/1ibmt8500. a 




SUBROUTINE MTRAY 


NAME: mt ray(a_name) 

DESCRIPTION: This subroutine sets the eight channel labels on the sMT8500 
chart recorder according to the array of labels in a_name. 


CALL: 

FORTRAN: 

integer*4 mtray, lupe 
dimension a_name(8) 
character a_name*100 
lupe = mtray(a_name) 

INPUT ARGUMENTS: 

a_name = array of 0 to 8 label strings 

OUTPUT ARGUMENTS: 

return: 1 = successful call 

0 = unsuccessful call 

REQUIRES: 

Each label string must end with a "0" said be < 100 
characters in length. 

LIBRARY: /Iib/1ibmt8500.a 


SUBROUTINE MTRLS 


NAME: mtrlsO 

DESCRIPTION: This subroutine releases the MT8500 chart recorder 

from the real-time mode and leaves it in a ready state. 


CALL: 

FORTRAN: 

integer*4 mtrls, 1 
1 = mtrlsO 

C : 

1 = mtrlsO; 

INPUT ARGUMENTS: 

None 


OUTPUT ARGUMENTS: 

return value: 1 = successful call 

0 = unsuccessful call 


REQUIRES: 

Nothing 

LIBRARY: /lib/1ibmt8500.a 
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SUBROUTINE KTSTOP 


NAME: mtstopC) 

DESCRIPTION: Stop waveform printing from the MT8500 chart recorder 


CALL: 

FORTRAN: 

integerM mtstop, Jove 
Jove = mtstopC) 

C : 

i = mtstopC); 

INPUT ARGUMENTS: 

None 


OUTPUT ARGUMENTS: 

return value: 1 = successful call 

0 = unsuccessful call 


REQUIRES: 

Nothing 

LIBRARY: /Iib/1ibmt8500.a 
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HATH LIBRARY ROUTINES 


SECTION VII 
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SUBROUTINE ADD 


NAME: ADD() 

DESCRIPTION: This subroutine performs matrix addition of Input matrices [A] 
and lB]; sum matrix output In [C]. 

CALL: call add(a,na,b,nb,c,nc,m,n) 

INPUT ARGUMENTS: 

A,B matrices packed by columns In one dimensional arrays 
NA t NB dimension vectors of [A] and [B] respectively 
NA(1)=NB(1)=number of rows of [A] and [B] 
NA(2)=NB(2)=number of columns of [A] and [B] 

M number of rows of [A] and [B] 

N number of columns of [A] and [B] 

OUTPUT ARGUMENTS: 

C [C]=[A]+[B]; sum matrix packed by columns in one 

dimensional array 
NC dimension vector of [C] 

NC(1)=NA(1)=NB(1)=M 
NC(2)=NA(2)=NB(2)=N 

REQUIRES: 

NA(1)=NB(1)=NC(1)=M 
NA(2)=NB(2)=NC(2)=N 

LIBRARY: 

/llb/1ibmatx. a 


SUBROUTINE DISCRT 


NAME: DISCRT() 

DESCRIPTION: This subroutine computes exp(a*tau) and its integral using the 
series expansion definition for up to nterms. 

CALL: call discrt(a,na,phi,nphi,gam,ngam,n, tau,nterms) 


INPUT ARGUMENTS: 

A matrix packed by columns in one dimensional array 

NA dimension vector of fA] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

N number of rows and columns of [A] 

NTERMS number of terms to expand 


OUTPUT ARGUMENTS: 

PHI matrix packed by columns in one dimensional array 
at 

[PHI]=e 
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GAM matrix packed by columns In one dimensional array 

as 

[GAM]=lntegral e ds from 0 to tau 
NPHI, dimension vectors [PHI] and [GAM] respectively 
NGAM 

NPHI(1)=NGAM(1)=NA(1)=N 
NPHI(2)=NGAM(2)=NA(2)=N 

REQUIRES: 

NPHI(1)=NGAM(1)=NA(1)=N 
NPHI(2)=NGAM(2)=NA(2)=N 

LIBRARY: /llb/1lbmatx. a 


FUNCTION DOT 


NAME: DOTH 

DESCRIPTION: This function calculates the dot product of two vectors, 
a and b. 

CALL: caldot = dot (a,b,n) 

INPUT ARGUMENTS: 

A,B column vectors 

N dimension of vector [A] 

OUTPUT ARGUMENTS: 

DOT dot product of Input vectors 

REQUIRES: 

Nothing 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE EQUATE 


NAME: EQUATE() 


DESCRIPTION: This subroutine equates two matrices; [B]=[A]. 
CALL: call equate(a,na,b,nb,m,n) 


INPUT ARGUMENTS: 

A matrix packed by columns In one dimensional array 

NA dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

M number of rows of [A] 

N number of columns of [A] 


OUTPUT ARGUMENTS: 

B [B]=[A]; equivalent matrix packed by columns in one 
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NB 


dimensional array 
dimension vector of [B] 
NB(1)=M 
NB(2)=N 


REQUIRES: 

NA(1)=NB(1)=M 
NA(2)=NB(2)=N 

LIBRARY: /llb/1ibmatx.a 


FUNCTION FMAG 

NAME: FMAG() 

DESCRIPTION: This function finds the magnitude of a vector A. 

CALL: calmag * fmag(a,n) 

INPUT ARUGMENTS: 

A column vector 

N dimension of vector A 

OUTPUT ARUGMENTS: 

fmag magnitude of vector A 

REQUIRES: 

Nothing 

LIBRARY: /llb/1Ibmatx.a 


SUBROUTINE IDENT 


NAME: IDENT() 


DESCRIPTION: This subroutine creates an Identity matrix; [A]=[I]. 


CALL: call ldent(a,na,n) 


INPUT ARGUMENTS: 

A matrix packed b y columns In one dimensional array 

NA dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)“number of columns of [A] 

N number of rows and columns of [A] 


OUTPUT ARGUMENTS: 

A [A]=[I]; Identity matrix packed by columns in one 

dimensional array 


REQUIRES: 


NA(1)=NA(2)=N 
call to NULL 


680 





LIBRARY: 


/lib/1lbmatx. a 


SUBROUTINE MULT 


NAME: MULTO 

DESCRIPTION: Thin mibroutlnn performs matrix mult I pi lent Ion oT Input 
matrices [A] and [B]; product matrix output In [Cl. 

CALL: call mult(a,na,b,nb,c,nc,1,m,n) 

INPUT ARGUMENTS: 

A,B matrices packed by columns In one dimensional arrays 
NA,NB dimension vectors of [A] and [B] respectively 
NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

NB(l)=number of rows of [B] 

NB(2)=number of columns of [B] 

L number of rows of [A] and [C] 

M number of columns of [A] and rows of [B] 

N number of columns of [B] and [C] 

OUTPUT ARGUMENTS: 

C [C]=[A]*[B]; product matrix packed by columns in 

one dimensional array 
NC dimension vector of [C] 

NC(1)=NA(1)=L 
NC(2)=NB(2)=N 

REQUIRES: 

NA(1)=NC(1)=L 
NA(2)=NB(1)=M 
NB(2)=NC(2)=N 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE MULT3 


NAME: MULT3() 


DESCRIPTION: This subroutine performs matrix multiplication of Input 
matrices [A],[B], and [C]; product matrix output In [D]. 

CALL: call mult3(a, b, c, d, ar, ac, be, cc) 


INPUT ARGUMENTS: 

A, B, C 

AR 

AC 

BC 

CC 


matrices packed by columns In one dimensional arrays 

number of rows of [A] and [D] 

number of columns of [A] and rows of [B] 

number of columns of [B] and rows of [C] 

number of columns of [C] and [D] 
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OUTPUT ARGUMENTS: 

D [D]=[A]*[B]*[C]; product matrix packed by columns in 

one dimensional array 

REQUIRES: 

Nothing 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE NEC. 

NAME: NEG() 

DESCRIPTION: This subroutine renders the negative of a matrix ; [B]=-l.0*[A]. 

CALL: call neg(a,na,b,nb,m,n) 

INPUT ARGUMENTS: 

A 


M 
N 

OUTPUT ARGUMENTS: 

B tB]=—1.0*[A]; scaled matrix packed by columns In one 

dimensional array 
NB dimension vector of [B] 

NB(1)=M 
NB(2)=N 

REQUIRES: 

NA(1)=NB(1)=M 
NA(2)=NB(2)=N 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE NULL 

NAME: NULL() 

DESCRIPTION: This subroutine nulls a matrix; [A]=[0]. 

CALL: call nul1(a,na,n,m) 

INPUT ARGUMENTS: 

A matrix packed by columns in one dimensional array 

NA dimension vector of [A] 

NA(1)=number of rows of [A] 

NA(2)=number of columns of [A] 

N number of rows of [A] 

M number of columns of {A] 


matrix packed by columns In one dimensional array 
dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 
number of rows of [A] 
number of columns of [A] 
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OUTPUT ARGUMENTS: 

A [A]=[0]; null matrix packed by columns In one 

dimensional array 

REQUIRES: 

NA(1)=N 
NA(2)=M 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE PRNT 

NAME: PRNTO 

DESCRIPTION: This subroutine prints a matrix for output to a printer; 

132 columns are used; matrix Is printed 1 row by 9 columns. 

CALL: call prnt(a,na,nrow,mcol,name) 

INPUT ARGUMENTS: 

A matrix packed by columns In one dimensional array 

NA dimension vector of lA] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

NROW number of rows of [A] 

MCOL number of columns of [A] 

NAME descriptor of matrix to be printed; may be passed as 
a string or as a variable declared as a character 

OUTPUT ARGUMENTS: 

NONE 

REQUIRES: 

NA(1)=NROW 
NA(2)=MC0L 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE PRNT1 

NAME: PRNT1() 

DESCRIPTION: This subroutine prints a matrix for output to a terminal. 
CALL: call prntl(a,na,n,m,name) 

INPUT ARGUMENTS: 

A matrix packed by columns in one dimensional array 

NA dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

N number of rows of [A] 
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M number of columns of lA] 

NAME descriptor of matrix to be printed; may be passed as 
a string or as a variable declared as a character 

OUTPUT ARGUMENTS: 

NONE 

REQUIRES: 

NA(1)=N 
NA(2)=M 

LIBRARY: /llb/1ibmatx. a 


SUBROUTINE PRNTF 

NAME: PRNTF() 

DESCRIPTION: This subroutine prints a matrix to a file. 

CALL: call prntf(a,na,nrow,mcol,name) 

INPUT ARGUMENTS: 

A matrix packed by columns In one dimensional array 

NA dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

NROW number of rows of [AJ 
MCOL number of columns of [A] 

NAME descriptor of matrix to be printed; may be passed as 
a string or as a variable declared as a character 

OUTPUT ARGUMENTS: 

NONE 

REQUIRES: 

NA(1)=NROW 
NA(2)=MC0L 

Output file must be opened and closed in calling program; unit 
number 6 is used. 

LIBRARY: /llb/1ibmatx. a 


SUBROUTINE SCALE 

NAME: SCALE!) 

DESCRIPTION: This subroutine scales a matrix by a scalar; [B]=c[A]. 
CALL: call scale(a,na,b,nb,m,n,c) 

INPUT ARGUMENTS: 

A matrix packed by columns in one dimensional array 

NA dimension vector of [A] 
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NA(l)=number of rows of [Aj 
NA(2)=number of columns of [A] 

M number of rows of [A] 

N number of columns of lA] 

C scale factor passed by value or as variable 


OUTPUT ARGUMENTS: 

B [B]=c[A]; scaled matrix packed by columns In one 

dimensional array 
NB dimension vector of IB] 

NB(1)=M 
NB(2)=N 


REQUIRES: 

NA(1)=NB(1)=M 
NA(2)=NB(2)=N 

LIBRARY: /llb/1lbmatx.a 


SUBROUTINE SUBTRACT 


NAME: SUBTRACT() 

DESCRIPTION: This subroutine performs matrix subtraction of Input 

matrices [A] and IB]; difference matrix output In tC]. 

CALL: call subtract(a,na,b,nb,c,nc,m,n) 

INPUT ARGUMENTS: 

A,B matrices packed by columns In one dimensional arrays 
NA,NB dimension vectors of [A] and [B] respectively 
NA(1)=NB(1)=number of rows of [A] and [B] 
NA(2)=NB(2)=number of columns of [A] and [B] 

M number of rows of [A] and [B] 

N number of columns of [A] and [B] 

OUTPUT ARGUMENTS: 

C (C]=[A]-[B]; difference matrix packed by columns in 

one dimensional array 
NC dimension vector of [C] 

NC(1)=NA(1)=NB(1)=M 
NC(2)=NA(2)=NB(2)=N 

REQUIRES: 

NA(1)=NB(1)=NC(1)=M 
NA(2)=NB(2)=NC(2)=N 

LIBRARY: /llb/1lbmatx.a 
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FUNCTION TRACE 


NAME: TRACE() 

DESCRIPTION: This function finds the trace of a matrix. 

CALL: t = trace(a,na,n) 

INPUT ARGUMENTS: 

A matrix packed by columns in one dimensional array 

NA dimension vector of [A] 

NA(l)=number of rows of [A] 

NA(2)=number of columns of [A] 

N number of rows and columns of [A] 

OUTPUT ARGUMENTS: 

TRACE trace of the input matrix 

REQUIRES: 

NA(1)=NA(2)=N 

LIBRARY: /IIb/1ibmatx.a 


SUBROUTINE TRANP 


NAME: TRANP() 

DESCRIPTION: This subroutine takes the transpose of a matrix; 


T 

[B]=[A] 


CALL: call tranp(a,na,b,nb,m,n) 


INPUT ARGUMENTS: 

matrix packed by columns In one dimensional array 
NA dimension vector of [A] 

NA(1)=number of rows of [A] 

NA(2)=number of columns of [A] 

M number of rows of [A] 

N number of columns of [A] 


OUT PITT ARGUMENTS: 

B 

NB 


T 

[B]=[A]; transposed matrix packed by columns in one 
dimensional array 
dimension vector of [B] 

NB(1)=N 
NB(2)=M 


REQUIRES: 

NA(1)=NB(2)=M 
NA(2)=NB(1)=N 
call to EQUATE 


LIBRARY: /lib/1ibmatx.a 
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SECTION VIII 
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APPENDIX A 

SCOLE INSTRUMENT CABINET PATCHING 




CABINET 1 


CABINET 2 


CABINET 3 


CMC 3: 


T-Motor 26-not connected 





T-Motor 27-not connected 





Tach 28 

-> 

neff 25 

-> 

DT-ADC 11 

Tach 29 

-> 

neff 26 

-> 

DT-ADC 12 

Resl 30 

-> 

neff 27 

-> 

DT-ADC 5 

Resl 31 

-> 

neff 28 

-> 

DT-ADC 6 

Res3 32 

“> 

neff 29 

-> 

DT-ADC 7 

Res3 33 

-> 

neff 30 

-> 

DT-ADC 8 



neff 31-not 

connected 




neff 32-not 

connected 


CMC 1: 





T-Motor 34-not connected 





T-Motor 35-not connected 





Tach 36 

-> 

neff 18 

-> 

DT-ADC 9 

Tach 37 

-> 

neff 19 

-> 

DT-ADC 10 

Resl 38 

"> 

neff 20 

-> 

DT-ADC 1 

Resl 39 

-> 

neff 21 

-> 

DT-ADC 2 

Res3 40 

-> 

neff 22 

-> 

DT-ADC 3 

Res3 41 

“> 

neff 23 

-> 

DT-ADC 4 



neff 24-not 

connected 


ACCELEROMETERS: 





1 (shaccx) 

-> 

neff 9 

-> 

ADC 7 

2 (shaccy) 

-> 

neff 10 

-> 

ADC 8 

3 (shaccz) 

-> 

neff 11 

-> 

ADC 9 

4 (mlaccx) 

-> 

neff 12 

-> 

ADC 10 

5 (mlaccy) 

-> 

neff 13 

-> 

ADC 11 

6 (m2accx) 

-> 

neff 14 

-> 

ADC 12 

7 (m2accy) 

-> 

neff 15 

-> 

ADC 13 

8 (refacx) 

-> 

neff 16 

-> 

ADC 14 

9 (refacy) 

-> 

neff 17 

-> 

ADC 15 

RATE TRANSDUCERS: 





Shuttle 10 (shurgr) 

-> 

neff 1 

-> 

ADC 1 

Shuttle 11 (shurgp) 

-> 

neff 2 

-> 

ADC 2 

Shuttle 12 (shurgy) 

-> 

neff 3 

-> 

ADC 3 

Mast 13 (refrgr) 

-> 

neff 4 

-> 

ADC 4 

Mast 14 (refrgp) 

-> 

neff 5 

-> 

ADC 5 

Mast 15 (refrgy) 

-> 

neff 6 

-> 

ADC 6 

PROOF INERTIA 





16 (mlwhx) 

“> 


-> 

DAC 4 

17 (mlwhy) 

-> 


-> 

DAC 5 

18 (m2whx) 

-> 


-> 

DAC 6 

19 (m2why) 

-> 


-> 

DAC 7 

20 (refwhx) 

-> 


-> 

DAC 1 

21 (refwhy) 

-> 


-> 

DAC 2 

22 (refwhz) 

-> 


-> 

DAC 3 


23- not connected 

24- not connected 
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25-not connected 


BUFFER 




1 

-> 

-> 

MT-8500 1 

2 

-> 

-> 

MT-8500 2 

3 

-> 

-> 

HT-8500 3 

4 

-> 

-> 

HT-8500 4 

5 

-> 

-> 

MT-8500 5 

6 

-> 

-> 

MT-8500 6 

7 

-> 

-> 

MT-8500 7 

8 

-> 

-> 

MT-8500 8 
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SUMMARY 

The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 
synthesis have been evaluated using computer simulation. In several cases, 
ground experiments have also been used to validate system performance under 
more realistic conditions. There remains a need, however, to test 

additional control laws for flexible spacecraft and to directly compare 
competing design techniques. In this paper an NASA program is discussed 
which has been initiated to make direct comparisons of control laws for, 
first, a mathematical problem, then an experimental test article is being 
assembled under the cognizance of the Spacecraft Control Branch at the NASA 
Langley Research Center with the advice and counsel of the IEEE Subcom¬ 
mittee on Large Space Structures. The physical apparatus will consist of a 
softly supported dynamic model of an antenna attached to the Shuttle by a 
flexible beam. The control objective will include the task of directing 
the line-of-sight of the Shuttle/antenna configuration toward a fixed 
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target, under conditions of noisy data, limited control authority and 
random disturbances. The open competition started in the early part 

of 1984. Interested researchers are provided information intended to 
facilitate the analysis and control synthesis tasks. A workshop is planned 
for early December at the NASA Langley Research Center to discuss and 

compare results. 

INTRODUCTION 

Many future spacecraft will be large and consequently quite flexible. 
As the size of antennae is increased, the frequencies of the first flex¬ 

ible modes will decrease and overlap the pointing system bandwidth. It 
will no longer be possible to use low gain systems with simple notch 
filters to provide the required control performance. Multiple sensors and 
actuators, and sophisticated control laws will be necessary to ensure 
stability, reliability and the pointing accuracy required for large, 

flexible spacecraft. 

Control of such spacecraft has been studied with regard given to 
modeling, order reduction, fault management, stability and dynamic system 
performance. Numerous example applications have been used to demonstrate 
specific approaches to pertinent control problems. Both computer simula¬ 
tions and laboratory experiment results have been offered as evidence of 
the validity of the approaches to control large, flexible spacecraft. 
Concerns remain, however, because of the chronic difficulties in control¬ 
ling these lightly damped large-scale systems. Because of these concerns 
and because of the desire to offer a means of comparing technical 
approaches directly, an NASA/IEEE Design Challenge is being offered. An 
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experimental test article is being assembled under the cognizance of the 
Spacecraft Control Branch at the NASA Langley Research Center with the 
advice and counsel of the IEEE (COLSS) Subcommittee on Large Space 
Structures. This Spacecraft Control Laboratory Experiment (SCOLE) will 
serve as the focus of a design challenge for the purpose of comparing 
directly different approaches to control synthesis, modeling, order 
reduction, state estimation and system Identification. 

The configuration of the SCOLE will represent a large antenna attached 
to the Space Shuttle orbiter by a flexible beam. This configuration was 
chosen because of its similarity to proposed space flight experiments and 
proposed space-based antenna systems. This paper will discuss the "Design 
Challenge” in terms of both a mathematical problem and a physical experi¬ 
mental apparatus. The SCOLE program is not part of any flight program. 


a 

A 

c 

d 

e 

E 

f 

n 

g 

GI 

I 


SYMBOLS 

acceleration vector ft/sec^ 
beam cross section area 
observation matrix 

noise contaminating direction cosine matrix measurements 

line-of-sight error 

modulus of elasticity 

concentrated force expressions 

force vector 

concentrated moment expressions 
torsional rigidity 

moment of inertia matrix for entire Shuttle/antenna configuration 
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Il 

14 

!♦ 

10 

If 

L 

Mi 

M4 

M D 

m 

m l 

104 

P 


V 1 

v 4 

U <J> 

u y 

X,Y,Z 

A 

5 


moment of inertia matrix. Shuttle body 

moment of inertia matrix, reflector body 

beam cross section moment of inertia, roll bending 

beam cross section moment of inertia, pitch bending 

beam polar moment of inertia, yaw torsion 

length of the reflector mast, beam 

control moment applied to the Shuttle body 

control moment applied to the reflector body 

disturbance moment applied to the Shuttle body 

mass of entire Shuttle/antenna configuration 

mass of Shuttle body 

mass of reflector body 

mass density of beam 

beam position variable 

direction cosine matrix, Shuttle body () . = T, ()„, , , 

9 J earth 1 Shuttle body 

direction cosine matrix, reflector body () - * T, () 

9 J earth 4 reflector 

body 

inertial velocity, Shuttle body 

inertial velocity, reflector body 

lateral deflection of beam bending in y-z plane 

lateral deflection of beam bending in x-z plane 

angular deflection of beam twisting about z axis 

position variables 

displacement of proof-mass actuator 

line-of-sight pointing requirement 

noise contaminating angular velocity measurements 







pitch, roll, heading 
C damping ratio 

x noise contaminating acceleration measurements 

0)^ angular velocity of Shuttle body 

u)^ angular velocity of reflector body 

DISCUSSION 

The objective of the NASA-IEEE Design Challenge concerning the control 
of flexible spacecraft is to promote direct comparison of different 
approaches to control, state estimation and systems identification* The 
design challenge has principal parts, the first using a mathematical model, 
and the second using laboratory experimental apparatus. The specific parts 
of the Spacecraft Control Laboratory Experiment (SCOLE) program will be 
discussed in detail. 


Control Objectives 

The primary control task is to rapidly slew or change the line-of- 
sight of an antenna attached to the space Shuttle orbiter, and to settle or 
damp the structural vibrations to the degree required for precise pointing 
of the antenna. The objective will be to minimize the time required to 
slew and settle, until the antenna line-of-sight remains within the 
angle <5. A secondary control task is to change direction during the 
"on-target" phase to prepare for the next slew maneuver. The objective Is 
to change attitude and stabilize as quickly as possible, while keeping the 
line-of-sight error less than <5. 
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Math Model Dynamics 


The initial phase of the design challenge will use a mathematical 
model of the Shuttle orbiter/antenna configuration. It is necessary to 
obtain a balance, of course, between complex formulations which might be 
more accurate and simplified formulations which ease the burden of 
analysis. 

The dynamics are described by a distributed parameter beam equation 
with rigid bodies, each having mass and inertia at either end. One body 
represents Space Shuttle orbiter; the other body is the antenna reflector. 
The equations for the structural dynamics and Shuttle motion are formed by 
adding to the rigid-body equations of motion, beam-bending and torsion 
equations. The boundary conditions at the ends of the beam contain the 
forces and moments of the rigid Shuttle and reflector bodies. The 
nonlinear kinetmatics couples the otherwise uncoupled beam equations. 
Additional terras represent the action of two, 2-axis proof-mass actuators 
at locations on the beam chosen by the designer. 

The rigid-body equations of motion for the Shuttle body are given by: 

“i - - ^'vri + M i + % ♦ Yi 5 
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Similarly, for the reflector body. 


“4 * - * M 4 + V> 


r 4 1 F B,4 
m. 


The direction cosine matrices defining the attitudes of the Shuttle and 
reflector bodies are given by: 

• T T 

T » - uj T 

i ii 


• T T 

T. - - tu T 
4 4 4 


The direction cosine matrices defining the attitudes of the Shuttle and the 
reflector bodies are related to the beam end conditions. 


1 0 0 

0 cosA<|> -sinA<t> 

0 sinA$ cosAiJi 


cosA0 0 sinAG 

0 1 0 

-sinAG 0 cosAG 


cosAY -sinAY 0 

sinAY cosAY 0 

0 0 1 


where: 


AY 


"t 


- U„ 


s=L 


A0 


A<J> 


3u 6 

3s 

3u, 

_ 4 

3s 


s=0 


3u e 

3s 


s=L 


s s L 


!j! 

3s 


s=0 


s*0 
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The equations of motion for the flexible beam-like truss connecting the 
reflector and Shuttle bodies consist of standard beam bending and torsion 
partial differential equations with energy dissapative terms which enable 
damped modes with constant characteristics for fixed, though dynamic, end 
conditions. The system of equations can be viewed as driven by changing 
end conditions and forces applied at the locations of the proof-mass 
actuators. 


ROLL BEAM BENDING: 


8 2 u, 


a 3 

3 u. 


PA ~ -T- + 2 ^d> •'PA El. — y- + El, —- l [f. 6(s-s ) + g. (s-s )] 

at 2 * * 3s 2 9t * 3s 4 n*l *’ n " *’ n 38 


PITCH BEAM BENDING: 


a 2 a 3 a 4 / 

9 9 9 u~ 4 

' ■ —- + EI 0 — - A, [ f e,„ S < s - s n> + *e,n U (s - s n>] 


PA —Y + 2C 0 /PA El 


3t 


3s 3t 


3s n=l 


YAW BEAM TORSION: 


a 2 a 3 a 2 

3 u», 3 Uy 3 u Y 

pi y I 2 + 2c 'i' I H' 7T77 + GI, f 


3t 


3s 3t 


3s n=l 


g »,n 5(s 


s n > 


where: 


.2 

3 u. 


m. 


*.i 1 at 2 


s=0 


(SHUTTLE BODY FORCE} 


3 2 «. 




♦» 2 2 3 1 2 


3 2 A 


+ m, 




2 3t 2 


8=S, 


{PROOF-MASS ACTUATOR FORCE} 
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{PROOF-MASS ACTUATOR} 


'<{>,3 


.2 

3 at 2 


3 2 A 


+ m 


hi 


3 2 

J 3t 


s*s. 


3 % 

f *’ 4 ' ra4 at 2 


a 2 

3 u,. 


.2 

3 u 

- I l —T~ /32 * 5 + F w 
zz, 4 at 2 y 


s=130 


m. 


0,1 1 Jt 2 


S^S, 


^2 

o u r 


m,. 


0 2 2 2 

u,z z 3t ^ 


+ *\j 

+ " 2 _ T^ 


s=s. 


a 2 

3 u,. 


0,3 


m 


3 2 
J 3t 


3 2 A 


+ in. 


hi 


3 2 

3t 


s=s. 


{REFLECTOR BODY FORCE} 


{SHUTTLE BODY FORCE} 


{PROOF-MASS ACTUATOR FORCE} 


{PROOF-MASS ACTUATOR FORCE} 


a 2 
9 u 


m, 


0 


0 * 4 4 at 2 


3 2 u„, 

- ^,4 ^ /18 * 75 - F x 


{REFLECTOR BODY FORCE} 


s=130 
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(SHUTTLE BODY, MOMENTS} 


/ g<M \ 

I g 0,l I = Vl + Vl“l + M 1 + *D 

' g Y,l / 



(PROOF-MASS ACTUATOR, MOMENT} 



(PROOF-MASS ACTUATOR, MOMENT} 



(REFLECTOR BODY, MOMENT} 


The angular velocity of the reflector body is related to the Shuttle body 
by: 
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The line-of-sight error described in figure 2 is affected by both the 
pointing error of the Shuttle body and the misalignment of the reflector 
due to the deflection of the beam supporting the reflector. The line-of- 
sight is defined by a ray from the feed which is reflected at the center of 
the reflector. Its direction in the Shuttle body coordinates is given by: 


r los 


h + "r + 2 

** 

• 

1 

S-/ 

V 

1 - 1 

CM 

1 

J- 

1 

E- 

r a r r " r f 

• r a 


where: 


Rp is the feed location (3.75, 0, 0) 

Rr is the location of the center of the reflector (18.75, -32.5, 
-130) 

R^ is a unit vector in the direction of the reflector axis in 
Shuttle body coordinates 


The vector R^ can be related to the direction cosine attitude matrices 
for the Shuttle body, Tj, and the reflector body, , by 



The relative alignment of the reflector to the Shuttle body is given by 
T 

T.T, which is a function of the structural deformations of the beam. 

1 4 
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The line-of-sight error, e, is the angular difference between the 
target direction, given by the unit vector, Dj, and the line-of-sight 
direction in Earth axes, TjRlqs* 

e - ARCSIS |d t X T^l or ARCSIN 

Computer programs are available which generate time histories of the 
rigid body and the mode shapes and frequencies for the body-beam-body 
configuration for "pitch" bending, "roll" bending and "yaw" twisting. 
Since the modes are based on solving explicitly the distributed parameter 
equations (without damping and without kinematic coupling) there is no 
limit to the number of modal characteristic sets that can be generated by 
the program. It will be the analyst's decision as to how many modes need 
to be considered. 


Laboratory Experiment Description 

The second part of the design challenge is to validate in the 

laboratory, the system performance of the more promising control system 
designs of the first part. The experimental apparatus will consist of a 
dynamic model of the Space Shuttle orbiter with a large antenna reflector 
attached by means of a flexible beam. The dynamic model will be exten¬ 
sively instrumented and will have attached force and moment generating 
devices for control and for disturbance generation. A single, flexible 
tether will be used to suspend the dynamic model, allowing complete angular 
freedom in yaw, and limited freedom in pitch and roll. An inverted 

position will be used to let the reflector mast to hang so that gravity 

effects on mast bending will be minimized. The dynamics of the laboratory 
model will of necessity be different from the mathematical model discussed 
earlier. 
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Design Challenge, Part One 


For part one of the design challenge, the following mathematical 
problem is addressed. Given the dynamic equations of the Shuttle/antenna 
configuration, what control policy minimizes the time to slew to a target 
and to stabilize so that the line-of-sight (LOS) error is held, for a time, 
within a specified amount, 6. During the time that the LOS error is 
within 6, the attitude must change 90° to prepare for the next slew 
maneuver. This was previously referred to as the sescondary control task. 
The maximum moment and force generating capability will be limited. Advan¬ 
tage may be taken of selecting the most suitable initial alignment of the 
Shuttle/antenna about its assigned initial RF axis, line-of-sight. 
Random, broad band-pass disturbances will be applied to the configuration. 
Two proof-mass, force actuators may be positioned anywhere along the beam. 
The design guidelines are summarized below: 

1. The initial line-of-sight error is 20 degrees. 

e(o) = 20 degrees 

2. The initial target direction is straight down. 



3. The initial alignment about the line-of-sight is free to be chosen 
by the designer. Advantage may be taken of the low value of 
moment of inertia in roll. The Shuttle/antenna is at rest 
initially. 

4. The objective is to point the line-of-sight of the antenna and 
stabilize to within 0.02 degree of the target as quickly as 
possible. 

6 = 0.02 degree 
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5. Control moments can be applied at 100 Hz sampling rate to both the 
Shuttle and reflector bodies of 10,000 ft-lb for each axis. The 

commanded moment for each axis is limited to 10,000 ft-lb. The 

actual control moment’s response to the commanded value is 
first-order with a time constant of 0.1 second. 

For the rolling moment applied to the Shuttle body: 


- 1 ° 4 i B X,l,=o™,andi *° 4 


- 0.1 


- 0 . 1 , 


Yl (n + » - e Vl <n) + (1 ' 6 > Yl.conunand^ 


Equations for other axes and for the reflector body are similar. 

6. Control forces can be applied at the center of the reflector in 
the X and Y directions only. The commanded force in a 
particular direction is limited to 800 lbs. The actual control 
force's response to the commanded value is first-order with a 
response time of 0.1 second. 

For the side for applied to the reflector body: 


-800 < F v . < 800 

— Y,command — 


F Y (n + 1) - r 0 - 1 F Y (n) + (1 - a" 0 ’ 1 ) F Y , command <"> 


Equations for X-axis are similar. 

7. Control forces using two proof-mass actuators (each having both 
X and Y axes) can be applied at two points on the beam. The 
strokes are limited to ± 1 ft, and the masses weight 10 lbs each. 
The actual stroke follows a first-order response to limited 
commanded values. 
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For the X-axis of the proof-mass actuator at S 2 ‘- 


* — A X, 2,command — * 


A X,2 (n + 1} " A X,2 (n) + ° " e " 0,1) A X,2,command^ n) 


Equations for other axes and locations are similar. 

8. The inertial attitude direciton cosine matrix for the Shuttle body 
lags in time the actual values by 0.01 second and are made at a 


noise has zero mean. 


T .(n + 1) * T „ (n) + 

s,measured s,true 


where: 


e{ dij <n>1 * 0 

E{d 1J (n)d kL (n)} - 0 

E{dij(n)djj (n + k)} = 0 


Each 

element of 

: the direction 

s contaminated 

by additive. 

j an rms value 

of 0.001. The 

t u < n > 

d 12 (n> 

d l 3 (n) 


'jiU) 

d 22 Cn> 

d 2 3 (n) 


! 31 (n) 

d 32<"> 

d 3 3 (n) 


for i 

* k or j 

* L 



for k * 0 


[. 001 ] 


for k = 0 


706 




9. The angular velocity measurements for both the Shuttle and 
reflector bodies pass through a first-order filter with 0.05 sec 
time constant and lag in time the actual values by 0.01 second and 
are made at a rate of 100 samples per second. Each rate 

measurement is contaminated by additive, Gaussian, uncorrelated 
noise having an rms value of 0.02 degree per second. The noise 
has zero mean. 

For example: 


(0 


1 


(n + 1) a (Jj 

,X,measured^ 1,X,filtered 

E ^ G i x (n) e i x (n + = 0 

- (. 02) 2 


(n) + c 1 >x (n) 
for k * 0 

for k = 0 


where 


1,X,filtered 


20 U 1,X,filtered + 2( ^ ^l.X.true 


10. Three-axis accelerometers are located on the Shuttle body at the 
base of the mast and on the reflector body at its center. Two- 
axes (X and Y) accelerometers are located at intervals of 
10 feet along the mast. The acceleration measurements pass 
through a first-order filter with a 0.05 second time constant and 
lag in time the actual values by 0.01 second, and are made at a 
rate of 100 samples per second. Each measurement is contaminated 
by Gaussian additive, uncorrelated noise having an rms value of 
0.05 ft/sec^. 
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For example: 


3 1 ,X, measured^ 11 + ^ 3 1,X,filtered^ 1 ^ + T 1 ,X^ n ^ 

e{t (n) t (n + k)} = 0 for k * 0 

1 y A 1 )A 

* (.05)2 for k = 0 

where: 

a = - 20 3 + 20 w 

d l,X,filtered u ,X,filtered l,X,true 

11. Gaussian, uncorrelated step-like disturbances are applied 
100 times per second to the Shuttle body in the form of 3-axes 
moments, having rms values of 100 ft-lbs. These disturbances 
have zero mean. 

For example: 

E {m d x (n) M d x (n + k)} = 0 for k * 0 

= (100) 2 for k - 0 

In summary, the designer's task for part one is to: (1) derive a 
control law for slewing and stabilization, coded in FORTRAN; (2) select an 
initial attitude in preparation for slewing 20 degrees; and (3) select two 
positions for the 2-axes proof-mass actuators. An official system 
performance assessment computer program will be used to establish the time 
required to slew and stabilize the Shuttle/antenna configuration. 
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Design Challenge, Part Two 


As in part one, the task is to minimize the time to slew and stabilize 
a Shuttle/antenna configuration. The difference is that in part two of the 
design challenge, a physical laboratory model will be used instead of the 
dynamic equations of part one. The constraints on total moment and force 
generation capability will apply to part two, as for part one. Again, the 
analyst may select the initial alignment about the assigned initial RF 
line-of-sight. Disturbances will be injected into the Shuttle/antenna 
model. The designer’s task will be similar to that for part one. 

CONCLUDING REMARKS 

A Design Challenge, in two parts, has been offered for the purpose of 
comparing directly different approach to controlling a flexible 
Shuttle/antenna configuration. The first part of the design challenge uses 
only mathematical equations of the vehicle dynamics; the second part uses a 
physical laboratory model of the same configuration. The Spacecraft 
Control Laboratory Experiment (SCOLE) program is being conducted under the 
cognizance of the Spacecraft Control Branch at the NASA Langley Research 
Center. The NASA/IEEE Design Challenge has the advice and counsel of the 
IEEE“COLSS Subcommittee on Large Space Structures. Workshops will be held 
to enable investigators to compare results of their research. 
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MASS CHARACTERISTICS 
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The moment of inertia becomes: 
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Figure 1. Drawing of the Shuttle/Antenna Configuration. 
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Roll bending mode # 1. Frequency= .32 Hz 



Roll bending mode § 2. Frequency= 1.29 Hz 



Roll bending mode # 3. Frequency= 4.80 Hz 



Roll bending mode # 4. Frequency= 12.29 Hz 



Roll bending mode # 5. Frequency= 23.68 Hz 



Roll bending mode # 6. Frequency= 38.89 Hz 



Roll bending mode # 7. Frequency^ 57.90 Hz 



Roll bending mode § 8. Frequency= 80.72 Hz 


Figure 4«i.- Plots of normalized roll bending mode shapes 
for SCOLE configuration. 

7 / 5 * 


INTOrU fl. 1NF0PLT 7. ItTCTU 6- l^OTLT 6. tPTOPU 4. IffCPLT 3. MFOPLT 2- WOPLT 




Pitch bending mode # 1. Frequency= .29 Hz 



Pitch bending mode # 2. Frequency^ 1.65 Hz 



Pitch bending mode # 3. Frequency= 4.97 Hz 



Pitch bending mode # 4 Frequency= 12.36 Hz 



Pitch bending mode # 5. Frequency= 23.72 Hz 



Pitch bending mode # 8. Frequency= 80.73 Hz 


Figure 4b,- Plots of normalized pitch bending mode shapes 
for SCOLE configuration. 
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Figure 4c.- Plots of normalized torsional mode shapes for 
SCOLE configuration. 
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IN THE PRESENTATION ... 

* INTRODUCTION TO THE SCOLE EXPERIMENT 

* LATEST SYSTEM UPGRADE AND CURRENT STATUS 

OF THE SCOLE FACILITY 

* RESEARCH ACTIVITY (IN HOUSE & GI PROGRAM) 

* SOME EXPERIMENT RESULTS 

Movie: The Optimal Attitude Control Testing on SCOLE 

* FUTURE PLANS 

* SUMMARY 
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INTRODUCTION TO THE SCOLE EXPERIMENT 


* THE NASA-IEEE DESIGN CHALLENGE (Taylor-Balakrishnan, 1984) 

Control of flexible spacecraft to promote direct comparison of 
different approaches to: 

— Control 

— State Estimation 

— System Identification 
Using: 1. A mathematical model 

2. A laboratory experimental apparatus 

* CONTROL OBJECTIVES 

— Rapid slew or change of line—of—sight of an antenna attached 
to the space Shuttle orbiter 

— Damp the structural vibrations to a prescribed degree 

— Minimize the time for slew and settle to a given orientation 

— Change direction during the "on—target" phase to prepare for 
the next slew maneuver 
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INTRODUCTION (CONT.) 


* THE SPACECRAFT CONTROL LABORATORY EXPERIMENT - 

SCOLE 

— Experiment conceived to provide a common design challenge 
for interested investigators 

— Develop and validate control methods for large flexible 
spacecraft configuration 

— Demonstration of slewing and pointing control problems 
along with vibration suppression 

— Utilizing antenna—like structural configuration and inertial 
sensors and actuators 

— Variety of inertial quality sensor and actuator types 

+ Accelerometers, rate sensors, optical position sensors 
+ Thrusters, control moment gyros, reaction wheels 

— Real time computing capability 
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LATEST SYSTEM UPGRADE AND CURRENT STATUS 

OF THE SCOLE FACILITY 

LATEST UPGRADE: 

* DUAL CAMERA ATTITUDE MEASUREMENT SYSTEM 

* CAMERA POSITION IDENTIFICATION AND SCOLE 

ATTITUDE TRACKING ALGORITHMS 

* POWER AMPLIFIERS TO ACCOMODATE THE HI-TORQUE 

REACTION WHEELS 

* REAL TIME COMPUTER UPGRADE - HARDWARE AND 

SOFTWARE 

* SCOLE DOCUMENTATION - HARDWARE AND SOFTWARE 

* DEVELOPMENT OF CMGs STEERING LAW - IN PROGRESS 
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LATEST SYSTEM UPGRADE AND CURRENT STATUS 
OF THE SCOLE FACILITY (CONT.) 

UNDER THE CURRENT STATUS, SCOLE IS CAPABLE OF 
CONDUCTING THE FOLLOWING TESTS: 

* GENERATE VIBRATION DATA FILES FOR DISTRIBUTION 

* ON-LINE PARAMETER IDENTIFICATION 

* VIBRATION SUPPRESSION - CANTILEVER MODE TESTING 

* FIXED PIVOT POINT - 3 DOF TESTING 

* FREE PIVOT POINT - 5 DOF TESTING 
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RESEARCH ACTIVITY (IN-HOUSE AND GI PROGRAM) 

* NONLINEAR CONTROL DESIGN AND RELIABILITY ISSUES 

* DEVELOP AND TEST A DISTRIBUTED ADAPTIVE CONTROL 

SYSTEM 

* DECENTRALIZED CLOSED LOOP MANEUVER 

* DISTRIBUTED OUTPUT MODEL-FOLLOWING CONTROL LAW 

* PARAMETER IDENTIFICATION AND PERTURBATION METHOD 

CONTROLLER 

* ANALYTIC REDUNDANCY MANAGEMENT FOR SYSTEM WITH 

APPRECIABLE STRUCTURAL DYNAMICS 

* RIGID BODY ATTITUDE TRACKING ALGORITHMS 

* CONTROL MOMENT GYRO STEERING LAW 
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SOME EXPERIMENTAL RESULTS 

* KALMAN FILTER ESTIMATE 

* PARAMETER IDENTIFICATION ON THE SCOLE MAST 

* ON-LINE IDENTIFICATION AND ATTITUDE CONTROL 

FOR SCOLE 

* MOVIE: SCOLE ATTITUDE CONTROL 
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KALMAN FILTER ESTIMATE 
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THE REAL TIME SCOLE COMPUTER 


* CHARLES RIVER DATA SYSTEM (CRDS) - UNIVERS 32 

MICRO-COMPUTER BASED ON THE MOTOROLA M68020 
MICROPROCESSOR AND M68881 FLOATING POINT 
CO-PROCESSOR 

* RANDOM ACCESS MEMORY 4.0 M-BYTE 

* TWO 80 M-BYTE HARD DISKS 

* OPERATING SYSTEM IS BASED ON AT&T UNIX SYSTEM V 

* COMPILERS - FORTRAN, C, PASCAL 
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CRDS HARDWARE 
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CRDS SCOLE SOFTWARE 


* MENU DRIVEN SCOLE EXPERIMENT OPERATING ENVIRONMENT 

* USER-FRIENDLY FORTRAN k C SUBROUTINES AND FUNCTIONS 

TO CONTROL EXPERIMENT HARDWARE 

* CAPABILITY TO STORE REALTIME SENSOR k ACTUATOR DATA 

IN MATRIX-X, MATLAB AND OTHER POST PROCESSING 
SOFTWARE PACKAGES 
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CRDS OPERATING ENVIRONMENT 

* TCP/IP NETWORKING CAPABILITY 

* USER CHOICE TO OPERATE IN UNIX V BOURNE SHELL 

OR UNOS VIII OPERATING SYSTEM 

* FORTRAN, C, AND PASCAL COMPILERS AVAILABLE 

* GKS GRAFPAK GRAPHICS 


* SYSTEM ADMINISTRATOR HELP 
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FUTURE PLANS 


* DEVELOPMENT OF CMGs STEERING LAW TO ACCOMODATE: 

— Maneuver bounds: yaw +/- 45 degrees 

pitch -f/— 20 degrees 

roll +/— 20 degrees 

_ Generation of CMG gimbal rate command that results 

in the desired control torque 

* BENCHMARK PROGRAMS THAT ILLUMINATE: 

— Computer architectural issues that affect 
the speed of software 

— Hardware performance of the SCOLE system 

— Control effectiveness parameters 

* REAL TIME COMPUTER UPGRADE (HARDWARE-SOFTWARE) 

ACCORDING TO FUTURE NEEDS 

* DOCUMENTATION 

— Complete software user manual 

— Experimental setup diagrams 

— Various wiring diagrams to accomodate the in—house user 
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SUMMARY 


* THE SCOLE LABORATORY FACILITY IS PROVIDED AS A 
TEST-BED FOR EVALUATION OF CONTROL LAWS FOR 
LARGE FLEXIBLE STRUCTURES 

* THE SCOLE APPARATUS IS EQUIPPED WITH SENSORS, 
ACTUATORS AND A REAL TIME COMPUTER THAT MAY 
ACCOMODATE VARIOUS CONTROL LAWS AND TESTS IN 
THE RESEARCH OF LARGE FLEXIBLE STRUCTURES 

* SEVERAL CONTROL LAWS TESTED ON THE SCOLE SHOWED 
GOOD AGREEMENT WITH THEORY AND COMPUTER 
SIMULATIONS 

* GUEST INVESTIGATORS ARE INVITED NOT ONLY TO 
DERIVE CONTROL LAWS BUT ALSO TO TEST THEM 
IN-HOUSE ON THE SCOLE FACILITY 
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